/** * copilot suggests: */ #include #include double complex csin(double complex z) { double complex iz = I * z; return (cexp(iz) - cexp(-iz)) / (2 * I); }