#include "f2c.h" #ifdef KR_headers double log(), f__cabs(), atan2(); VOID z_log(resx, z) doublecomplex *resx, *z; #else #undef abs #include extern double f__cabs(double, double); void z_log(doublecomplex *resx, doublecomplex *z) #endif { doublecomplex res; res.i = atan2(z->i, z->r); res.r = log( f__cabs( z->r, z->i ) ); resx->r = res.r; resx->i = res.i; }