How to save ode45 outputs in a loop?
3 views (last 30 days)
Show older comments
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
0 Comments
Accepted Answer
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.
0 Comments
More Answers (1)
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);
0 Comments
See Also
Categories
Find more on Loops and Conditional Statements in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!