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

def f(p):
  return p*p

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 "area=",
print  area