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