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