[SciPy-User] calculate definite integral of sampled data
Charles R Harris
Wed Aug 10 20:26:32 CDT 2011
On Tue, Aug 9, 2011 at 11:59 PM, Manuel Graune <email@example.com> wrote:
> Hi everyone,
> to calculate the definite integral of a function or an array of sampled
> data scipy provides (among others) the quad and trapz functions.
> So it is possible to compute e. g. the definite integral of cos(t) over
> some interval by doing
> definite_integral= scipy.integrate.quad(cos,lower_limit,upper_limit)
> definite_integral= scipy.integrate.trapz(some_array).
> Now, if I want to plot cos(t) and the integral of cos(t) from 0 to t in
> a graph, the necessary array can be calculated by:
> def intfunc(fnc,upper_limit):
> return scipy.integrate.quad(fnc,0.0,upper_limit)
> definite_inegral= intfunc(cos,t)
> which seems (whithout knowing the actual code) a bit wasteful and slow
> but is relatively concise.
> Now for my question: scipy provides e. g. the trapz-function to
> calculate definite integral of a complete array of sampled data.
> However, I have no idea how to get achieve the same as above for
> sampled data (apart from manually iterating in a for-loop). Is there
> a function somewhere which delivers an array of the definite integrals
> for each of the data-points in an array?
How many data points do you have?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User