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