/* triquad_int.h  1 <= n <=6 at this time                      */
/*                nn returned is n*n                           */
/*                integral = sum i=0 to nn-1 W[i]*f(X[i],Y[i]) */

int triquad_int(int n, double x1, double y1,
                       double x2, double y2,
                       double x3, double y3,
                       double X[], double y[], double W[]);