Plot [(√((2*x+1/3)*π)*x^x/e^x - Γ(x+1))/Γ(x+1)], x=1..60