double fact(int n) { int i; double product; product = 1; for(i = 1; i<=n; i++) product=product*i; return product; }