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