# Index exceeds the number of array elements (1)

1 view (last 30 days)
Heya :) on 6 Jan 2021
Commented: Heya :) on 6 Jan 2021
I need to plot x1, y1, x2 and y2 for different values of a1 and a2. Please help me to resolve the issue. Thanks in advance.
clear all; close all; clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%+
% Simulation %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
OMEGA1=1:0.5:14.5;
OMEGA2=1:0.5:14.5;
for num1=1:length(OMEGA1)
for num2=1:length(OMEGA2)
%value of constants
omega1=OMEGA1(num1);
omega2=OMEGA2(num2);
a1=0.1;
a2=0.2;
a3=0.3;
a4=0.4;
a5=0.5;
G=1;C12=0.01;C21=0.02;
dt=0.01; %step size
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x11(1)=0.5;
y11(1)=0.5;
x12(1)=0.5;
y12(1)=0.5;
x13(1)=0.5;
y13(1)=0.5;
x14(1)=0.5;
y14(1)=0.5;
x15(1)=0.5;
y15(1)=0.5;
x21(1)=1;
y21(1)=1;
x22(1)=1;
y22(1)=1;
x23(1)=1;
y23(1)=1;
x24(1)=1;
y24(1)=1;
x25(1)=1;
y25(1)=1;
for i=2:10000
x11(i)=x11(i-1)+((a1-x11(i-1)^2-y11(i-1)^2)*x11(i-1)-omega1*y11(i-1)+G*C12*(x21(i-1)-x11(i-1)))*dt;
x12(i)=x12(i-1)+((a2-x12(i-1)^2-y12(i-1)^2)*x12(i-1)-omega1*y12(i-1)+G*C12*(x22(i-1)-x12(i-1)))*dt;
x13(i)=x13(i-1)+((a3-x13(i-1)^2-y13(i-1)^2)*x13(i-1)-omega1*y13(i-1)+G*C12*(x23(i-1)-x13(i-1)))*dt;
x14(i)=x14(i-1)+((a4-x14(i-1)^2-y14(i-1)^2)*x14(i-1)-omega1*y14(i-1)+G*C12*(x24(i-1)-x14(i-1)))*dt;
x15(i)=x15(i-1)+((a5-x15(i-1)^2-y15(i-1)^2)*x15(i-1)-omega1*y15(i-1)+G*C12*(x25(i-1)-x15(i-1)))*dt;
y11(i)=y11(i-1)+((a1-x11(i-1)^2-y11(i-1)^2)*y11(i-1)+omega1*x11(i-1)+G*C12*(y21(i-1)-y11(i-1)))*dt;
y12(i)=y12(i-1)+((a2-x12(i-1)^2-y12(i-1)^2)*y12(i-1)+omega1*x12(i-1)+G*C12*(y22(i-1)-y12(i-1)))*dt;
y13(i)=y13(i-1)+((a3-x13(i-1)^2-y13(i-1)^2)*y13(i-1)+omega1*x13(i-1)+G*C12*(y23(i-1)-y13(i-1)))*dt;
y14(i)=y14(i-1)+((a4-x14(i-1)^2-y14(i-1)^2)*y14(i-1)+omega1*x14(i-1)+G*C12*(y24(i-1)-y14(i-1)))*dt;
y15(i)=y15(i-1)+((a5-x15(i-1)^2-y15(i-1)^2)*y15(i-1)+omega1*x15(i-1)+G*C12*(y25(i-1)-y15(i-1)))*dt;
x21(i)=x21(i-1)+((a1-x21(i-1)^2-y21(i-1)^2)*x21(i-1)-omega2*y21(i-1)+G*C21*(x11(i-1)-x21(i-1)))*dt;
x22(i)=x22(i-1)+((a2-x22(i-1)^2-y22(i-1)^2)*x22(i-1)-omega2*y22(i-1)+G*C21*(x12(i-1)-x22(i-1)))*dt;
x23(i)=x23(i-1)+((a3-x23(i-1)^2-y23(i-1)^2)*x23(i-1)-omega2*y23(i-1)+G*C21*(x13(i-1)-x23(i-1)))*dt;
x24(i)=x24(i-1)+((a4-x24(i-1)^2-y24(i-1)^2)*x24(i-1)-omega2*y24(i-1)+G*C21*(x14(i-1)-x24(i-1)))*dt;
x25(i)=x25(i-1)+((a5-x25(i-1)^2-y25(i-1)^2)*x25(i-1)-omega2*y25(i-1)+G*C21*(x15(i-1)-x25(i-1)))*dt;
y21(i)=y21(i-1)+((a1-x21(i-1)^2-y21(i-1)^2)*y21(i-1)+omega2*x21(i-1)+G*C21*(y11(i-1)-y21(i-1)))*dt;
y22(i)=y22(i-1)+((a2-x22(i-1)^2-y22(i-1)^2)*y22(i-1)+omega2*x22(i-1)+G*C21*(y12(i-1)-y22(i-1)))*dt;
y23(i)=y23(i-1)+((a3-x23(i-1)^2-y23(i-1)^2)*y23(i-1)+omega2*x23(i-1)+G*C21*(y13(i-1)-y23(i-1)))*dt;
y24(i)=y24(i-1)+((a4-x24(i-1)^2-y24(i-1)^2)*y24(i-1)+omega2*x24(i-1)+G*C21*(y14(i-1)-y24(i-1)))*dt;
y25(i)=y25(i-1)+((a5-x25(i-1)^2-y25(i-1)^2)*y25(i-1)+omega2*x25(i-1)+G*C21*(y15(i-1)-y25(i-1)))*dt;
end
end
end
figure
hold on
plot(x11,'r')
plot(x12,'g')
plot(x13,'r')
plot(x14,'g')
plot(x15,'r')
Ive J on 6 Jan 2021
Except for your misleading script indentation, your script has no errors. You need to exactly clarify what you mean by issue.
Heya :) on 6 Jan 2021
The problem is resolved. Thanks