plot abs(sqrt(2pi*n)*(n/e)^n*(1+1/(12*n))-n!)/n! with n from 1 to 1000