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