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