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