void zero_muller ( double complex func ( double complex x ), double fatol,
  int itmax, double complex x1, double complex x2, double complex x3, 
  double xatol, double xrtol, double complex *xnew, double complex *fxnew );