Hyperbolic Area Cotangent.
double acoth(double value) { return 0.5 * stdmath.log((value + 1) / (value - 1)); }