'''
Created on 2012-10-18
@author: User
'''
from scipy.integrate import quad
from pylab import *
# define a integrand
def intgrnd1(x):
'''
return an integrands
'''
return cos(x)**2
#//def
# define a integral as a function of its upper limit and vectorizes it
def intrgl(f_up_limit,f_low_limit=0):
'''
return an integral as a functjion of its upper limit
the lower limit is 0
'''
return quad(intgrnd1, f_low_limit, f_up_limit)
# vectorize function intrgl
v_intrgl=vectorize(intrgl)
#================== output ======================
# get results
lst_up_limits=np.arange(0,10,1)
results, errs = v_intrgl(lst_up_limits)
#plot results
figure()
plot(lst_up_limits, results)
show()
print results
Reference: