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[] );
void r8vec_copy ( int n, double a1[], double a2[] );
int rule_order ( int p );
void rule00 ( int n, double a[], double b[], double c[], double w[] );
void rule01 ( int n, double a[], double b[], double c[], double w[] );
void rule02 ( int n, double a[], double b[], double c[], double w[] );
void rule03 ( int n, double a[], double b[], double c[], double w[] );
void rule04 ( int n, double a[], double b[], double c[], double w[] );
void rule05 ( int n, double a[], double b[], double c[], double w[] );
void rule06 ( int n, double a[], double b[], double c[], double w[] );
void rule07 ( int n, double a[], double b[], double c[], double w[] );
void rule08 ( int n, double a[], double b[], double c[], double w[] );
void rule09 ( int n, double a[], double b[], double c[], double w[] );
void rule10 ( int n, double a[], double b[], double c[], double w[] );
void rule11 ( int n, double a[], double b[], double c[], double w[] );
void rule12 ( int n, double a[], double b[], double c[], double w[] );
void rule13 ( int n, double a[], double b[], double c[], double w[] );
void rule14 ( int n, double a[], double b[], double c[], double w[] );
void rule15 ( int n, double a[], double b[], double c[], double w[] );
void rule16 ( int n, double a[], double b[], double c[], double w[] );
void rule17 ( int n, double a[], double b[], double c[], double w[] );
void rule18 ( int n, double a[], double b[], double c[], double w[] );
void rule19 ( int n, double a[], double b[], double c[], double w[] );
void rule20 ( int n, double a[], double b[], double c[], double w[] );
void rule21 ( int n, double a[], double b[], double c[], double w[] );
void rule22 ( int n, double a[], double b[], double c[], double w[] );
void rule23 ( int n, double a[], double b[], double c[], double w[] );
void rule24 ( int n, double a[], double b[], double c[], double w[] );
void rule25 ( int n, double a[], double b[], double c[], double w[] );
void rule26 ( int n, double a[], double b[], double c[], double w[] );
void rule27 ( int n, double a[], double b[], double c[], double w[] );
void rule28 ( int n, double a[], double b[], double c[], double w[] );
void rule29 ( int n, double a[], double b[], double c[], double w[] );
void rule30 ( int n, double a[], double b[], double c[], double w[] );
void rule31 ( int n, double a[], double b[], double c[], double w[] );
void rule32 ( int n, double a[], double b[], double c[], double w[] );
void rule33 ( int n, double a[], double b[], double c[], double w[] );
void rule34 ( int n, double a[], double b[], double c[], double w[] );
void rule35 ( int n, double a[], double b[], double c[], double w[] );
void rule36 ( int n, double a[], double b[], double c[], double w[] );
void rule37 ( int n, double a[], double b[], double c[], double w[] );
void rule38 ( int n, double a[], double b[], double c[], double w[] );
void rule39 ( int n, double a[], double b[], double c[], double w[] );
void rule40 ( int n, double a[], double b[], double c[], double w[] );
void rule41 ( int n, double a[], double b[], double c[], double w[] );
void rule42 ( int n, double a[], double b[], double c[], double w[] );
void rule43 ( int n, double a[], double b[], double c[], double w[] );
void rule44 ( int n, double a[], double b[], double c[], double w[] );
void rule45 ( int n, double a[], double b[], double c[], double w[] );
void rule46 ( int n, double a[], double b[], double c[], double w[] );
void rule47 ( int n, double a[], double b[], double c[], double w[] );
void rule48 ( int n, double a[], double b[], double c[], double w[] );
void rule49 ( int n, double a[], double b[], double c[], double w[] );
void rule50 ( int n, double a[], double b[], double c[], double w[] );
double triangle_unit_area ( );
double triangle_unit_monomial_integral ( int expon[] );
void triangle_symq_rule ( int p, int n, double a[], double b[], 
  double c[], double w[] );