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