#include "vs.h"
int main() {
	double x[3] = {1.0, 1.5, 3.0};
	Quadrature qp(1, 3);
   H1 zai(qp),
      X = -(1-zai)*zai/2.0 * x[0]
          +(1-zai)*(1+zai) * x[1]
          +zai*(1+zai)/2.0 * x[2];
   J d_l(d(X));
   C0 v = (((H0)X)-1).pow(2) | d_l;
   cout << v << endl;
   return 0;
}