#include double acosh(double x) { return log(x + sqrt(x*x - 1)); }