[ReiniBugman]

double fact( unsigned long n )
{
double retval;
retval=1.0;
while (n>0)
retval=retval*n--;
return retval;
}