The planar spring-pendulum is modeled by the set of equations $$m\ddot{r}=mr\dotθ^2 + mg\cos\theta − k(r − L)$$ and, $$r^2\ddot{\theta}= −2r\dot r\dotθ − gr\sin\theta$$ First we should define the angular momentum by $p_\theta = mr^2\dot\theta$ and the radial momentum by $p_r = m\dot r$ and rewrite the spring-pendulum system as a set of four first-order ODEs for $x = (r, \theta , p_r, p_\theta )$. So I did like this, but I don't know if it's correct: $$\dot r = p_r/m, ~\dot\theta=p_\theta/mr^2$$ and, $$\dot p_r = mr\dot\theta^2 + mg\cos\theta - k(r-L)$$ and, $$\dot p_\theta = m(-2r\dot r\dot\theta - gr\sin\theta)$$ Now I have to find the equilibrium solution(s), $x_{eq}$ , of the equations, i.e., those solutions for which $x$ is constant.
How do I find these solutions?