Clear Filters
Clear Filters

How to save ode45 outputs in a loop?

3 views (last 30 days)
EB
EB on 9 Oct 2019
Edited: EB on 9 Oct 2019
I am trying to save my T and Y values in a loop, but it gives me the error "Subscripted assignment dimension mismatch."
I originally had [T,Y] = ode45(dydt, t, [y0, dy0], options); but it would only save when j = 7.
Help please?
for j = 1:7
[T(j),Y(j)] = ode45(dydt, t, [y0, dy0]);
end

Accepted Answer

Star Strider
Star Strider on 9 Oct 2019
The easiest way would be to save them as cella rrays, then sort those out later:
for j = 1:7
[T{j},Y{j}] = ode45(dydt, t, [y0, dy0]);
end
Note the curly brackets {} indicating cell array indexing.
It will be easier to work with them if you define ‘tspan’ as a vector with more than two elements:
t = linspace(0, 5, 50); % time span
This will result in all the arrays having equal numbers of rows.

More Answers (1)

James Tursa
James Tursa on 9 Oct 2019
Edited: James Tursa on 9 Oct 2019
Don't use a loop. Just call ode45 once and it will give you the entire results in T and Y.
[T,Y] = ode45(dydt, t, [y0, dy0]);
If you want outputs at specific times, set t accordingly. E.g.,
t = linspace(0,5,100);

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Tags

Products


Release

R2016a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!