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