#include double sinh(double x) { return (exp(x) - exp(-x)) / 2.0; }