Na verdade foi erro meu, copiei o mesmo valor.
Mathematica:
Gamma[1.5, 0, 0.10] = 0.019861
R:
> pgamma(1.5, .10, lower=FALSE)*gamma(.10)
[1] 0.1079214
Vendo a formula achei que o certo seria: pgamma(1.5, .10, lower=FALSE)*gamma(1.5) mesmo assim não é a mesma resposta do mathematica.
[ ],
Pedro Rafael Diniz Marinho.
---------- Forwarded message ----------
From:
beniltoncarvalho [via R-br] <ml-node+s2285057n4656817h4@n4.nabble.com>
Date: 2012/11/1
Subject: Re: [R-br] [Dúvida] Função Gama Incompleta.
To: Pedro Rafael Diniz Marinho <
pedro.rafael.marinho@gmail.com>
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?