#include "vs.h"
int main() {
	const double PI = 3.141592654;
   const double alpha = 0.0;
   const double beta = 0.0;
	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),
#if defined(__COS)
			f1 = cos(PI*z1), f2 = cos(PI*z2);
#else
         f1 = sin(PI*z1), f2 = sin(PI*z2);
#endif
      C0 integ_1, integ_2;
      if(i != 0) integ_1 &= (1-x) * ( (z1 * f1) | J(x / 10.0)); else integ_1 &= C0(0.0);
      if(i != 10) integ_2 &= x * ( (1-z2) * f2) | J((1.0-x) / 10.0); else integ_2 &= C0(0.0);
      double u = (double) (integ_1 + integ_2) + (1-x) * alpha + x * beta;
   	cout << "u(" << x << "): " << u << endl;
   	if(i != 10) x += 0.1;
   }
	return 0;
}