// test_gaulegf.cpp  test  gaulegf.c 
#include <iostream>
using namespace std;
#include "gaulegf.h"

static double f(double p)
{
  return p*p;
}

int main(int argc, char * argv[])
{
  int n = 4;
  double x[5];
  double w[5];
  double area;

  cout << "test_gaulegf.cpp running" << endl;
  cout << "calling gaulegf(1.0, 2.0, x, w, n)" << endl;
  gaulegf(1.0, 2.0, x, w, n);
  for(int i=1; i<=n; i++)
    cout << "x[" << i << "]=" << x[i] << ", w[" << i << "]=" << w[i] << endl;
	
  area = 0.0;
  for(int i=1; i<=n; i++)
  {
    area += w[i]*f(x[i]);
  }
  cout << "area= " << area << endl;
  cout << "error=" << area-(7.0/3.0) << endl;
  return 0;
} // end test_gaulegf.cpp