void legendre_3d_exactness ( double a[], double b[], int n, double x[], 
  double y[], double z[], double w[], int t );
double legendre_3d_monomial_integral ( double a[], double b[], int p[] );
double r8vec_dot_product ( int n, double a1[], double a2[] );
void timestamp ( );