#if REAL > REAL_LONG_DOUBLE
Real error_real = Integrate<Real>(f<Real>,Real(0.0L), Real(1.0L), intervals) - Real(PI);
clock_t clock_real = clock() - start;
#if REAL > REAL_LONG_DOUBLE
Real error_real = Integrate<Real>(f<Real>,Real(0.0L), Real(1.0L), intervals) - Real(PI);
clock_t clock_real = clock() - start;