#include double cosh (double x) { return (exp(x) + exp(-x)) / 2; }