void roots_rc ( int n, double x[], double fx[], double &ferr, double xnew[], 
  double q[] );
void r8mat_fs ( int lda, int n, double a[], double x[] );
void timestamp ( );