#include #include double cabs(double complex z) { return hypot(creal(z), cimag(z)); }