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