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