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