Flute acoustics: measurement, modelling and design - School of ...

APPENDIX B. PROGRAM LISTINGS 235
}
return divz(sinz(z1), cosz(z1));
complex cotz(complex z1) {
return divz(cosz(z1), sinz(z1));
}
complex logz(complex z1) {
complex result;
result.Re = log(modz(z1));
result.Im = argz(z1);
return result;
}
complex arctanz(complex z1) {
complex two = {2.0, 0.0};
return multz(divz(j, two), logz(divz(addz(j, z1), subz(j, z1))));
}
complex parallel(complex z1, complex z2) {
return divz(multz(z1, z2), addz(z1, z2));
}
complex sqrtz(complex z1) {
double A = modz(z1);
double theta = argz(z1);
A = sqrt(A);
theta = theta / 2;
return multz(real(A), expjz(real(theta)));
}
void printComplex(complex z) {
printf("%e%+ei", z.Re, z.Im);
}