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