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