#include "libm.h"
+/* sin(z) = -i sinh(i z) */
+
double complex csin(double complex z)
{
- double complex r;
-
- r = csinh(cpack(-cimag(z), creal(z)));
- return cpack(cimag(r), -creal(r));
+ z = csinh(cpack(-cimag(z), creal(z)));
+ return cpack(cimag(z), -creal(z));
}