factln.cpp

此程序为VC++常用数值算法这本书附赠的光盘中包含了本书中全部的源代码

double factln(int n)
{
	int i;
    double a[101],temp;
    for( i = 1;i<=100;i++)
        a[i] = -1.0;
    if (n < 0)
	{
        cout<<"negative factorial";
        exit(1);
	}
    if (n <= 99)
	{
        if( a[n + 1] < 0.0)  a[n + 1] = gammln(n + 1.0);
        temp = a[n + 1];
	}
    else
        temp = gammln(n + 1.0);
	return temp;
}