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