I am given the following ode:
$$x''=mg\sin(\theta)+k(L-\sqrt{x^2+h^2})\left(\frac{x}{\sqrt{x^2+h^2}}\right)-bx'$$ where $x$ is a function of $t.$ How do I convert this equation to a system of first order odes? I tried the usual technique of substituting $x_1=x$ and $x_2=x'$ but it does not work since the ode is non-linear. Perhaps there is another way?