780 template<class Complex, class F>
781 Complex complex_newton(F g, Complex guess, int max_iterations = std::numeric_limits<typename Comple…  in complex_newton()
783    typedef typename Complex::value_type Real;  in complex_newton()
788    Complex z0 = guess + Complex(1, 0);  in complex_newton()
789    Complex z1 = guess + Complex(0, 1);  in complex_newton()
790    Complex z2 = guess;  in complex_newton()
797          Complex q = (z2 - z1) / (z1 - z0);  in complex_newton()
800          Complex qp1 = static_cast<Complex>(1) + q;  in complex_newton()
801          Complex A = q * (pair.first - qp1 * P1.first + q * P0.first);  in complex_newton()
803 …Complex B = (static_cast<Complex>(2) * q + static_cast<Complex>(1)) * pair.first - qp1 * qp1 * P1.…  in complex_newton()
804          Complex C = qp1 * pair.first;  in complex_newton()
805          Complex rad = sqrt(B * B - static_cast<Complex>(4) * A * C);  in complex_newton()
806          Complex denom1 = B + rad;  in complex_newton()
807          Complex denom2 = B - rad;  in complex_newton()
808          Complex correction = (z1 - z2) * static_cast<Complex>(2) * C;  in complex_newton()