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