double polynomial_root_bound ( int n, complex <double> c[] );
void bisection_rc ( double &a, double &b, double &x, double fx, int &job );
double polyval ( int m, double c[], double x );