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

代码语言

代码分类

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

作者:Dezai.CN / 发布于2013/10/22/ 520


	简单的说就是整数阶乘的推广，它有一个积分的表达式：
	&Gamma;(x)=&int;e^(-t)*t^(x-1)dt (积分的下限式0,上限式+&infin;)

	#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));
	}

评论列表

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

试试其它关键字

　伽马函数　

同语言下

可能有用的

Dezai.CN贡献的其它代码(4037)