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 );