#include <iostream>

#include <Function.H>
#include <GaussKronod.H>

struct Function {

    typedef double Value;

    Function(const double p=0.0): a(p) { }

    double operator()(const double x) const {
        return x*x-a;
    }

    double a;
};

int
main() {

    typedef Integration::GaussKronod<21>   Integrator;
    typedef Integrator::Integral<Function> Integral;

    Function f(0.5);
    Integrator integrate;
    Integral integral = integrate(f,-1.0,1.0);

    std::cout << "Integral     = " << integral.integral() << std::endl
              << "Abs integral = " << integral.abs()      << std::endl
              << "Asc integral = " << integral.asc()      << std::endl
              << "Err integral = " << integral.error()    << std::endl;
}