solve dx/dt = -(a + bx^2)