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