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