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