Acho q nao deu resultado diferente nao... ou eu estou enxergando coisas estranhas... vc postou no seu email original:
No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921
e eu postei:
pgamma(1.5, .10, lower=FALSE)*gamma(.10)
que vc confirmou como:
pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214
acho q nao consigo ver a diferenca?