Calculate the trapezoidal sum of a function over an interval. An external y function needs to be created. Uses numpy

trapezoidalSum.py

def trapezoidalSum(xRange,n):
	tot = 0
	width = 0
	width = (xRange[1]-xRange[0])/n
	x = [ interval for interval in arange(xRange[0],xRange[1]+width, width) ]
	yRange = [ y(n) for n in x ]
	tot += yRange[0]+yRange[n]
	tot += sum([ 2*yRange[i] for i in range(1,n) ])
	tot = tot*width/2
	return tot