Discussing Stirling's approximation, Wolfram Mathworld article mentions a modification of it due to Gosper: instead of the usual

$$n!\approx n^ne^{-n}\sqrt{2n\pi},$$

we have a tiny addition of $\frac13$ to $2n$ under the radical sign:

$$n!\approx n^ne^{-n}\sqrt{\left(2n+\frac13\right)\pi}.$$

This radically improves the precision of approximation for all $n\ge0$. But all the places I've seen discussing it don't explain how it was obtained. Wolfram Mathworld just says about Gosper's modification

...a better approximation to n! (i.e., one which approximates the terms in Stirling's series instead of truncating them)...

but it doesn't make me understand how this change was derived.

So, how to derive this modification? And is there a further improvement to the whole Stirling series, or is it just "whole series stuffed into one expression"?