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