计算浮点数的阶乘

因为学习的原因需要计算浮点数的阶乘，通常说的阶乘都是相对于整数来说的。

通过搜索相关的资料，得知可以通过推广 Gamma公式到实数范围，得到求伽玛函数的一种数值解。代码如下：

Code
 1         /// <summary>
 2         /// 近似计算伽玛函数。
 3         /// </summary>
 4         /// <param name="xx">待计算阶乘的数</param>
 5         /// <returns>计算结果</returns>
 6         static double Gamma(double xx)
 7         {
 8             double x, y, tmp, ser;
 9             double[] cof=new double[] {76.18009172947146,-86.50532032941677,24.01409824083091,-1.231739572450155,0.1208650973866179e-2,-0.5395239384953e-5};
10             int j;
11             x = xx;
12             y = x;
13             tmp = x + 5.5;
14             tmp -= (x + 0.5) * Math.Log(tmp);
15             ser = 1.000000000190015;
16             for (j = 0; j <= 5; j++)
17             {
18                 ser += cof[j] / (++y);
19             }
20             return xx*Math.Exp((-tmp + Math.Log(2.5066282746310005 * ser / x)));
21 

参考资料：

• Calculators and the Gamma Function http://www.rskey.org/gamma.htm
• 《C数值计算》 电子工业出版社 第二版 P156-157

转载于:https://www.cnblogs.com/hawking106123/archive/2009/09/26/1574604.html