int cg_rc ( int n, double b[], double x[], double r[], double z[], 
  double p[], double q[], int job );
void r8mat_mv ( int m, int n, double a[], double x[], double ax[] );
double *r8vec_uniform_01_new ( int n, int &seed );
void timestamp ( );
double *wathen ( int nx, int ny, int n );
int wathen_order ( int nx, int ny );