f[1]=1

  for i=2 to n do

    f[i]=f[i-1]*i