/*
  Numerical integration of f(x) on [a,b] using quanc8.cpp
  AG: February 2007
*/
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
#include "quanc8.cpp"

using namespace std;

double f(double);

int main ()
{
    double a, b;
    double result, errest, flag;
    int nofun;
    const double pi = 3.1415926;
    const double abserr=0.0, relerr=1.0e-10;  //see quanc8.cpp

    cout.precision(6);
    cout.setf(ios::fixed | ios::showpoint);
    
    a = 0.0;             // left  endpoint
    b = 1.0*pi;          // right endpoint

    quanc8(f, a, b, abserr, relerr, result, errest, nofun, flag);

    cout << setw(12) << "nofun = "    << setw(10) << nofun  << endl;
    cout << setw(12) << "integral = " << setw(10) << result << endl;
    cout << setw(12) << "est.error = "<< setw(10) << errest << endl;
    cout << setw(12) << "flag = "     << setw(10) << flag   << endl;

    system ("pause");
    return 0;
}

/*
 *  Function f(x)
*/
    double f(double x)
{
    double y;
    y = sin(x);
    return y;
}

/* Test output:
         
    nofun =         33
 integral =   2.000000
est.error =   0.000000
     flag =   0.000000

*/