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