void comp_next ( int n, int k, int a[], int *more, int *h, int *t );
double *monomial_value ( int m, int n, int e[], double x[] );
double quadrilateral_unit_area ( );
double quadrilateral_unit_monomial_integral ( int expon[2] );
void r8vec_copy ( int n, double a1[], double a2[] );
double r8vec_sum ( int n, double a[] );
int rule_order ( int degree );
void rule00 ( int n, double x[], double y[], double w[] );
void rule01 ( int n, double x[], double y[], double w[] );
void rule02 ( int n, double x[], double y[], double w[] );
void rule03 ( int n, double x[], double y[], double w[] );
void rule04 ( int n, double x[], double y[], double w[] );
void rule05 ( int n, double x[], double y[], double w[] );
void rule06 ( int n, double x[], double y[], double w[] );
void rule07 ( int n, double x[], double y[], double w[] );
void rule08 ( int n, double x[], double y[], double w[] );
void rule09 ( int n, double x[], double y[], double w[] );
void rule10 ( int n, double x[], double y[], double w[] );
void rule11 ( int n, double x[], double y[], double w[] );
void rule12 ( int n, double x[], double y[], double w[] );
void rule13 ( int n, double x[], double y[], double w[] );
void rule14 ( int n, double x[], double y[], double w[] );
void rule15 ( int n, double x[], double y[], double w[] );
void rule16 ( int n, double x[], double y[], double w[] );
void rule17 ( int n, double x[], double y[], double w[] );
void rule18 ( int n, double x[], double y[], double w[] );
void rule19 ( int n, double x[], double y[], double w[] );
void rule20 ( int n, double x[], double y[], double w[] );
void square_symq ( int degree, int n, double x[], double y[], double w[] );
void square_symq_gnuplot ( int n, double x[], double y[], char *header );