integral.py

from random import uniform
#calculate intergral
def integral(fun,a,b,N=10000):
	sum=0
	for i in range(N):
		sum=sum+fun(uniform(a,b))
	return sum/N*(b-a)

# test
print integral(lambda x : x + 1,0,1)