# test_gauleg.py  test  gauleg.py  function gaulegf
from gauleg import gaulegf

def f(p):
  return p*p

print "test_gauleg.py integrate f(p)=p^2 from 1 to 2"
print "calling x,w = gaulegf(1.0, 2.0, 4)"
n = 4
x,w = gaulegf(1.0, 2.0, n)
print "x=",
print x
print "w=",
print w

area = 0.0
for i in range(1, n+1):
  area += w[i]*f(x[i])

print "exact area=7/3"
print "computed area= ",
print area
print "error= ",
print area-7.0/3.0