 Solution

With the proton fixed, $[-{\hslash }^2{\nabla }^2/(2m_{\mu })-e^2/(4\pi {\epsilon }_{0}r)]\psi =E\psi$: kinetic energy plus Coulomb attraction, assuming nonrelativistic point particles and neglecting spin, radiation, and proton motion (or replacing $m_{\mu }$ by the reduced mass). Substitution of the $2s$ form gives $E=-{\hslash }^2{\alpha }^2/(8m_{\mu })$ and $\alpha =m_{\mu }{e}^2/(4\pi {\epsilon }_0{\hslash }^2)$. Direct differentiation verifies it. The curve starts positive, crosses at $r=2/\alpha$, then approaches zero from below; ordinary hydrogen has the same shape but the smaller electron mass gives a larger radial scale.