#include double complex catanh(double complex z) { return -1.0 * I * catan(z * I); }