#include // for csqrt, cacosh, clog, complex double complex cacosh(double complex z) { return clog(z + csqrt(z + 1) * csqrt(z - 1)); }