extern double sqrt(), f__cabs();
VOID
c_sqrt(r
, z
) complex *r
, *z
;
extern double f__cabs(double, double);
void c_sqrt(complex *r
, complex *z
)
if( (mag
= f__cabs(z
->r
, z
->i
)) == 0.)
r
->r
= t
= sqrt(0.5 * (mag
+ z
->r
) );
t
= sqrt(0.5 * (mag
- z
->r
) );