double factorial ( int n )
{
  int i;
  double value;

  value = 1.0;
  for ( i = 1; i <= n; i++ )
  {
    value = value * i;
  }
  return value;
}