資源簡介
一個關于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=“<?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
評論
共有 條評論