Bifurcation diagram for a non autonomous system.

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);

Answer 1

[Jyotish Robin]

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.

Read more about this at: https://www.mathworks.com/help/matlab/matlab_prog/debugging-process-and-features.html

Hope this helps.

Regards,

Jyotish

[LALE ASIK]

Thank you so much.