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