factorial function

double factorial (double n)
{
    double x = n;
    for(int i = n; i > 2; i--)
        {
            x *= (i-1);
        } // end for
    return x;
} // end function