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