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