dsolve(diff(y(x),x) = 1/sinh(x)*cosh(x)+_F1(y(x)-ln(sinh(x))),y(x))