How to integrate an array?

[Tom Kaz]

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.
-- 
http://mail.python.org/mailman/listinfo/python-list

Re: How to integrate an array?

[Tom Kaz]

> Do each function separately. The routine that scipy.integrate.quad uses
> adapts to the local conditions of the function (when the function is
> flat, it uses fewer samples; when steep, more).

It's not so easy  to do it separately. I want to integrate function that
includes  linalg.expm - as you see calculation the function is quite
expensive. 

T.Kaz.
-- 
http://mail.python.org/mailman/listinfo/python-list