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