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