The hyperbolic sine function using angle
angle
BigDouble sinh(BigDouble angle) { return (exp(angle) - exp(-angle)) / 2.big; }