dsolve(diff(y(x),x) = (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+x*y(x))/x^2,y(x))