4H2012 lsumtable.py

#!/usr/bin/python
a=2
b=9
def f(x):
    return x**3

#print "h = ", h
#print range(n)
#print [i for i in range(n)]
#print [i*h for i in range(n)]
#print [a+i*h for i in range(n)]
#print [f(a+i*h) for i in range(n)]
#print [f(a+i*h)*h for i in range(n)]
#print sum([f(a+i*h)*h for i in range(n)])

print "n\tlsum(n)"
for x in range(6):
    h=1.0*(b-a)/10**x
    print "10**",x,"\t",sum([f(a+i*h)*h for i in range(10**x)])