Gamma函数的解法

一个关于Gamma函数的数值解法，很实用，准确度很高1

#include
#include

using namespace std;

double gammln（const xx）
{
	int j;
	double xytmpser;
	static const double cof[6]={76.18009172947146
	-86.5053203294167724.0140982408391-1.231739572450155
	0.1208650713866179e-2-0.5395239384953e-5};

	y=x=xx;
	tmp=x+5.5;
	tmp-=（x+0.5）*log（tmp）;
	ser=1.000000000190015;
	for（j=0;j<6;j++） ser+=cof[j]/++y;
	return -tmp+log（2.5066282746310005*ser/x）;

}

double main（）
{
	int c=0.25;
	double xx=6y;
	y=exp（gammln（xx））;
	cout<<“y=“<< y< return 0;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        535  2009-02-20 10:26  Gamma\Gamma.dsw

     文件      41984  2009-02-20 21:32  Gamma\Gamma.ncb

     文件      74752  2009-02-20 21:32  Gamma\Debug\vc60.idb

     文件     110592  2009-02-20 21:32  Gamma\Debug\vc60.pdb

     文件    2009712  2009-02-20 12:43  Gamma\Debug\Gamma.pch

     文件     557102  2009-02-20 21:32  Gamma\Debug\Gamma.exe

     文件    1090560  2009-02-20 21:32  Gamma\Debug\Gamma.pdb

     文件     148390  2009-02-20 21:32  Gamma\Debug\Gamma.obj

     文件     792400  2009-02-20 21:32  Gamma\Debug\Gamma.ilk

     文件       1232  2009-02-20 21:32  Gamma\Gamma.plg

     文件       4632  2009-02-20 10:50  Gamma\Gamma.dsp

     文件        546  2009-02-20 21:32  Gamma\Gamma.cpp

     文件      48640  2009-02-20 21:32  Gamma\Gamma.opt

     目录          0  2009-02-20 10:26  Gamma\Debug

     目录          0  2009-02-20 10:26  Gamma

----------- ---------  ---------- -----  ----

              4881077                    15