dsolve(diff(diff(y(x),x),x) = -(-a^2*sinh(x)^2-n*(n-1))/sinh(x)^2*y(x),y(x))