伽马函数计算(可求小数)代码,C/C++代码库,德仔网

代码语言

代码分类

【C/C++】伽马函数计算(可求小数)

作者:睿杰 / 发布于2016/5/11/ 517

简单的说就是整数阶乘的推广，它有一个积分的表达式： Γ(x)=∫e^(-t)*t^(x-1)dt (积分的下限式0,上限式+∞) 算法源自《常用算法程序集》徐士良

#include "stdio.h"

double Gamma(x)
double x;
{
    int i;
    double y,t,s,u;
    static double a[11]={ 0.0000677106,-0.0003442342, 0.0015397681,
                         -0.0024467480, 0.0109736958,-0.0002109075,
                          0.0742379071, 0.0815782188, 0.4118402518,
                          0.4227843370, 1.0};
    if (x<=0.0)
    {
        printf("err**x<=0!\n");
        return(-1.0);
    }

    y=x;
    if (y<=1.0)
    {
        t=1.0/(y*(y+1.0));
        y=y+2.0;
    }
    else
        if(y<=2.0)
        {
            t=1.0/y;
            y=y+1.0;
        }
        else
            if(y<=3.0)
            {
                t=1.0;
            }
            else
            {
                t=1.0;
                while (y>3.0)
                {
                    y=y-1.0;
                    t=t*y;
                }
            }
    s=a[0];
    u=y-2.0;
    for (i=1; i<=10; i++)
    {
        s=s*u+a[i];
    }
    s=s*t;
    return(s);
}

double DecimalFactorial(double x)
{
    return Gamma(x+1);
}

void main()
{
    printf("%f\n",DecimalFactorial(3.5));
}

评论列表

本站所提供的代码，版权归原作者所有，若有侵犯作者版权，请与我们联系，我们将立即删除或修改。谢谢!
本站所有代码发布及提供者。

试试其它关键字

　伽马　

同语言下

可能有用的

睿杰贡献的其它代码(15)