#include "vs.h"
int main() {
	double x = 0.0,
          w[11] = {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};
   for(int i = 0; i < 11; i++) {
   	Quadrature q1(w, 0.0, x, 11), q2(w, x, 1.0, 11);
      H0 z1(q1), z2(q2); C0 f1, f2;
      if(i != 0) f1 &= ( (z1 * log(1+z1)) | J(x / 10.0)); else f1 &= C0(0.0);
      if(i != 10) f2 &= log(1+x) * ( z2 | J((1.0-x) / 10.0)); else f2 &= C0(0.0);
      double u = (double) (f1 + f2);
   	cout << "u(" << x << "): " << u << "analytical: " << (-(x*x)/4+x/2) << endl;
   	if(i != 10) x += 0.1;
   }
	return 0;
}