dsolve(diff(y(x),x) = F(-(-y(x)^2+b)/x^2)*x/y(x),y(x))