怎么计算伽玛函数的近似值？

Ysu2008

它那无穷上限怎么数值积分？

ccmmjj

记不太清楚了，是用分部积分法和递推法吧。对几类特殊的阿法值有很漂亮的结果，但对一般的值通常都是查表或用某种数学软件。

Ysu2008

顺便问一下，上哪儿去找高精度的标准正态分布表或者是值？

我想自己计算高精度正态分布表。我现在是利用龙贝格公式计算上下限有限的数值积分，比如计算：

如果控制精度设定为 1E-20 ，则算得圆周率为 3.141592653589793 ，精度似乎还不错。就是不知算正态分布怎么样，想要一张高精度的表比对一下。

Ysu2008

ccmmjj 发表于 2016-6-15 17:05
记不太清楚了，是用分部积分法和递推法吧。对几类特殊的阿法值有很漂亮的结果，但对一般的值通常都是查表或 ...

我就是想知道数学软件是怎么算的，我想自己算。

ccmmjj

我的年代是学查表的。数学软件的使用不是很清楚。这个问题要问elim.他使用软件计算能力挺强。
至于计算软件怎么计算，我以前听说是用级数展开来算的。你是在编计算程序吗？
陆老师是概率统计教授，高精分布表请问他有没有。你需要几位小数的？

elim

大概搜索了一下：

luyucheng1

excel中有很多内置函数。我在做皮尔逊三数值积分计算时，伽马函数就是调用 GAMALN 计算的。有的你可以在其中找找看。

永远

顶一下……～～

coolboy

FUNCTION gammln(xx)
REAL gammln,xx
INTEGER j
DOUBLE PRECISION ser,stp,tmp,x,y,cof(6)
SAVE cof,stp
DATA cof,stp/76.18009172947146d0,-86.50532032941677d0,
& 24.01409824083091d0,-1.231739572450155d0,.1208650973866179d-2,
& -.5395239384953d-5,2.5066282746310005d0/
x=xx
y=x
tmp=x+5.5d0
tmp=(x+0.5d0)*log(tmp)-tmp
ser=1.000000000190015d0
do j=1,6
y=y+1.d0
ser=ser+cof(j)/y
enddo
gammln=tmp+log(stp*ser/x)
return
END