不完全伽马函数的计算

VC中实现的计算不完全伽马函数的数值计算程序。

//   基本思想是用切比雪夫多项式来逼近
#include “stdafx.h“
#include “gamma1.h“        
#include “stdio.h“

double gamma1（ double x）
{
	int i;
	double ytsu;
	static double a[11]={0.0000677106-0.0003442342
						0.0015397681-0.00244674800.0109736958
						-0.00021090750.07423790710.0815782188
						0.41184025180.42278433701.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）;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3741  2010-07-08 21:58  Gamma_Incomplete\Debug\gamma1.obj

     文件       4910  2010-07-08 22:03  Gamma_Incomplete\Debug\gamma2.obj

     文件     192585  2010-07-08 22:03  Gamma_Incomplete\Debug\Gamma_Incomplete.exe

     文件     204116  2010-07-08 22:03  Gamma_Incomplete\Debug\Gamma_Incomplete.ilk

     文件       3039  2010-07-08 22:03  Gamma_Incomplete\Debug\Gamma_Incomplete.obj

     文件     187128  2010-07-08 21:58  Gamma_Incomplete\Debug\Gamma_Incomplete.pch

     文件     443392  2010-07-08 22:03  Gamma_Incomplete\Debug\Gamma_Incomplete.pdb

     文件       1853  2010-07-08 21:58  Gamma_Incomplete\Debug\StdAfx.obj

     文件      41984  2010-07-22 11:01  Gamma_Incomplete\Debug\vc60.idb

     文件      53248  2010-07-08 22:03  Gamma_Incomplete\Debug\vc60.pdb

     文件        720  2010-07-08 19:57  Gamma_Incomplete\gamma1.cpp

     文件         42  2010-07-22 11:01  Gamma_Incomplete\gamma1.h

     文件       1086  2010-07-08 21:59  Gamma_Incomplete\gamma2.cpp

     文件         60  2010-07-08 21:39  Gamma_Incomplete\gamma2.h

     文件        505  2010-07-08 22:03  Gamma_Incomplete\Gamma_Incomplete.cpp

     文件       4904  2010-07-08 22:33  Gamma_Incomplete\Gamma_Incomplete.dsp

     文件        557  2010-07-08 21:34  Gamma_Incomplete\Gamma_Incomplete.dsw

     文件      50176  2010-07-22 11:01  Gamma_Incomplete\Gamma_Incomplete.ncb

     文件      49664  2010-07-22 11:01  Gamma_Incomplete\Gamma_Incomplete.opt

     文件        266  2010-07-22 11:01  Gamma_Incomplete\Gamma_Incomplete.plg

     文件       1268  2010-07-08 21:34  Gamma_Incomplete\ReadMe.txt

     文件        303  2010-07-08 21:34  Gamma_Incomplete\StdAfx.cpp

     文件        667  2010-07-08 21:34  Gamma_Incomplete\StdAfx.h

     目录          0  2010-07-08 22:03  Gamma_Incomplete\Debug

     目录          0  2010-07-22 11:01  Gamma_Incomplete

----------- ---------  ---------- -----  ----

              1246214                    25