
Divide both sides by
to get


Substitute
, so that
. Then



The remaining ODE is separable. Separating the variables gives

Integrate both sides. On the left, split up the integrand into partial fractions.




Then

On the right, we have

Solving for
explicitly is unlikely to succeed, so we leave the solution in implicit form,

and finally solve in terms of
by replacing
:


