extern void sig_die(char*,int);
void c_div(complex *c
, complex *a
, complex *b
)
sig_die("complex division by zero", 1);
ratio
= (double)b
->r
/ b
->i
;
den
= b
->i
* (1 + ratio
*ratio
);
c
->r
= (a
->r
*ratio
+ a
->i
) / den
;
c
->i
= (a
->i
*ratio
- a
->r
) / den
;
ratio
= (double)b
->i
/ b
->r
;
den
= b
->r
* (1 + ratio
*ratio
);
c
->r
= (a
->r
+ a
->i
*ratio
) / den
;
c
->i
= (a
->i
- a
->r
*ratio
) / den
;