extended_simpson.cpp

算法的一些集合

#include "vs.h"
int main() {
   const int N = 11;
   double c[N] = {1.0/3.0, 4.0/3.0, 2.0/3.0, 4.0/3.0, 2.0/3.0, 4.0/3.0,
   					2.0/3.0, 4.0/3.0, 2.0/3.0, 4.0/3.0, 1.0/3.0},
                X[2] = {1.0, 2.0},
                h = (X[1]-X[0])/N;
   Quadrature qp(c, X[0], X[1], N);
   H0 x(qp);
   C0 area = (x.pow(2)-2*x+1) | J(h);
   cout << area << endl;
	return 0;
}