I need to compute integral of some array function, something like: from scipy import * f = lambda x: array( [sin(x),cos(x)] ) integrate.quad(f, 0, 1) Unfortunately integrate.quad requires Float type. Any ideas how to perform this? Thanks, T.Kaz.