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