Bifurcation diagram for a non autonomous system.

6 views (last 30 days)
I am trying to draw a bifurcation diagram for a non-autonomous system. The code:
limit=[1400:1:1500];
ymin=[];
ymax=[];
for Ke=0:0.01:3;
[~,U]=ode45(@Onetwo,[0:1:1500],[0.5;0.25;1.5]);
u2=U(limit,2);
The system gives me the following error:
Index exceeds matrix dimensions.
Error in LKEbif (line 24)
u2=U(limit,2);

Accepted Answer

Jyotish Robin
Jyotish Robin on 10 Oct 2017
Hi!
The reason for the error "Index exceeds matrix dimensions." in your piece of code is because U is an array of size 51 x 3 but you are trying to use 'limit' which range from 1400 to 1500 to index into U.
U(a,b) % will be valid only if a<52 and b<4 where a and b are positive integers OR Logicals.
In cases of trivial issues like this, I would suggest you to try debugging using breakpoints in MATLAB.
Hope this helps.
Regards,
Jyotish

More Answers (0)

Community Treasure Hunt

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

Start Hunting!