void power_method ( int n, double a[], double y[], int it_max, double tol,
  double *lambda, int *it_num );
void r8mat_mv ( int m, int n, double a[], double x[], double ax[] );
double r8vec_dot ( int n, double a1[], double a2[] );
double r8vec_norm_l2 ( int n, double a[] );