C#版本计算Gamma函数

清龙

1.     static class Calc
   
2.     {
   
3.         static public double NullD = 1e300;
   
4.         static public double GammaFun(double x)
   
5.         {
   
6.             if (x == Math.Floor(x))
   
7.             {
   
8.                 if (x > 0)
   
9.                 {
   
10.                     int i = 1;
    
11.                     for (int j = 2; j < x; j++)
    
12.                         i *= j;
    
13.                     return Convert.ToDouble(i);
    
14.                 }
    
15.             }
    
16.             else
    
17.             {
    
18.                 double z = 0, r = 0, m = 0;
    
19.                 if (Math.Abs(x) > 1)
    
20.                 {
    
21.                     z = Math.Abs(x);
    
22.                     m = Math.Floor(x);
    
23.                     r = 1.0;
    
24.                     for (int k = 1; k <= m; k++)
    
25.                         r = r * (z - k);
    
26.                     z = z - m;
    
27.                 }
    
28.                 else
    
29.                     z = x;
    
30. 
    
31.                 double[] G = new double[26]{
    
32. 1.0,0.5772156649015329,-0.6558780715202538,
    
33. -0.420026350340952e-1,0.1665386113822915e0,-0.421977345555443e-1,
    
34. -0.96219715278770e-2, 0.72189432466630e-2,
    
35. -0.11651675918591e-2, -0.2152416741149e-3,
    
36. 0.1280502823882e-3, -0.201348547807e-4,
    
37. -0.12504934821e-5, 0.11330272320e-5,
    
38. -0.2056338417e-6, 0.61160950e-8,
    
39. 0.50020075e-8, -0.11812746e-8,
    
40. 0.1043427e-9, 0.77823e-11,
    
41. -0.36968e-11, 0.51e-12,
    
42. -0.206e-13, -0.54e-14,
    
43. 0.14e-14, 0.1e-15};
    
44.                 double Gr = G[25];
    
45.                 for (int k = 24; k >= 0; k--)
    
46.                     Gr = Gr * z + G[k];
    
47.                 double Ga = 1.0 / (Gr * z);
    
48.                 if (Math.Abs(x) > 1.0)
    
49.                 {
    
50.                     Ga = Ga * r;
    
51.                     if (x < 0)
    
52.                         Ga = -1 * Math.PI / (x * Ga * Math.Sin(Math.PI * x));
    
53.                 }
    
54.                 return Ga;
    
55.             }
    
56.             return NullD;
    
57.         }
    
58. 
    
59.     }
    
60. 
    

    static class Calc<br />     {<br />         static public double NullD = 1e300;<br />         static public double GammaFun(double x)<br />         {<br />             if (x == Math.Floor(x))<br />             {<br />                 if (x > 0)<br />                 {<br />                     int i = 1;<br />                     for (int j = 2; j < x; j++)<br />                         i *= j;<br />                     return Convert.ToDouble(i);<br />                 }<br />             }<br />             else<br />             {<br />                 double z = 0, r = 0, m = 0;<br />                 if (Math.Abs(x) > 1)<br />                 {<br />                     z = Math.Abs(x);<br />                     m = Math.Floor(x);<br />                     r = 1.0;<br />                     for (int k = 1; k <= m; k++)<br />                         r = r * (z - k);<br />                     z = z - m;<br />                 }<br />                 else<br />                     z = x;<br /> <br />                 double[] G = new double[26]{<br /> 1.0,0.5772156649015329,-0.6558780715202538,<br /> -0.420026350340952e-1,0.1665386113822915e0,-0.421977345555443e-1,<br /> -0.96219715278770e-2, 0.72189432466630e-2,<br /> -0.11651675918591e-2, -0.2152416741149e-3,<br /> 0.1280502823882e-3, -0.201348547807e-4,<br /> -0.12504934821e-5, 0.11330272320e-5,<br /> -0.2056338417e-6, 0.61160950e-8,<br /> 0.50020075e-8, -0.11812746e-8,<br /> 0.1043427e-9, 0.77823e-11,<br /> -0.36968e-11, 0.51e-12,<br /> -0.206e-13, -0.54e-14, <br /> 0.14e-14, 0.1e-15};<br />                 double Gr = G[25];<br />                 for (int k = 24; k >= 0; k--)<br />                     Gr = Gr * z + G[k];<br />                 double Ga = 1.0 / (Gr * z);<br />                 if (Math.Abs(x) > 1.0)<br />                 {<br />                     Ga = Ga * r;<br />                     if (x < 0)<br />                         Ga = -1 * Math.PI / (x * Ga * Math.Sin(Math.PI * x));<br />                 }<br />                 return Ga;<br />             }<br />             return NullD;<br />         }<br /> <br />     }<br /> <br />

直接调用Calc.GammaFun(x)就出来结果。用的是静态类的写法，像Math类一样 Calc.NullD 表示double类型的缺测数据

xiaolu001

感谢分享，学习。

鱼香Jack

感谢楼主分享，学习了

y18735399

有附带简单说明吗  看别人编的纯代码真是件头大的事

qxtlyf

要加以说明或介绍

lai85211

看看   来  下  额

清龙

时间久远，我都不记得我写了这个了。