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