#include "vs.h"
int main() {
   double f[4] = {1.0, 2.0, 2.0, 3.0};
   Quadrature qp(2, 4);
   H0 z(2, (double*)0, qp), zai, eta;
   zai &= z[0], eta &= z[1];
   C0 volume = ((1-zai)*(1-eta)/4*f[0]
               +(1+zai)*(1-eta)/4*f[1]
               +(1+zai)*(1+eta)/4*f[2]
               +(1-zai)*(1+eta)/4*f[3]) | J(0.25);
   cout << volume << endl;
	return 0;
}