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