1 #include "map.h" 2 3 void 4 ccubrt(double zr, double zi, double *wr, double *wi) 5 { 6 double r, theta; 7 theta = atan2(zi,zr); 8 r = cubrt(hypot(zr,zi)); 9 *wr = r*cos(theta/3); 10 *wi = r*sin(theta/3); 11 } 12