solve dz(x) / dx-(z(x))^2 = -a^2 + b/x