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?