4H2012 lrmte.py

#!/usr/bin/python
from matplotlib import *
from pylab import *

a=int(sys.argv[1])
b=int(sys.argv[2])
n=int(sys.argv[3])
f=sys.argv[4]

def findArea(o,a,b,n,f):
    h=1.0*(b-a)/n
    x=[a+(i+o)*h for i in range(n)]
    y=[f(a+(i+o)*h) for i in range(n)]
    areas=[f(a+(i+o)*h)*h for i in range(n)]
    area=sum(areas)
    return area

print 'n%10s %10s %10s %10s %10s' %('l', 'r', 'm', 't', 'e')
for x in range(5):
    l=findArea(0,a,b,n**x,lambda x:eval(f))
    r=findArea(1,a,b,n**x,lambda x:eval(f))
    m=findArea(0.5,a,b,n**x,lambda x:eval(f))
    t=(r+l)/2
    e=(r-l)/2
    print 'n**%2d %10.4f %10.4f %10.4f %10.4f %10.4f' % (x, l, r, m, t, e)