/** * copilot suggests: */ #include #include float complex csinf(float complex z) { float complex iz = I * z; return (cexpf(iz) - cexpf(-iz)) / (2 * I); }