( \sqrt{\pi \left ( 8x^3+4x^2+x+1/30 \right )^{1/3}}\left ( \frac{x}{e} \right ) ^x )/gamma(x+1)