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