1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
#include "f2c.h" #ifdef KR_headers double log(), f__cabs(), atan2(); VOID z_log(r, z) doublecomplex *r, *z; #else #undef abs #include "math.h" extern double f__cabs(double, double); void z_log(doublecomplex *r, doublecomplex *z) #endif { double zi = z->i, zr = z->r; r->i = atan2(zi, zr); r->r = log( f__cabs( zr, zi ) ); }