#include double atanh(double x) { return log(1 + x) / 2.0; }