plot([(e*sqrt(x))*(x/e)^x,x!],(x,0,5))