Trapezoidal Rule for definite integrals

definttrap.py

a = 0
b = 1
N = 500
area = 0

delta_x = (b - a)*1./N

def f(t):
	return t**2


for i in range(N):

	x = i*delta_x + a
	area=area + 0.5*(f(x)+f(x+delta_x))*delta_x
	
	
print area