fact(n)=my(p=1);for(k=2,n,p*=k);p