solve ODE (y/y')^2+y^2=b^2(x-(y/y'))^2