1 2 3 4 5 6 7 8 9
#include "libm.h" /* sin(z) = -i sinh(i z) */ double complex csin(double complex z) { z = csinh(CMPLX(-cimag(z), creal(z))); return CMPLX(cimag(z), -creal(z)); }