How to integrate given the two data columns

Hello. I have used ode45 to figure out current wrt time.

I have my matrix tz and my matrix CURRENT for values of current with time from tz(1)=0 to tz(end). How can I intagrate it? Any clever way if I do not know the time step in tz?

Answer 1

[Riccardo Scorretti]

Most probably like this:

CHARGE = sum(CURRENT(1:end-1) .* diff(tz));

[Riccardo Scorretti]

Something better can be done, but if you have a large number of points, it will give more or less the same result.

Answer 2

[Torsten]

Add the differential equation

dz/dt = CURRENT, z(t0) = 0

to your system of ODEs.

z will automatically be the integrated CURRENT over time.