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