3 extern void z_div (doublecomplex *, doublecomplex *, doublecomplex *);
5 pow_zi (doublecomplex * p, doublecomplex * a, integer * b) /* p = a**b */
11 static doublecomplex one = { 1.0, 0.0 };
34 t = q.r * x.r - q.i * x.i;
35 q.i = q.r * x.i + q.i * x.r;
40 t = x.r * x.r - x.i * x.i;