/* test_gaulegf.c  test  gaulegf.c */
#include "gaulegf.h"
#include <stdio.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;
  int i;

  printf("calling gaulegf(1.0, 2.0, x, w, n) \n");
  gaulegf(1.0, 2.0, x, w, n);
  for(i=1; i<=n; i++)
    printf("x[%d]=%f, w[%d]=%f \n", i, x[i], i, w[i]);
	
  area = 0.0;
  for(i=1; i<=n; i++)
  {
    area += w[i]*f(x[i]);
  }
  printf("area=%f \n", area);
  printf("error=%e \n", area-(7.0/3.0));
  return 0;
} /* end test_gaulegf.c */