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