Please help me to spot the error.

Question

Haya Ali 2023-1-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1899830-please-help-me-to-spot-the-error

回答： Voss 2023-1-24

I have made a code for 23 coupled oscillators. Matlab is showing error but its not telling the what is the error. Below is my code.

clear all; close all; clc;
% value of constants%
G      = 1;
dt     = 0.01;   % step size
%Oscillator 1%
a1 = 0.10; omega1 = 20; C1_2 =0; C1_3 =0.010; C1_4 =0.0101; C1_5 =0; C1_6 =0.0102; C1_7 =0.0103; C1_8 =0; C1_9 =0.0104; C1_10 =0.0105; C1_11 =0; C1_12 =0.0106; C1_13 =0.0107; C1_14 =0; C1_15 =0.0108; C1_16 =0.0109; C1_17 =0; C1_18 =0.01010; C1_19 =0.01011; C1_20 =0; C1_21 =0.01012; C1_22 =0.01013; C1_23 =0;
%Oscillator 2%
a2 = 0.3; omega2 = 10; C2_1 =0.03; C2_3 =0; C2_4 =0.031; C2_5 =0.032; C2_6 =0; C2_7 =0.033; C2_8 =0.034; C2_9 =0; C2_10 =0.035; C2_11 =0.036; C2_12 =0; C2_13 =0.037; C2_14 =0.038; C2_15 =0; C2_16 =0.039; C2_17 =0.0310; C2_18 =0; C2_19 =0.0311; C2_20 =0.0312; C2_21 =0; C2_22 =0.0313; C2_23 =0.0314;
%Oscillator 3%
a3 = 0.8; omega3 = 40; C3_1 =0.08; C3_2 =0.081; C3_4 =0; C3_5 =0.082; C3_6 =0.083; C3_7 =0; C3_8 =0.084; C3_9 =0.085; C3_10 =0; C3_11 =0.086; C3_12 =0.087; C3_13 =0; C3_14 =0.088; C3_15 =0.089; C3_16 =0; C3_17 =0.0810; C3_18 =0.0811; C3_19 =0; C3_20 =0.0812; C3_21 =0.0813; C3_22 =0; C3_23 =0.0814;
%Oscillator 4%
a4 = 0.5; omega4 = 30; C4_1 =0; C4_2 =0.05; C4_3 =0.051; C4_5 =0; C4_6 =0.052; C4_7 =0.053; C4_8 =0; C4_9 =0.054; C4_10 =0.055; C4_11 =0; C4_12 =0.056; C4_13 =0.057; C4_14 =0; C4_15 =0.058; C4_16 =0.059; C4_17 =0; C4_18 =0.0510; C4_19 =0.0511; C4_20 =0.0512; C4_21 =0.0513; C4_22 =0.0514; C4_23 =0;
%Oscillator 5%
a5 = 0.6; omega5 = 60; C5_1 =0.06; C5_2 =0; C5_3 =0.061; C5_4 =0.062; C5_6 =0; C5_7 =0.063; C5_8 =0.064; C5_9 =0; C5_10 =0.065; C5_11 =0.066; C5_12 =0.067; C5_13 =0.068; C5_14 =0.069; C5_15 =0; C5_16 =0.0610; C5_17 =0.0611; C5_18 =0; C5_19 =0.0612; C5_20 =0.0613; C5_21 =0; C5_22 =0.0614; C5_23 =0.0615;
%Oscillator 6%
a6 = 0.1; omega6 = 50; C6_1 =0.01; C6_2 =0.011; C6_3 =0; C6_4 =0.012; C6_5 =0.013; C6_7 =0; C6_8 =0.014; C6_9 =0.015; C6_10 =0; C6_11 =0.016; C6_12 =0.017; C6_13 =0; C6_14 =0.018; C6_15 =0.019; C6_16 =0; C6_17 =0.0110; C6_18 =0.0111; C6_19 =0; C6_20 =0.0112; C6_21 =0.0113; C6_22 =0; C6_23 =0.0114;
%Oscillator 7%
a7 = 0.4; omega7 = 80; C7_1 =0; C7_2 =0.04; C7_3 =0.041; C7_4 =0; C7_5 =0.042; C7_6 =0.043; C7_8 =0; C7_9 =0.044; C7_10 =0.045; C7_11 =0; C7_12 =0.046; C7_13 =0.047; C7_14 =0; C7_15 =0.048; C7_16 =0.049; C7_17 =0; C7_18 =0.0410; C7_19 =0.0411; C7_20 =0; C7_21 =0.0412; C7_22 =0.0413; C7_23 =0;
%Oscillator 8%
a8 = 0.2; omega8 = 70; C8_1 =0.02; C8_2 =0; C8_3 =0.021; C8_4 =0.022; C8_5 =0; C8_6 =0.023; C8_7 =0.024; C8_9 =0; C8_10 =0.025; C8_11 =0.026; C8_12 =0; C8_13 =0.027; C8_14 =0.028; C8_15 =0; C8_16 =0.029; C8_17 =0.0210; C8_18 =0; C8_19 =0.0211; C8_20 =0.0212; C8_21 =0; C8_22 =0.0213; C8_23 =0.0214;
%Oscillator 9%
a9 = 0.7; omega9 = 60; C9_1 =0.07; C9_2 =0.071; C9_3 =0; C9_4 =0.072; C9_5 =0.073; C9_6 =0; C9_7 =0.074; C9_8 =0.075; C9_10 =0; C9_11 =0.076; C9_12 =0.077; C9_13 =0; C9_14 =0.078; C9_15 =0.079; C9_16 =0; C9_17 =0.0710; C9_18 =0.0711; C9_19 =0; C9_20 =0.0712; C9_21 =0.0713; C9_22 =0; C9_23 =0.0714;
%Oscillator 10%
a10 = 0.10; omega10 = 90; C10_1 =0; C10_2 =0.09; C10_3 =0.091; C10_4 =0; C10_5 =0.092; C10_6 =0.093; C10_7 =0; C10_8 =0.094; C10_9 =0.095; C10_11 =0; C10_12 =0.096; C10_13 =0.097; C10_14 =0; C10_15 =0.098; C10_16 =0.099; C10_17 =0; C10_18 =0.0910; C10_19 =0.0911; C10_20 =0; C10_21 =0.0912; C10_22 =0.0913; C10_23 =0;
%Oscillator 11%
a11 = 0.20; omega11 = 20; C11_1 =0.020; C11_2 =0; C11_3 =0.0201; C11_4 =0.0202; C11_5 =0; C11_6 =0.0203; C11_7 =0.0204; C11_8 =0; C11_9 =0.0205; C11_10 =0.0206; C11_12 =0; C11_13 =0.0207; C11_14 =0.0208; C11_15 =0; C11_16 =0.0209; C11_17 =0.02010; C11_18 =0; C11_19 =0.02011; C11_20 =0.02012; C11_21 =0; C11_22 =0.02013; C11_23 =0.02014;
%Oscillator 12%
a12 = 0.13; omega12 = 10; C12_1 =0.013; C12_2 =0.0131; C12_3 =0; C12_4 =0.0132; C12_5 =0.0133; C12_6 =0; C12_7 =0.0134; C12_8 =0.0135; C12_9 =0; C12_10 =0.0136; C12_11 =0.0137; C12_13 =0; C12_14 =0.0138; C12_15 =0.0139; C12_16 =0; C12_17 =0.01310; C12_18 =0.01311; C12_19 =0; C12_20 =0.01312; C12_21 =0.01313; C122_2 =0; C12_23 =0.01314;
%Oscillator 13%
a13 = 0.18; omega13 = 40; C13_1 =0.018; C13_2 =0.0181; C13_3 =0.0182; C13_4 =0; C13_5 =0.0183; C13_6 =0.0184; C13_7 =0; C13_8 =0.0185; C13_9 =0.0186; C13_10 =0; C13_11 =0.0187; C13_12 =0.0188; C13_14 =0; C13_15 =0.0189; C13_16 =0.01810; C13_17 =0; C13_18 =0.01811; C13_19 =0.01812; C13_20 =0; C13_21 =0.01813; C13_22 =0.01814; C13_23 =0;
%Oscillator 14%
a14 = 0.23; omega14 = 30; C14_1 =0.023; C14_2 =0.0231; C14_3 =0; C14_4 =0.0232; C14_5 =0.0233; C14_6 =0; C14_7 =0.0234; C14_8 =0.0235; C14_9 =0; C14_10 =0.0236; C14_11 =0; C14_12 =0.0237; C14_13 =0.0238; C14_15 =0; C14_16 =0.0239; C14_17 =0.02310; C14_18 =0; C14_19 =0.02311; C14_20 =0.02312; C14_21 =0; C14_22 =0.02313; C14_23 =0.02314;
%Oscillator 15%
a15 = 0.22; omega15 = 60; C15_1 =0.022; C15_2 =0; C15_3 =0.0221; C15_4 =0; C15_5 =0.0222; C15_6 =0.0223; C15_7 =0; C15_8 =0.0224; C15_9 =0.0225; C15_10 =0; C15_11 =0.0226; C15_12 =0; C15_13 =0.0227; C15_14 =0.0228; C15_16 =0; C15_17 =0.0229; C15_18 =0.02210; C15_19 =0; C15_20 =0.02211; C15_21 =0.02212; C15_22 =0; C15_23 =0.02213;
%Oscillator 16%
a16 = 0.11; omega16 = 50; C16_1 =0; C16_2 =0.011; C16_3 =0.0111; C16_4 =0; C16_5 =0.0112; C16_6 =0.0113; C16_7 =0; C16_8 =0.0114; C16_9 =0.0115; C16_10 =0; C16_11 =0.0116; C16_12 =0.0117; C16_13 =0; C16_14 =0.0118; C16_15 =0.0119; C16_17 =0; C16_18 =0.01110; C16_19 =0.01111; C16_20 =0; C16_21 =0.01112; C16_22 =0.01113; C16_23 =0;
%Oscillator 17%
a17 = 0.14; omega17 = 80; C17_1 =0.014; C17_2 =0; C17_3 =0.0141; C17_4 =0.0142; C17_5 =0; C17_6 =0.0143; C17_7 =0.0144; C17_8 =0; C17_9=0.0145; C17_10 =0.0146; C17_11 =0; C17_12 =0.0147; C17_13 =0.0148; C17_14 =0; C17_15 =0.0149; C17_16 =0.01410; C17_18 =0; C17_19 =0.01411; C17_20 =0.01412; C17_21 =0; C17_22 =0.01413; C17_23 =0.01414;
%Oscillator 18%
a18 = 0.21; omega18 = 70; C18_1 =0.021; C18_2 =0.0211; C18_3 =0; C18_4 =0.0212; C18_5 =0.0212; C18_6 =0; C18_7 =0.0215; C18_8 =0.0216; C18_9 =0; C18_10 =0.0217; C18_11 =0.0218; C18_12 =0; C18_13 =0.0219; C18_14 =0.02110; C18_15 =0; C18_16 =0.02111; C18_17 =0.02112; C18_19 =0; C18_20 =0.02113; C18_21 =0.02114; C18_22 =0; C18_23 =0.02115;
%Oscillator 19%
a19 = 0.17; omega19 = 10; C19_1 =0; C19_2 =0.017; C19_3 =0.0171; C19_4 =0; C19_5 =0.0172; C19_6 =0.0173; C19_7 =0; C19_8 =0.0174; C19_9 =0.0175; C19_10 =0; C19_11 =0.0176; C19_12 =0.0177; C19_13 =0; C19_14 =0.0178; C19_15 =0.0179; C19_16 =0; C19_17 =0.01710; C19_18 =0.01711; C19_20 =0; C19_21 =0.01712; C19_22 =0.01713; C19_23 =0;
%Oscillator 20%
a20 = 0.19; omega20 = 90; C20_1 =0.019; C20_2 =0; C20_3 =0.0191; C20_4 =0.0192; C20_5 =0; C20_6 =0.0193; C20_7 =0.0194; C20_8 =0; C20_9 =0.0195; C20_10 =0.0196; C20_11 =0; C20_12 =0.0197; C20_13 =0.0198; C20_14 =0; C20_15 =0.0199; C20_16 =0.01910; C20_17 =0; C20_18 =0.01911; C20_19 =0.01912; C20_21 =0; C20_22 =0.01913; C20_23 =0.01914;
%Oscillator 21%
a21 = 0.15; omega21 = 20; C21_1 =0.015; C21_2 =0.0151; C21_3 =0; C21_4 =0.0152; C21_5 =0.0153; C21_6 =0; C21_7 =0.0154; C21_8 =0.0155; C21_9 =0; C21_10 =0.0156; C21_11 =0.0157; C21_12 =0; C21_13 =0.0158; C21_14 =0.0159; C21_15 =0; C21_16 =0.01510; C21_17 =0.01511; C21_18 =0; C21_19 =0.01512; C21_20 =0.01513; C21_22 =0; C21_23 =0.01514;
%Oscillator 22%
a22 = 0.12; omega22 = 10; C22_1 =0; C22_2 =0.012; C22_3 =0.0121; C22_4 =0; C22_5 =0.0122; C22_6 =0.0123; C22_7 =0; C22_8 =0.0124; C22_9 =0.0125; C22_10 =0; C22_11 =0.0126; C22_12 =0.0127; C22_13 =0; C22_14 =0.0128; C22_15 =0.0129; C22_16 =0; C22_17 =0.01210; C22_18 =0.01211; C22_19 =0; C22_20 =0.01212; C22_21 =0.01213; C22_23 =0;
%Oscillator 23%
a23 = 0.16; omega23 = 40; C23_1 =0.016; C23_2 =0; C23_3 =0.0161; C23_4 =0.0162; C23_5 =0; C23_6 =0.0163; C23_7 =0.0164; C23_8 =0; C23_9 =0.0165; C23_10 =0.0166; C23_11 =0; C23_12 =0.0167; C23_13 =0.0168; C23_14 =0; C23_15 =0.0169; C23_16 =0.01610; C23_17 =0; C23_18 =0.01611; C23_19 =0.01612; C23_20 =0; C23_21 =0.01613; C23_22 =0.01614;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x1(1) = 0.5; y1(1) = 0.5; x2(1) = 0.5; y2(1) = 0.5; x3(1) = 0.5; y3(1) = 0.5; x4(1) = 0.5; y4(1) = 0.5; x5(1) = 0.5; y5(1) = 0.5; x6(1) = 0.5; y6(1) = 0.5; x7(1) = 0.5; y7(1) = 0.5; x8(1) = 0.5; y8(1) = 0.5; x9(1) = 0.5; y9(1) = 0.5; x10(1) = 0.5; y10(1) = 0.5; 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; x16(1) = 0.5; y16(1) = 0.5; x17(1) = 0.5; y17(1) = 0.5; x18(1) = 0.5; y18(1) = 0.5; x19(1) = 0.5; y19(1) = 0.5; x20(1) = 0.5; y20(1) = 0.5; x21(1) = 0.5; y21(1) = 0.5; x22(1) = 0.5; y22(1) = 0.5; x23(1) = 0.5; y23(1) = 0.5;
for i = 2:1000
    %Oscillator 1
    x1(i) = x1(i-1) + ( ( a1 - x1(i-1)^2 - y1(i-1)^2 )*x1(i-1) - omega1*y1(i-1) + G*C1_2*( x2(i-1) - x1(i-1) ) + G*C1_3*( x3(i-1) - x1(i-1) ) + G*C1_4*( x4(i-1) ) - x1(i-1) + G*C1_5*( x5(i-1) - x1(i-1) ) + G*C1_6*( x6(i-1) - x1(i-1) ) + G*C1_7*( x7(i-1) - x1(i-1) ) + G*C1_8*( x8(i-1) - x1(i-1) ) + G*C1_9*( x9(i-1) - x1(i-1) ) + G*C1_10*( x10(i-1) - x1(i-1) ) + G*C1_11*( x11(i-1) - x1(i-1) ) + G*C1_12*( x12(i-1) - x1(i-1) ) + G*C1_13*( x13(i-1) - x1(i-1) ) + G*C1_14*( x14(i-1) - x1(i-1) ) + G*C1_15*( x15(i-1) - x1(i-1) ) + G*C1_16*( x16(i-1) - x1(i-1) ) + G*C1_17*( x17(i-1) - x1(i-1) ) + G*C1_18*( x18(i-1) - x1(i-1) ) + G*C1_19*( x19(i-1) - x1(i-1) ) + G*C1_20*( x20(i-1) - x1(i-1) ) + G*C1_21*( x21(i-1) - x1(i-1) ) + G*C1_22*( x22(i-1) - x1(i-1) ) + G*C1_23*( x23(i-1) - x1(i-1) ) )*dt;
    y1(i) = y1(i-1) + ( ( a1 - x1(i-1)^2 - y1(i-1)^2 )*y1(i-1) + omega1*x1(i-1) + G*C1_2*( y2(i-1) - y1(i-1) ) + G*C1_3*( y3(i-1) - y1(i-1) ) + G*C1_4*( y4(i-1) ) - y1(i-1) + G*C1_5*( y5(i-1) - y1(i-1) ) + G*C1_6*( y6(i-1) - y1(i-1) ) + G*C1_7*( y7(i-1) - y1(i-1) ) + G*C1_8*( y8(i-1) - y1(i-1) ) + G*C1_9*( y9(i-1) - y1(i-1) ) + G*C1_10*( y10(i-1) - y1(i-1) ) + G*C1_11*( y11(i-1) - y1(i-1) ) + G*C1_12*( y12(i-1) - y1(i-1) ) + G*C1_13*( y13(i-1) - y1(i-1) ) + G*C1_14*( y14(i-1) - y1(i-1) ) + G*C1_15*( y15(i-1) - y1(i-1) ) + G*C1_16*( y16(i-1) - y1(i-1) ) + G*C1_17*( y17(i-1) - y1(i-1) ) + G*C1_18*( y18(i-1) - y1(i-1) ) + G*C1_19*( y19(i-1) - y1(i-1) ) + G*C1_20*( y20(i-1) - y1(i-1) ) + G*C1_21*( y21(i-1) - y1(i-1) ) + G*C1_22*( y22(i-1) - y1(i-1) ) + G*C1_23*( y23(i-1) - y1(i-1) ) )*dt;
    %Oscillator 2
    x2(i) = x2(i-1) + ( ( a2 - x2(i-1)^2 - y2(i-1)^2 )*x2(i-1) - omega2*y2(i-1) + G*C2_1*( x1(i-1) - x2(i-1) ) + G*C2_3*( x3(i-1) - x2(i-1) ) + G*C2_4*( x4(i-1) ) - x2(i-1) + G*C2_5*( x5(i-1) - x2(i-1) ) + G*C2_6*( x6(i-1) - x2(i-1) ) + G*C2_7*( x7(i-1) - x2(i-1) ) + G*C2_8*( x8(i-1) - x2(i-1) ) + G*C2_9*( x9(i-1) - x2(i-1) ) + G*C2_10*( x10(i-1) - x2(i-1) ) + G*C2_11*( x11(i-1) - x2(i-1) ) + G*C2_12*( x12(i-1) - x2(i-1) ) + G*C2_13*( x13(i-1) - x2(i-1) ) + G*C2_14*( x14(i-1) - x2(i-1) ) + G*C2_15*( x15(i-1) - x2(i-1) ) + G*C2_16*( x16(i-1) - x2(i-1) ) + G*C2_17*( x17(i-1) - x2(i-1) ) + G*C2_18*( x18(i-1) - x2(i-1) ) + G*C2_19*( x19(i-1) - x2(i-1) ) + G*C2_20*( x20(i-1) - x2(i-1) ) + G*C2_21*( x21(i-1) - x2(i-1) ) + G*C2_22*( x22(i-1) - x2(i-1) ) + G*C2_23*( x23(i-1) - x2(i-1) ) )*dt;
    y2(i) = y2(i-1) + ( ( a2 - x2(i-1)^2 - y2(i-1)^2 )*y2(i-1) + omega2*x2(i-1) + G*C2_1*( y1(i-1) - y2(i-1) ) + G*C2_3*( y3(i-1) - y2(i-1) ) + G*C2_4*( y4(i-1) ) - y2(i-1) + G*C2_5*( y5(i-1) - y2(i-1) ) + G*C2_6*( y6(i-1) - y2(i-1) ) + G*C2_7*( y7(i-1) - y2(i-1) ) + G*C2_8*( y8(i-1) - y2(i-1) ) + G*C2_9*( y9(i-1) - y2(i-1) ) + G*C2_10*( y10(i-1) - y2(i-1) ) + G*C2_11*( y11(i-1) - y2(i-1) ) + G*C2_12*( y12(i-1) - y2(i-1) ) + G*C2_13*( y13(i-1) - y2(i-1) ) + G*C2_14*( y14(i-1) - y2(i-1) ) + G*C2_15*( y15(i-1) - y2(i-1) ) + G*C2_16*( y16(i-1) - y2(i-1) ) + G*C2_17*( y17(i-1) - y2(i-1) ) + G*C2_18*( y18(i-1) - y2(i-1) ) + G*C2_19*( y19(i-1) - y2(i-1) ) + G*C2_20*( y20(i-1) - y2(i-1) ) + G*C2_21*( y21(i-1) - y2(i-1) ) + G*C2_22*( y22(i-1) - x2(i-1) ) + G*C2_23*( y23(i-1) - x2(i-1) ) )*dt;
    %Oscillator 3
    x3(i) = x3(i-1) + ( ( a3 - x3(i-1)^2 - y3(i-1)^2 )*x3(i-1) - omega3*y3(i-1) + G*C3_1*( x1(i-1) - x3(i-1) ) + G*C3_2*( x2(i-1) - x3(i-1) ) + G*C3_4*( x4(i-1) ) - x3(i-1) + G*C3_5*( x5(i-1) - x3(i-1) ) + G*C3_6*( x6(i-1) - x3(i-1) ) + G*C3_7*( x7(i-1) - x3(i-1) ) + G*C3_8*( x8(i-1) - x3(i-1) ) + G*C3_9*( x9(i-1) - x3(i-1) ) + G*C3_10*( x10(i-1) - x3(i-1) ) + G*C3_11*( x11(i-1) - x3(i-1) ) + G*C3_12*( x12(i-1) - x3(i-1) ) + G*C3_13*( x13(i-1) - x3(i-1) ) + G*C3_14*( x14(i-1) - x3(i-1) ) + G*C3_15*( x15(i-1) - x3(i-1) ) + G*C3_16*( x16(i-1) - x3(i-1) ) + G*C3_17*( x17(i-1) - x3(i-1) ) + G*C3_18*( x18(i-1) - x3(i-1) ) + G*C3_19*( x19(i-1) - x3(i-1) ) + G*C3_20*( x20(i-1) - x3(i-1) ) + G*C3_21*( x21(i-1) - x3(i-1) ) + G*C3_22*( x22(i-1) - x3(i-1) ) + G*C3_23*( x23(i-1) - x3(i-1) ) )*dt;
    y3(i) = y3(i-1) + ( ( a3 - x3(i-1)^2 - y3(i-1)^2 )*y3(i-1) + omega3*x3(i-1) + G*C3_1*( y1(i-1) - y3(i-1) ) + G*C3_2*( y2(i-1) - y3(i-1) ) + G*C3_4*( y4(i-1) ) - y3(i-1) + G*C3_5*( y5(i-1) - y3(i-1) ) + G*C3_6*( y6(i-1) - y3(i-1) ) + G*C3_7*( y7(i-1) - y3(i-1) ) + G*C3_8*( y8(i-1) - y3(i-1) ) + G*C3_9*( y9(i-1) - y3(i-1) ) + G*C3_10*( y10(i-1) - y3(i-1) ) + G*C3_11*( y11(i-1) - y3(i-1) ) + G*C3_12*( y12(i-1) - y3(i-1) ) + G*C3_13*( y13(i-1) - y3(i-1) ) + G*C3_14*( y14(i-1) - y3(i-1) ) + G*C3_15*( y15(i-1) - y3(i-1) ) + G*C3_16*( y16(i-1) - y3(i-1) ) + G*C3_17*( y17(i-1) - y3(i-1) ) + G*C3_18*( y18(i-1) - y3(i-1) ) + G*C3_19*( y19(i-1) - y3(i-1) ) + G*C3_20*( y20(i-1) - y3(i-1) ) + G*C3_21*( y21(i-1) - y3(i-1) ) + G*C3_22*( y22(i-1) - y3(i-1) ) + G*C3_23*( y23(i-1) - y3(i-1) ) )*dt;
    %Oscillator 4
    x4(i) = x4(i-1) + ( ( a4 - x4(i-1)^2 - y4(i-1)^2 )*x4(i-1) - omega4*y4(i-1) + G*C4_1*( x1(i-1) - x4(i-1) ) + G*C4_2*( x2(i-1) - x4(i-1) ) + G*C4_3*( x3(i-1) ) - x4(i-1) + G*C4_5*( x5(i-1) - x4(i-1) ) + G*C4_6*( x6(i-1) - x4(i-1) ) + G*C4_7*( x7(i-1) - x4(i-1) ) + G*C4_8*( x8(i-1) - x4(i-1) ) + G*C4_9*( x9(i-1) - x4(i-1) ) + G*C4_10*( x10(i-1) - x4(i-1) ) + G*C4_11*( x11(i-1) - x4(i-1) ) + G*C4_12*( x12(i-1) - x4(i-1) ) + G*C4_13*( x13(i-1) - x4(i-1) ) + G*C4_14*( x14(i-1) - x4(i-1) ) + G*C4_15*( x15(i-1) - x4(i-1) ) + G*C4_16*( x16(i-1) - x4(i-1) ) + G*C4_17*( x17(i-1) - x4(i-1) ) + G*C4_18*( x18(i-1) - x4(i-1) ) + G*C4_19*( x19(i-1) - x4(i-1) ) + G*C4_20*( x20(i-1) - x4(i-1) ) + G*C4_21*( x21(i-1) - x4(i-1) ) + G*C4_22*( x22(i-1) - x4(i-1) ) + G*C4_23*( x23(i-1) - x4(i-1) ) )*dt;
    y4(i) = y4(i-1) + ( ( a4 - x4(i-1)^2 - y4(i-1)^2 )*y4(i-1) + omega4*x4(i-1) + G*C4_1*( y1(i-1) - y4(i-1) ) + G*C4_2*( y2(i-1) - y4(i-1) ) + G*C4_3*( y3(i-1) ) - y4(i-1) + G*C4_5*( y5(i-1) - y4(i-1) ) + G*C4_6*( y6(i-1) - y4(i-1) ) + G*C4_7*( y7(i-1) - y4(i-1) ) + G*C4_8*( y8(i-1) - y4(i-1) ) + G*C4_9*( y9(i-1) - y4(i-1) ) + G*C4_10*( y10(i-1) - y4(i-1) ) + G*C4_11*( y11(i-1) - y4(i-1) ) + G*C4_12*( y12(i-1) - y4(i-1) ) + G*C4_13*( y13(i-1) - y4(i-1) ) + G*C4_14*( y14(i-1) - y4(i-1) ) + G*C4_15*( y15(i-1) - y4(i-1) ) + G*C4_16*( y16(i-1) - y4(i-1) ) + G*C4_17*( y17(i-1) - y4(i-1) ) + G*C4_18*( y18(i-1) - y4(i-1) ) + G*C4_19*( y19(i-1) - y4(i-1) ) + G*C4_20*( y20(i-1) - y4(i-1) ) + G*C4_21*( y21(i-1) - y4(i-1) ) + G*C4_22*( y22(i-1) - y4(i-1) ) + G*C4_23*( y23(i-1) - y4(i-1) ) )*dt;
    %Oscillator 5
    x5(i) = x5(i-1) + ( ( a5 - x5(i-1)^2 - y5(i-1)^2 )*x5(i-1) - omega5*y5(i-1) + G*C5_1*( x1(i-1) - x5(i-1) ) + G*C5_2*( x2(i-1) - x5(i-1) ) + G*C5_3*( x3(i-1) ) - x5(i-1) + G*C5_4*( x4(i-1) - x5(i-1) ) + G*C5_6*( x6(i-1) - x5(i-1) ) + G*C5_7*( x7(i-1) - x5(i-1) ) + G*C5_8*( x8(i-1) - x5(i-1) ) + G*C5_9*( x9(i-1) - x5(i-1) ) + G*C5_10*( x10(i-1) - x5(i-1) ) + G*C5_11*( x11(i-1) - x5(i-1) ) + G*C5_12*( x12(i-1) - x5(i-1) ) + G*C5_13*( x13(i-1) - x5(i-1) ) + G*C5_14*( x14(i-1) - x5(i-1) ) + G*C5_15*( x15(i-1) - x5(i-1) ) + G*C5_16*( x16(i-1) - x5(i-1) ) + G*C5_17*( x17(i-1) - x5(i-1) ) + G*C5_18*( x18(i-1) - x5(i-1) ) + G*C5_19*( x19(i-1) - x5(i-1) ) + G*C5_20*( x20(i-1) - x5(i-1) ) + G*C5_21*( x21(i-1) - x5(i-1) ) + G*C5_22*( x22(i-1) - x5(i-1) ) + G*C5_23*( x23(i-1) - x5(i-1) ) )*dt;
    y5(i) = y5(i-1) + ( ( a5 - x5(i-1)^2 - y5(i-1)^2 )*y5(i-1) + omega5*x5(i-1) + G*C5_1*( y1(i-1) - y5(i-1) ) + G*C5_2*( y2(i-1) - y5(i-1) ) + G*C5_3*( y3(i-1) ) - y5(i-1) + G*C5_4*( y4(i-1) - y5(i-1) ) + G*C5_6*( y6(i-1) - y5(i-1) ) + G*C5_7*( y7(i-1) - y5(i-1) ) + G*C5_8*( y8(i-1) - y5(i-1) ) + G*C5_9*( y9(i-1) - y5(i-1) ) + G*C5_10*( y10(i-1) - y5(i-1) ) + G*C5_11*( y11(i-1) - y5(i-1) ) + G*C5_12*( y12(i-1) - y5(i-1) ) + G*C5_13*( y13(i-1) - y5(i-1) ) + G*C5_14*( y14(i-1) - y5(i-1) ) + G*C5_15*( y15(i-1) - y5(i-1) ) + G*C5_16*( y16(i-1) - y5(i-1) ) + G*C5_17*( y17(i-1) - y5(i-1) ) + G*C5_18*( y18(i-1) - y5(i-1) ) + G*C5_19*( y19(i-1) - y5(i-1) ) + G*C5_20*( y20(i-1) - y5(i-1) ) + G*C5_21*( y21(i-1) - y5(i-1) ) + G*C5_22*( y22(i-1) - y5(i-1) ) + G*C5_23*( y23(i-1) - y5(i-1) ) )*dt;
    %Oscillator 6
    x6(i) = x6(i-1) + ( ( a6 - x6(i-1)^2 - y6(i-1)^2 )*x6(i-1) - omega6*y6(i-1) + G*C6_1*( x1(i-1) - x6(i-1) ) + G*C6_2*( x2(i-1) - x6(i-1) ) + G*C6_3*( x3(i-1) ) - x6(i-1) + G*C6_4*( x4(i-1) - x6(i-1) ) + G*C6_5*( x5(i-1) - x6(i-1) ) + G*C6_7*( x7(i-1) - x6(i-1) ) + G*C6_8*( x8(i-1) - x6(i-1) ) + G*C6_9*( x9(i-1) - x6(i-1) ) + G*C6_10*( x10(i-1) - x6(i-1) ) + G*C6_11*( x11(i-1) - x6(i-1) ) + G*C6_12*( x12(i-1) - x6(i-1) ) + G*C6_13*( x13(i-1) - x6(i-1) ) + G*C6_14*( x14(i-1) - x6(i-1) ) + G*C6_15*( x15(i-1) - x6(i-1) ) + G*C6_16*( x16(i-1) - x6(i-1) ) + G*C6_17*( x17(i-1) - x6(i-1) ) + G*C6_18*( x18(i-1) - x6(i-1) ) + G*C6_19*( x19(i-1) - x6(i-1) ) + G*C6_20*( x20(i-1) - x6(i-1) ) + G*C6_21*( x21(i-1) - x6(i-1) ) + G*C6_22*( x22(i-1) - x6(i-1) ) + G*C6_23*( x23(i-1) - x6(i-1) ) )*dt;
    y6(i) = y6(i-1) + ( ( a6 - x6(i-1)^2 - y6(i-1)^2 )*y6(i-1) + omega6*x6(i-1) + G*C6_1*( y1(i-1) - y6(i-1) ) + G*C6_2*( y2(i-1) - y6(i-1) ) + G*C6_3*( y3(i-1) ) - y6(i-1) + G*C6_4*( y4(i-1) - y6(i-1) ) + G*C6_5*( y5(i-1) - y6(i-1) ) + G*C6_7*( y7(i-1) - y6(i-1) ) + G*C6_8*( y8(i-1) - y6(i-1) ) + G*C6_9*( y9(i-1) - y6(i-1) ) + G*C6_10*( y10(i-1) - y6(i-1) ) + G*C6_11*( y11(i-1) - y6(i-1) ) + G*C6_12*( y12(i-1) - y6(i-1) ) + G*C6_13*( y13(i-1) - y6(i-1) ) + G*C6_14*( y14(i-1) - y6(i-1) ) + G*C6_15*( y15(i-1) - y6(i-1) ) + G*C6_16*( y16(i-1) - y6(i-1) ) + G*C6_17*( y17(i-1) - y6(i-1) ) + G*C6_18*( y18(i-1) - y6(i-1) ) + G*C6_19*( y19(i-1) - y6(i-1) ) + G*C6_20*( y20(i-1) - y6(i-1) ) + G*C6_21*( y21(i-1) - y6(i-1) ) + G*C6_22*( y22(i-1) - y6(i-1) ) + G*C6_23*( y23(i-1) - y6(i-1) ) )*dt;
    %Oscillator 7
    x7(i) = x7(i-1) + ( ( a7 - x7(i-1)^2 - y7(i-1)^2 )*x7(i-1) - omega7*y7(i-1) + G*C7_1*( x1(i-1) - x7(i-1) ) + G*C7_2*( x2(i-1) - x7(i-1) ) + G*C7_3*( x3(i-1) ) - x7(i-1) + G*C7_4*( x4(i-1) - x7(i-1) ) + G*C7_5*( x5(i-1) - x7(i-1) ) + G*C7_6*( x6(i-1) - x7(i-1) ) + G*C7_8*( x8(i-1) - x7(i-1) ) + G*C7_9*( x9(i-1) - x7(i-1) ) + G*C7_10*( x10(i-1) - x7(i-1) ) + G*C7_11*( x11(i-1) - x7(i-1) ) + G*C7_12*( x12(i-1) - x7(i-1) ) + G*C7_13*( x13(i-1) - x7(i-1) ) + G*C7_14*( x14(i-1) - x7(i-1) ) + G*C7_15*( x15(i-1) - x7(i-1) ) + G*C7_16*( x16(i-1) - x7(i-1) ) + G*C7_17*( x17(i-1) - x7(i-1) ) + G*C7_18*( x18(i-1) - x7(i-1) ) + G*C7_19*( x19(i-1) - x7(i-1) ) + G*C7_20*( x20(i-1) - x7(i-1) ) + G*C7_21*( x21(i-1) - x7(i-1) ) + G*C7_22*( x22(i-1) - x7(i-1) ) + G*C7_23*( x23(i-1) - x7(i-1) ) )*dt;
    y7(i) = y7(i-1) + ( ( a7 - x7(i-1)^2 - y7(i-1)^2 )*y7(i-1) + omega7*x7(i-1) + G*C7_1*( y1(i-1) - y7(i-1) ) + G*C7_2*( y2(i-1) - y7(i-1) ) + G*C7_3*( y3(i-1) ) - y7(i-1) + G*C7_4*( y4(i-1) - y7(i-1) ) + G*C7_5*( y5(i-1) - y7(i-1) ) + G*C7_6*( y6(i-1) - y7(i-1) ) + G*C7_8*( y8(i-1) - y7(i-1) ) + G*C7_9*( y9(i-1) - y7(i-1) ) + G*C7_10*( y10(i-1) - y7(i-1) ) + G*C7_11*( y11(i-1) - y7(i-1) ) + G*C7_12*( y12(i-1) - y7(i-1) ) + G*C7_13*( y13(i-1) - y7(i-1) ) + G*C7_14*( y14(i-1) - y7(i-1) ) + G*C7_15*( y15(i-1) - y7(i-1) ) + G*C7_16*( y16(i-1) - y7(i-1) ) + G*C7_17*( y17(i-1) - y7(i-1) ) + G*C7_18*( y18(i-1) - y7(i-1) ) + G*C7_19*( y19(i-1) - y7(i-1) ) + G*C7_20*( y20(i-1) - y7(i-1) ) + G*C7_21*( y21(i-1) - y7(i-1) ) + G*C7_22*( y22(i-1) - y7(i-1) ) + G*C7_23*( y23(i-1) - y7(i-1) ) )*dt;
    %Oscillator 8
    x8(i) = x8(i-1) + ( ( a8 - x8(i-1)^2 - y8(i-1)^2 )*x8(i-1) - omega8*y8(i-1) + G*C8_1*( x1(i-1) - x8(i-1) ) + G*C8_2*( x2(i-1) - x8(i-1) ) + G*C8_3*( x3(i-1) ) - x8(i-1) + G*C8_4*( x4(i-1) - x8(i-1) ) + G*C8_5*( x5(i-1) - x8(i-1) ) + G*C8_6*( x6(i-1) - x8(i-1) ) + G*C8_7*( x7(i-1) - x8(i-1) ) + G*C8_9*( x9(i-1) - x8(i-1) ) + G*C8_10*( x10(i-1) - x8(i-1) ) + G*C8_11*( x11(i-1) - x8(i-1) ) + G*C8_12*( x12(i-1) - x8(i-1) ) + G*C8_13*( x13(i-1) - x8(i-1) ) + G*C8_14*( x14(i-1) - x8(i-1) ) + G*C8_15*( x15(i-1) - x8(i-1) ) + G*C8_16*( x16(i-1) - x8(i-1) ) + G*C8_17*( x17(i-1) - x8(i-1) ) + G*C8_18*( x18(i-1) - x8(i-1) ) + G*C8_19*( x19(i-1) - x8(i-1) ) + G*C8_20*( x20(i-1) - x8(i-1) ) + G*C8_21*( x21(i-1) - x8(i-1) ) + G*C8_22*( x22(i-1) - x8(i-1) ) + G*C8_23*( x23(i-1) - x8(i-1) ) )*dt;
    y8(i) = y8(i-1) + ( ( a8 - x8(i-1)^2 - y8(i-1)^2 )*y8(i-1) + omega8*x8(i-1) + G*C8_1*( y1(i-1) - y8(i-1) ) + G*C8_2*( y2(i-1) - y8(i-1) ) + G*C8_3*( y3(i-1) ) - y8(i-1) + G*C8_4*( y4(i-1) - y8(i-1) ) + G*C8_5*( y5(i-1) - y8(i-1) ) + G*C8_6*( y6(i-1) - y8(i-1) ) + G*C8_7*( y7(i-1) - y8(i-1) ) + G*C8_9*( y9(i-1) - y8(i-1) ) + G*C8_10*( y10(i-1) - y8(i-1) ) + G*C8_11*( y11(i-1) - y8(i-1) ) + G*C8_12*( y12(i-1) - y8(i-1) ) + G*C8_13*( y13(i-1) - y8(i-1) ) + G*C8_14*( y14(i-1) - y8(i-1) ) + G*C8_15*( y15(i-1) - y8(i-1) ) + G*C8_16*( y16(i-1) - y8(i-1) ) + G*C8_17*( y17(i-1) - y8(i-1) ) + G*C8_18*( y18(i-1) - y8(i-1) ) + G*C8_19*( y19(i-1) - y8(i-1) ) + G*C8_20*( y20(i-1) - y8(i-1) ) + G*C8_21*( y21(i-1) - y8(i-1) ) + G*C8_22*( y22(i-1) - y8(i-1) ) + G*C8_23*( y23(i-1) - y8(i-1) ) )*dt;
    %Oscillator 9
    x9(i) = x9(i-1) + ( ( a9 - x9(i-1)^2 - y9(i-1)^2 )*x9(i-1) - omega9*y9(i-1) + G*C9_1*( x1(i-1) - x9(i-1) ) + G*C9_2*( x2(i-1) - x9(i-1) ) + G*C9_3*( x3(i-1) ) - x9(i-1) + G*C9_4*( x4(i-1) - x9(i-1) ) + G*C9_5*( x5(i-1) - x9(i-1) ) + G*C9_6*( x6(i-1) - x9(i-1) ) + G*C9_7*( x7(i-1) - x9(i-1) ) + G*C9_8*( x8(i-1) - x9(i-1) ) + G*C9_10*( x10(i-1) - x9(i-1) ) + G*C9_11*( x11(i-1) - x9(i-1) ) + G*C9_12*( x12(i-1) - x9(i-1) ) + G*C9_13*( x13(i-1) - x9(i-1) ) + G*C9_14*( x14(i-1) - x9(i-1) ) + G*C9_15*( x15(i-1) - x9(i-1) ) + G*C9_16*( x16(i-1) - x9(i-1) ) + G*C9_17*( x17(i-1) - x9(i-1) ) + G*C9_18*( x18(i-1) - x9(i-1) ) + G*C9_19*( x19(i-1) - x9(i-1) ) + G*C9_20*( x20(i-1) - x9(i-1) ) + G*C9_21*( x21(i-1) - x9(i-1) ) + G*C9_22*( x22(i-1) - x9(i-1) ) + G*C9_23*( x23(i-1) - x9(i-1) ) )*dt;
    y9(i) = y9(i-1) + ( ( a9 - x9(i-1)^2 - y9(i-1)^2 )*y9(i-1) + omega9*x9(i-1) + G*C9_1*( y1(i-1) - y9(i-1) ) + G*C9_2*( y2(i-1) - y9(i-1) ) + G*C9_3*( y3(i-1) ) - y9(i-1) + G*C9_4*( y4(i-1) - y9(i-1) ) + G*C9_5*( y5(i-1) - y9(i-1) ) + G*C9_6*( y6(i-1) - y9(i-1) ) + G*C9_7*( y7(i-1) - y9(i-1) ) + G*C9_8*( y8(i-1) - y9(i-1) ) + G*C9_10*( y10(i-1) - y9(i-1) ) + G*C9_11*( y11(i-1) - y9(i-1) ) + G*C9_12*( y12(i-1) - y9(i-1) ) + G*C9_13*( y13(i-1) - y9(i-1) ) + G*C9_14*( y14(i-1) - y9(i-1) ) + G*C9_15*( y15(i-1) - y9(i-1) ) + G*C9_16*( y16(i-1) - y9(i-1) ) + G*C9_17*( y17(i-1) - y9(i-1) ) + G*C9_18*( y18(i-1) - y9(i-1) ) + G*C9_19*( y19(i-1) - y9(i-1) ) + G*C9_20*( y20(i-1) - y9(i-1) ) + G*C9_21*( y21(i-1) - y9(i-1) ) + G*C9_22*( y22(i-1) - y9(i-1) ) + G*C9_23*( y23(i-1) - y9(i-1) ) )*dt;
    %Oscillator 10
    x10(i) = x10(i-1) + ( ( a10 - x10(i-1)^2 - y10(i-1)^2 )*x10(i-1) - omega10*y10(i-1) + G*C10_1*( x1(i-1) - x10(i-1) ) + G*C10_2*( x2(i-1) - x10(i-1) ) + G*C10_3*( x3(i-1) - x10(i-1) ) + G*C10_4*( x4(i-1) - x10(i-1) ) + G*C10_5*( x5(i-1) - x10(i-1) ) + G*C10_6*( x6(i-1) - x10(i-1) ) + G*C10_7*( x7(i-1) - x10(i-1) ) + G*C10_8*( x8(i-1) - x10(i-1) ) + G*C10_9*( x9(i-1) - x10(i-1) ) + G*C10_11*( x11(i-1) - x10(i-1) ) + G*C10_12*( x12(i-1) - x10(i-1) ) + G*C10_13*( x13(i-1) - x10(i-1) ) + G*C10_14*( x14(i-1) - x10(i-1) ) + G*C10_15*( x15(i-1) - x10(i-1) ) + G*C10_16*( x16(i-1) - x10(i-1) ) + G*C10_17*( x17(i-1) - x10(i-1) ) + G*C10_18*( x18(i-1) - x10(i-1) ) + G*C10_19*( x19(i-1) - x10(i-1) ) + G*C10_20*( x20(i-1) - x10(i-1) ) + G*C10_21*( x21(i-1) - x10(i-1) ) + G*C10_22*( x22(i-1) - x10(i-1) ) + G*C10_23*( x23(i-1) - x10(i-1) ) )*dt;
    y10(i) = y10(i-1) + ( ( a10 - x10(i-1)^2 - y10(i-1)^2 )*y10(i-1) + omega10*x10(i-1) + G*C10_1*( y1(i-1) - y10(i-1) ) + G*C10_2*( y2(i-1) - y10(i-1) ) + G*C10_3*( y3(i-1) - y10(i-1) ) + G*C10_4*( y4(i-1) - y10(i-1) ) + G*C10_5*( y5(i-1) - y10(i-1) ) + G*C10_6*( y6(i-1) - y10(i-1) ) + G*C10_7*( y7(i-1) - y10(i-1) ) + G*C10_8*( y8(i-1) - y10(i-1) ) + G*C10_9*( y9(i-1) - y10(i-1) ) + G*C10_11*( y11(i-1) - y10(i-1) ) + G*C10_12*( y12(i-1) - y10(i-1) ) + G*C10_13*( y13(i-1) - y10(i-1) ) + G*C10_14*( y14(i-1) - y10(i-1) ) + G*C10_15*( y15(i-1) - y10(i-1) ) + G*C10_16*( y16(i-1) - y10(i-1) ) + G*C10_17*( y17(i-1) - y10(i-1) ) + G*C10_18*( y18(i-1) - y10(i-1) ) + G*C10_19*( y19(i-1) - y10(i-1) ) + G*C10_20*( y20(i-1) - y10(i-1) ) + G*C10_21*( y21(i-1) - y10(i-1) ) + G*C10_22*( y22(i-1) - y10(i-1) ) + G*C10_23*( y23(i-1) - y10(i-1) ) )*dt;
    %Oscillator 11
    x11(i) = x11(i-1) + ( ( a11 - x11(i-1)^2 - y11(i-1)^2 )*x11(i-1) - omega11*y11(i-1) + G*C11_1*( x1(i-1) - x11(i-1) ) + G*C11_2*( x2(i-1) - x11(i-1) ) + G*C11_3*( x3(i-1) - x11(i-1) ) + G*C11_4*( x4(i-1) - x11(i-1) ) + G*C11_5*( x5(i-1) - x11(i-1) ) + G*C11_6*( x6(i-1) - x11(i-1) ) + G*C11_7*( x7(i-1) - x11(i-1) ) + G*C11_8*( x8(i-1) - x11(i-1) ) + G*C11_9*( x9(i-1) - x11(i-1) ) + G*C11_10*( x10(i-1) - x11(i-1) ) + G*C11_12*( x12(i-1) - x11(i-1) ) + G*C11_13*( x13(i-1) - x11(i-1) ) + G*C11_14*( x14(i-1) - x11(i-1) ) + G*C11_15*( x15(i-1) - x11(i-1) ) + G*C11_16*( x16(i-1) - x11(i-1) ) + G*C11_17*( x17(i-1) - x11(i-1) ) + G*C11_18*( x18(i-1) - x11(i-1) ) + G*C11_19*( x19(i-1) - x11(i-1) ) + G*C11_20*( x20(i-1) - x11(i-1) ) + G*C11_21*( x21(i-1) - x11(i-1) ) + G*C11_22*( x22(i-1) - x11(i-1) ) + G*C11_23*( x23(i-1) - x11(i-1) ) )*dt;
    y11(i) = y11(i-1) + ( ( a11 - x11(i-1)^2 - y11(i-1)^2 )*y11(i-1) + omega11*x11(i-1) + G*C11_1*( y1(i-1) - y11(i-1) ) + G*C11_2*( y2(i-1) - y11(i-1) ) + G*C11_3*( y3(i-1) - y11(i-1) ) + G*C11_4*( y4(i-1) - y11(i-1) ) + G*C11_5*( y5(i-1) - y11(i-1) ) + G*C11_6*( y6(i-1) - y11(i-1) ) + G*C11_7*( y7(i-1) - y11(i-1) ) + G*C11_8*( y8(i-1) - y11(i-1) ) + G*C11_9*( y9(i-1) - y11(i-1) ) + G*C11_10*( y10(i-1) - y11(i-1) ) + G*C11_12*( y12(i-1) - y11(i-1) ) + G*C11_13*( y13(i-1) - y11(i-1) ) + G*C11_14*( y14(i-1) - y11(i-1) ) + G*C11_15*( y15(i-1) - y11(i-1) ) + G*C11_16*( y16(i-1) - y11(i-1) ) + G*C11_17*( y17(i-1) - y11(i-1) ) + G*C11_18*( y18(i-1) - y11(i-1) ) + G*C11_19*( y19(i-1) - y11(i-1) ) + G*C11_20*( y20(i-1) - y11(i-1) ) + G*C11_21*( y21(i-1) - y11(i-1) ) + G*C11_22*( y22(i-1) - y11(i-1) ) + G*C11_23*( y23(i-1) - y11(i-1) ) )*dt;
    %Oscillator 12
    x12(i) = x12(i-1) + ( ( a12 - x12(i-1)^2 - y12(i-1)^2 )*x12(i-1) - omega12*y12(i-1) + G*C12_1*( x1(i-1) - x12(i-1) ) + G*C12_2*( x2(i-1) - x12(i-1) ) + G*C12_3*( x3(i-1) - x12(i-1) ) + G*C12_4*( x4(i-1) - x12(i-1) ) + G*C12_5*( x5(i-1) - x12(i-1) ) + G*C12_6*( x6(i-1) - x12(i-1) ) + G*C12_7*( x7(i-1) - x12(i-1) ) + G*C12_8*( x8(i-1) - x12(i-1) ) + G*C12_9*( x9(i-1) - x12(i-1) ) + G*C12_10*( x10(i-1) - x12(i-1) ) + G*C12_11*( x11(i-1) - x12(i-1) ) + G*C12_13*( x13(i-1) - x12(i-1) ) + G*C12_14*( x14(i-1) - x12(i-1) ) + G*C12_15*( x15(i-1) - x12(i-1) ) + G*C12_16*( x16(i-1) - x12(i-1) ) + G*C12_17*( x17(i-1) - x12(i-1) ) + G*C12_18*( x18(i-1) - x12(i-1) ) + G*C12_19*( x19(i-1) - x12(i-1) ) + G*C12_20*( x20(i-1) - x12(i-1) ) + G*C12_21*( x21(i-1) - x12(i-1) ) + G*C12_22*( x22(i-1) - x12(i-1) ) + G*C12_23*( x23(i-1) - x12(i-1) ) )*dt;
    y12(i) = y12(i-1) + ( ( a12 - x12(i-1)^2 - y12(i-1)^2 )*y12(i-1) + omega12*x12(i-1) + G*C12_1*( y1(i-1) - y12(i-1) ) + G*C12_2*( y2(i-1) - y12(i-1) ) + G*C12_3*( y3(i-1) - y12(i-1) ) + G*C12_4*( y4(i-1) - y12(i-1) ) + G*C12_5*( y5(i-1) - y12(i-1) ) + G*C12_6*( y6(i-1) - y12(i-1) ) + G*C12_7*( y7(i-1) - y12(i-1) ) + G*C12_8*( y8(i-1) - y12(i-1) ) + G*C12_9*( y9(i-1) - y12(i-1) ) + G*C12_10*( y10(i-1) - y12(i-1) ) + G*C12_11*( y11(i-1) - y12(i-1) ) + G*C12_13*( y13(i-1) - y12(i-1) ) + G*C12_14*( y14(i-1) - y12(i-1) ) + G*C12_15*( y15(i-1) - y12(i-1) ) + G*C12_16*( y16(i-1) - y12(i-1) ) + G*C12_17*( y17(i-1) - y12(i-1) ) + G*C12_18*( y18(i-1) - y12(i-1) ) + G*C12_19*( y19(i-1) - y12(i-1) ) + G*C12_20*( y20(i-1) - y12(i-1) ) + G*C12_21*( y21(i-1) - y12(i-1) ) + G*C12_22*( y22(i-1) - y12(i-1) ) + G*C12_23*( y23(i-1) - y12(i-1) ) )*dt;
    %Oscillator 13
    x13(i) = x13(i-1) + ( ( a13 - x13(i-1)^2 - y13(i-1)^2 )*x13(i-1) - omega13*y13(i-1) + G*C13_1*( x1(i-1) - x13(i-1) ) + G*C13_2*( x2(i-1) - x13(i-1) ) + G*C13_3*( x3(i-1) - x13(i-1) ) + G*C13_4*( x4(i-1) - x13(i-1) ) + G*C13_5*( x5(i-1) - x13(i-1) ) + G*C13_6*( x6(i-1) - x13(i-1) ) + G*C13_7*( x7(i-1) - x13(i-1) ) + G*C13_8*( x8(i-1) - x13(i-1) ) + G*C13_9*( x9(i-1) - x13(i-1) ) + G*C13_10*( x10(i-1) - x13(i-1) ) + G*C13_11*( x11(i-1) - x13(i-1) ) + G*C13_12*( x12(i-1) - x13(i-1) ) + G*C13_14*( x14(i-1) - x13(i-1) ) + G*C13_15*( x15(i-1) - x13(i-1) ) + G*C13_16*( x16(i-1) - x13(i-1) ) + G*C13_17*( x17(i-1) - x13(i-1) ) + G*C13_18*( x18(i-1) - x13(i-1) ) + G*C13_19*( x19(i-1) - x13(i-1) ) + G*C13_20*( x20(i-1) - x13(i-1) ) + G*C13_21*( x21(i-1) - x13(i-1) ) + G*C13_22*( x22(i-1) - x13(i-1) ) + G*C13_23*( x23(i-1) - x13(i-1) ) )*dt;
    y13(i) = y13(i-1) + ( ( a13 - x13(i-1)^2 - y13(i-1)^2 )*y13(i-1) + omega13*x13(i-1) + G*C13_1*( y1(i-1) - y13(i-1) ) + G*C13_2*( y2(i-1) - y13(i-1) ) + G*C13_3*( y3(i-1) - y13(i-1) ) + G*C13_4*( y4(i-1) - y13(i-1) ) + G*C13_5*( y5(i-1) - y13(i-1) ) + G*C13_6*( y6(i-1) - y13(i-1) ) + G*C13_7*( y7(i-1) - y13(i-1) ) + G*C13_8*( y8(i-1) - y13(i-1) ) + G*C13_9*( y9(i-1) - y13(i-1) ) + G*C13_10*( y10(i-1) - y13(i-1) ) + G*C13_11*( y11(i-1) - y13(i-1) ) + G*C13_12*( y12(i-1) - y13(i-1) ) + G*C13_14*( y14(i-1) - y13(i-1) ) + G*C13_15*( y15(i-1) - y13(i-1) ) + G*C13_16*( y16(i-1) - y13(i-1) ) + G*C13_17*( y17(i-1) - y13(i-1) ) + G*C13_18*( y18(i-1) - y13(i-1) ) + G*C13_19*( y19(i-1) - y13(i-1) ) + G*C13_20*( y20(i-1) - y13(i-1) ) + G*C13_21*( y21(i-1) - y13(i-1) ) + G*C13_22*( y22(i-1) - y13(i-1) ) + G*C13_23*( y23(i-1) - y13(i-1) ) )*dt;
    %Oscillator 14
    x14(i) = x14(i-1) + ( ( a14 - x14(i-1)^2 - y14(i-1)^2 )*x14(i-1) - omega14*y14(i-1) + G*C14_1*( x1(i-1) - x14(i-1) ) + G*C14_2*( x2(i-1) - x14(i-1) ) + G*C14_3*( x3(i-1) - x14(i-1) ) + G*C14_4*( x4(i-1) - x14(i-1) ) + G*C14_5*( x5(i-1) - x14(i-1) ) + G*C14_6*( x6(i-1) - x14(i-1) ) + G*C14_7*( x7(i-1) - x14(i-1) ) + G*C14_8*( x8(i-1) - x14(i-1) ) + G*C14_9*( x9(i-1) - x14(i-1) ) + G*C14_10*( x10(i-1) - x14(i-1) ) + G*C14_11*( x11(i-1) - x14(i-1) ) + G*C14_12*( x12(i-1) - x14(i-1) ) + G*C14_13*( x13(i-1) - x14(i-1) ) + G*C14_15*( x15(i-1) - x14(i-1) ) + G*C14_16*( x16(i-1) - x14(i-1) ) + G*C14_17*( x17(i-1) - x14(i-1) ) + G*C14_18*( x18(i-1) - x14(i-1) ) + G*C14_19*( x19(i-1) - x14(i-1) ) + G*C14_20*( x20(i-1) - x14(i-1) ) + G*C14_21*( x21(i-1) - x14(i-1) ) + G*C14_22*( x22(i-1) - x14(i-1) ) + G*C14_23*( x23(i-1) - x14(i-1) ) )*dt;
    y14(i) = y14(i-1) + ( ( a14 - x14(i-1)^2 - y14(i-1)^2 )*y14(i-1) + omega14*x14(i-1) + G*C14_1*( y1(i-1) - y14(i-1) ) + G*C14_2*( y2(i-1) - y14(i-1) ) + G*C14_3*( y3(i-1) - y14(i-1) ) + G*C14_4*( y4(i-1) - y14(i-1) ) + G*C14_5*( y5(i-1) - y14(i-1) ) + G*C14_6*( y6(i-1) - y14(i-1) ) + G*C14_7*( y7(i-1) - y14(i-1) ) + G*C14_8*( y8(i-1) - y14(i-1) ) + G*C14_9*( y9(i-1) - y14(i-1) ) + G*C14_10*( y10(i-1) - y14(i-1) ) + G*C14_11*( y11(i-1) - y14(i-1) ) + G*C14_12*( y12(i-1) - y14(i-1) ) + G*C14_13*( y13(i-1) - y14(i-1) ) + G*C14_15*( y15(i-1) - y14(i-1) ) + G*C14_16*( y16(i-1) - y14(i-1) ) + G*C14_17*( y17(i-1) - y14(i-1) ) + G*C14_18*( y18(i-1) - y14(i-1) ) + G*C14_19*( y19(i-1) - y14(i-1) ) + G*C14_20*( y20(i-1) - y14(i-1) ) + G*C14_21*( y21(i-1) - y14(i-1) ) + G*C14_22*( y22(i-1) - y14(i-1) ) + G*C14_23*( y23(i-1) - y14(i-1) ) )*dt;
    %Oscillator 15
    x15(i) = x15(i-1) + ( ( a15 - x15(i-1)^2 - y15(i-1)^2 )*x15(i-1) - omega15*y15(i-1) + G*C15_1*( x1(i-1) - x15(i-1) ) + G*C15_2*( x2(i-1) - x15(i-1) ) + G*C15_3*( x3(i-1) - x15(i-1) ) + G*C15_4*( x4(i-1) - x15(i-1) ) + G*C15_5*( x5(i-1) - x15(i-1) ) + G*C15_6*( x6(i-1) - x15(i-1) ) + G*C15_7*( x7(i-1) - x15(i-1) ) + G*C15_8*( x8(i-1) - x15(i-1) ) + G*C15_9*( x9(i-1) - x15(i-1) ) + G*C15_10*( x10(i-1) - x15(i-1) ) + G*C15_11*( x11(i-1) - x15(i-1) ) + G*C15_12*( x12(i-1) - x15(i-1) ) + G*C15_13*( x13(i-1) - x15(i-1) ) + G*C15_14*( x14(i-1) - x15(i-1) ) + G*C15_16*( x16(i-1) - x15(i-1) ) + G*C15_17*( x17(i-1) - x15(i-1) ) + G*C15_18*( x18(i-1) - x15(i-1) ) + G*C15_19*( x19(i-1) - x15(i-1) ) + G*C15_20*( x20(i-1) - x15(i-1) ) + G*C15_21*( x21(i-1) - x15(i-1) ) + G*C15_22*( x22(i-1) - x15(i-1) ) + G*C15_23*( x23(i-1) - x15(i-1) ) )*dt;
    y15(i) = y15(i-1) + ( ( a15 - x15(i-1)^2 - y15(i-1)^2 )*y15(i-1) + omega15*x15(i-1) + G*C15_1*( y1(i-1) - y15(i-1) ) + G*C15_2*( y2(i-1) - y15(i-1) ) + G*C15_3*( y3(i-1) - y15(i-1) ) + G*C15_4*( y4(i-1) - y15(i-1) ) + G*C15_5*( y5(i-1) - y15(i-1) ) + G*C15_6*( y6(i-1) - y15(i-1) ) + G*C15_7*( y7(i-1) - y15(i-1) ) + G*C15_8*( y8(i-1) - y15(i-1) ) + G*C15_9*( y9(i-1) - y15(i-1) ) + G*C15_10*( y10(i-1) - y15(i-1) ) + G*C15_11*( y11(i-1) - y15(i-1) ) + G*C15_12*( y12(i-1) - y15(i-1) ) + G*C15_13*( y13(i-1) - y15(i-1) ) + G*C15_14*( y14(i-1) - y15(i-1) ) + G*C15_16*( y16(i-1) - y15(i-1) ) + G*C15_17*( y17(i-1) - y15(i-1) ) + G*C15_18*( y18(i-1) - y15(i-1) ) + G*C15_19*( y19(i-1) - y15(i-1) ) + G*C15_20*( y20(i-1) - y15(i-1) ) + G*C15_21*( y21(i-1) - y15(i-1) ) + G*C15_22*( y22(i-1) - y15(i-1) ) + G*C15_23*( y23(i-1) - y15(i-1) ) )*dt;
    %Oscillator 16
    x16(i) = x16(i-1) + ( ( a16 - x16(i-1)^2 - y16(i-1)^2 )*x16(i-1) - omega16*y16(i-1) + G*C16_1*( x1(i-1) - x16(i-1) ) + G*C16_2*( x2(i-1) - x16(i-1) ) + G*C16_3*( x3(i-1) - x16(i-1) ) + G*C16_4*( x4(i-1) - x16(i-1) ) + G*C16_5*( x5(i-1) - x16(i-1) ) + G*C16_6*( x6(i-1) - x16(i-1) ) + G*C16_7*( x7(i-1) - x16(i-1) ) + G*C16_8*( x8(i-1) - x16(i-1) ) + G*C16_9*( x9(i-1) - x16(i-1) ) + G*C16_10*( x10(i-1) - x16(i-1) ) + G*C16_11*( x11(i-1) - x16(i-1) ) + G*C16_12*( x12(i-1) - x16(i-1) ) + G*C16_13*( x13(i-1) - x16(i-1) ) + G*C16_14*( x14(i-1) - x16(i-1) ) + G*C16_15*( x15(i-1) - x16(i-1) ) + G*C16_17*( x17(i-1) - x16(i-1) ) + G*C16_18*( x18(i-1) - x16(i-1) ) + G*C16_19*( x19(i-1) - x16(i-1) ) + G*C16_20*( x20(i-1) - x16(i-1) ) + G*C16_21*( x21(i-1) - x16(i-1) ) + G*C16_22*( x22(i-1) - x16(i-1) ) + G*C16_23*( x23(i-1) - x16(i-1) ) )*dt;
    y16(i) = y16(i-1) + ( ( a16 - x16(i-1)^2 - y16(i-1)^2 )*y16(i-1) + omega16*x16(i-1) + G*C16_1*( y1(i-1) - y16(i-1) ) + G*C16_2*( y2(i-1) - y16(i-1) ) + G*C16_3*( y3(i-1) - y16(i-1) ) + G*C16_4*( y4(i-1) - y16(i-1) ) + G*C16_5*( y5(i-1) - y16(i-1) ) + G*C16_6*( y6(i-1) - y16(i-1) ) + G*C16_7*( y7(i-1) - y16(i-1) ) + G*C16_8*( y8(i-1) - y16(i-1) ) + G*C16_9*( y9(i-1) - y16(i-1) ) + G*C16_10*( y10(i-1) - y16(i-1) ) + G*C16_11*( y11(i-1) - y16(i-1) ) + G*C16_12*( y12(i-1) - y16(i-1) ) + G*C16_13*( y13(i-1) - y16(i-1) ) + G*C16_14*( y14(i-1) - y16(i-1) ) + G*C16_15*( y15(i-1) - y16(i-1) ) + G*C16_17*( y17(i-1) - y16(i-1) ) + G*C16_18*( y18(i-1) - y16(i-1) ) + G*C16_19*( y19(i-1) - y16(i-1) ) + G*C16_20*( y20(i-1) - y16(i-1) ) + G*C16_21*( y21(i-1) - y16(i-1) ) + G*C16_22*( y22(i-1) - y16(i-1) ) + G*C16_23*( y23(i-1) - y16(i-1) ) )*dt;
    %Oscillator 17
    x17(i) = x17(i-1) + ( ( a17 - x17(i-1)^2 - y17(i-1)^2 )*x17(i-1) - omega17*y17(i-1) + G*C17_1*( x1(i-1) - x17(i-1) ) + G*C17_2*( x2(i-1) - x17(i-1) ) + G*C17_3*( x3(i-1) - x17(i-1) ) + G*C17_4*( x4(i-1) - x17(i-1) ) + G*C17_5*( x5(i-1) - x17(i-1) ) + G*C17_6*( x6(i-1) - x17(i-1) ) + G*C17_7*( x7(i-1) - x17(i-1) ) + G*C17_8*( x8(i-1) - x17(i-1) ) + G*C17_9*( x9(i-1) - x17(i-1) ) + G*C17_10*( x10(i-1) - x17(i-1) ) + G*C17_11*( x11(i-1) - x17(i-1) ) + G*C17_12*( x12(i-1) - x17(i-1) ) + G*C17_13*( x13(i-1) - x17(i-1) ) + G*C17_14*( x14(i-1) - x17(i-1) ) + G*C17_15*( x15(i-1) - x17(i-1) ) + G*C17_16*( x16(i-1) - x17(i-1) ) + G*C17_18*( x18(i-1) - x17(i-1) ) + G*C17_19*( x19(i-1) - x17(i-1) ) + G*C17_20*( x20(i-1) - x17(i-1) ) + G*C17_21*( x21(i-1) - x17(i-1) ) + G*C17_22*( x22(i-1) - x17(i-1) ) + G*C17_23*( x23(i-1) - x17(i-1) ) )*dt;
    y17(i) = y17(i-1) + ( ( a17 - x17(i-1)^2 - y17(i-1)^2 )*y17(i-1) + omega17*x17(i-1) + G*C17_1*( y1(i-1) - y17(i-1) ) + G*C17_2*( y2(i-1) - y17(i-1) ) + G*C17_3*( y3(i-1) - y17(i-1) ) + G*C17_4*( y4(i-1) - y17(i-1) ) + G*C17_5*( y5(i-1) - y17(i-1) ) + G*C17_6*( y6(i-1) - y17(i-1) ) + G*C17_7*( y7(i-1) - y17(i-1) ) + G*C17_8*( y8(i-1) - y17(i-1) ) + G*C17_9*( y9(i-1) - y17(i-1) ) + G*C17_10*( y10(i-1) - y17(i-1) ) + G*C17_11*( y11(i-1) - y17(i-1) ) + G*C17_12*( y12(i-1) - y17(i-1) ) + G*C17_13*( y13(i-1) - y17(i-1) ) + G*C17_14*( y14(i-1) - y17(i-1) ) + G*C17_15*( y15(i-1) - y17(i-1) ) + G*C17_16*( y16(i-1) - y17(i-1) ) + G*C17_18*( y18(i-1) - y17(i-1) ) + G*C17_19*( y19(i-1) - y17(i-1) ) + G*C17_20*( y20(i-1) - y17(i-1) ) + G*C17_21*( y21(i-1) - y17(i-1) ) + G*C17_22*( y22(i-1) - y17(i-1) ) + G*C17_23*( y23(i-1) - y17(i-1) ) )*dt;
    %Oscillator 18
    x18(i) = x18(i-1) + ( ( a18 - x18(i-1)^2 - y18(i-1)^2 )*x18(i-1) - omega18*y18(i-1) + G*C18_1*( x1(i-1) - x18(i-1) ) + G*C18_2*( x2(i-1) - x18(i-1) ) + G*C18_3*( x3(i-1) - x18(i-1) ) + G*C18_4*( x4(i-1) - x18(i-1) ) + G*C18_5*( x5(i-1) - x18(i-1) ) + G*C18_6*( x6(i-1) - x18(i-1) ) + G*C18_7*( x7(i-1) - x18(i-1) ) + G*C18_8*( x8(i-1) - x18(i-1) ) + G*C18_9*( x9(i-1) - x18(i-1) ) + G*C18_10*( x10(i-1) - x18(i-1) ) + G*C18_11*( x11(i-1) - x18(i-1) ) + G*C18_12*( x12(i-1) - x18(i-1) ) + G*C18_13*( x13(i-1) - x18(i-1) ) + G*C18_14*( x14(i-1) - x18(i-1) ) + G*C18_15*( x15(i-1) - x18(i-1) ) + G*C18_16*( x16(i-1) - x18(i-1) ) + G*C18_17*( x17(i-1) - x18(i-1) ) + G*C18_19*( x19(i-1) - x18(i-1) ) + G*C18_20*( x20(i-1) - x18(i-1) ) + G*C18_21*( x21(i-1) - x18(i-1) ) + G*C18_22*( x22(i-1) - x18(i-1) ) + G*C18_23*( x23(i-1) - x18(i-1) ) )*dt;
    y18(i) = y18(i-1) + ( ( a18 - x18(i-1)^2 - y18(i-1)^2 )*y18(i-1) + omega18*x18(i-1) + G*C18_1*( y1(i-1) - y18(i-1) ) + G*C18_2*( y2(i-1) - y18(i-1) ) + G*C18_3*( y3(i-1) - y18(i-1) ) + G*C18_4*( y4(i-1) - y18(i-1) ) + G*C18_5*( y5(i-1) - y18(i-1) ) + G*C18_6*( y6(i-1) - y18(i-1) ) + G*C18_7*( y7(i-1) - y18(i-1) ) + G*C18_8*( y8(i-1) - y18(i-1) ) + G*C18_9*( y9(i-1) - y18(i-1) ) + G*C18_10*( y10(i-1) - y18(i-1) ) + G*C18_11*( y11(i-1) - y18(i-1) ) + G*C18_12*( y12(i-1) - y18(i-1) ) + G*C18_13*( y13(i-1) - y18(i-1) ) + G*C18_14*( y14(i-1) - y18(i-1) ) + G*C18_15*( y15(i-1) - y18(i-1) ) + G*C18_16*( y16(i-1) - y18(i-1) ) + G*C18_17*( y17(i-1) - y18(i-1) ) + G*C18_19*( y19(i-1) - y18(i-1) ) + G*C18_20*( y20(i-1) - y18(i-1) ) + G*C18_21*( y21(i-1) - y18(i-1) ) + G*C18_22*( y22(i-1) - y18(i-1) ) + G*C18_23*( y23(i-1) - y18(i-1) ) )*dt;
    %Oscillator 19
    x19(i) = x19(i-1) + ( ( a19 - x19(i-1)^2 - y19(i-1)^2 )*x19(i-1) - omega19*y19(i-1) + G*C19_1*( x1(i-1) - x19(i-1) ) + G*C19_2*( x2(i-1) - x19(i-1) ) + G*C19_3*( x3(i-1) - x19(i-1) ) + G*C19_4*( x4(i-1) - x19(i-1) ) + G*C19_5*( x5(i-1) - x19(i-1) ) + G*C19_6*( x6(i-1) - x19(i-1) ) + G*C19_7*( x7(i-1) - x19(i-1) ) + G*C19_8*( x8(i-1) - x19(i-1) ) + G*C19_9*( x9(i-1) - x19(i-1) ) + G*C19_10*( x10(i-1) - x19(i-1) ) + G*C19_11*( x11(i-1) - x19(i-1) ) + G*C19_12*( x12(i-1) - x19(i-1) ) + G*C19_13*( x13(i-1) - x19(i-1) ) + G*C19_14*( x14(i-1) - x19(i-1) ) + G*C19_15*( x15(i-1) - x19(i-1) ) + G*C19_16*( x16(i-1) - x19(i-1) ) + G*C19_17*( x17(i-1) - x19(i-1) ) + G*C19_18*( x18(i-1) - x19(i-1) ) + G*C19_20*( x20(i-1) - x19(i-1) ) + G*C19_21*( x21(i-1) - x19(i-1) ) + G*C19_22*( x22(i-1) - x19(i-1) ) + G*C19_23*( x23(i-1) - x19(i-1) ) )*dt;
    y19(i) = y19(i-1) + ( ( a19 - x19(i-1)^2 - y19(i-1)^2 )*y19(i-1) + omega19*x19(i-1) + G*C19_1*( y1(i-1) - y19(i-1) ) + G*C19_2*( y2(i-1) - y19(i-1) ) + G*C19_3*( y3(i-1) - y19(i-1) ) + G*C19_4*( y4(i-1) - y19(i-1) ) + G*C19_5*( y5(i-1) - y19(i-1) ) + G*C19_6*( y6(i-1) - y19(i-1) ) + G*C19_7*( y7(i-1) - y19(i-1) ) + G*C19_8*( y8(i-1) - y19(i-1) ) + G*C19_9*( y9(i-1) - y19(i-1) ) + G*C19_10*( y10(i-1) - y19(i-1) ) + G*C19_11*( y11(i-1) - y19(i-1) ) + G*C19_12*( y12(i-1) - y19(i-1) ) + G*C19_13*( y13(i-1) - y19(i-1) ) + G*C19_14*( y14(i-1) - y19(i-1) ) + G*C19_15*( y15(i-1) - y19(i-1) ) + G*C19_16*( y16(i-1) - y19(i-1) ) + G*C19_17*( y17(i-1) - y19(i-1) ) + G*C19_18*( y18(i-1) - y19(i-1) ) + G*C19_20*( y20(i-1) - y19(i-1) ) + G*C19_21*( y21(i-1) - y19(i-1) ) + G*C19_22*( y22(i-1) - y19(i-1) ) + G*C19_23*( y23(i-1) - y19(i-1) ) )*dt;
    %Oscillator 20
    x20(i) = x20(i-1) + ( ( a20 - x20(i-1)^2 - y20(i-1)^2 )*x20(i-1) - omega20*y20(i-1) + G*C20_1*( x1(i-1) - x20(i-1) ) + G*C20_2*( x2(i-1) - x20(i-1) ) + G*C20_3*( x3(i-1) - x20(i-1) ) + G*C20_4*( x4(i-1) - x20(i-1) ) + G*C20_5*( x5(i-1) - x20(i-1) ) + G*C20_6*( x6(i-1) - x20(i-1) ) + G*C20_7*( x7(i-1) - x20(i-1) ) + G*C20_8*( x8(i-1) - x20(i-1) ) + G*C20_9*( x9(i-1) - x20(i-1) ) + G*C20_10*( x10(i-1) - x20(i-1) ) + G*C20_11*( x11(i-1) - x20(i-1) ) + G*C20_12*( x12(i-1) - x20(i-1) ) + G*C20_13*( x13(i-1) - x20(i-1) ) + G*C20_14*( x14(i-1) - x20(i-1) ) + G*C20_15*( x15(i-1) - x20(i-1) ) + G*C20_16*( x16(i-1) - x20(i-1) ) + G*C20_17*( x17(i-1) - x20(i-1) ) + G*C20_18*( x18(i-1) - x20(i-1) ) + G*C20_19*( x19(i-1) - x20(i-1) ) + G*C20_21*( x21(i-1) - x20(i-1) ) + G*C20_22*( x22(i-1) - x20(i-1) ) + G*C20_23*( x23(i-1) - x20(i-1) ) )*dt;
    y20(i) = y20(i-1) + ( ( a20 - x20(i-1)^2 - y20(i-1)^2 )*y20(i-1) + omega20*x20(i-1) + G*C20_1*( y1(i-1) - y20(i-1) ) + G*C20_2*( y2(i-1) - y20(i-1) ) + G*C20_3*( y3(i-1) - y20(i-1) ) + G*C20_4*( y4(i-1) - y20(i-1) ) + G*C20_5*( y5(i-1) - y20(i-1) ) + G*C20_6*( y6(i-1) - y20(i-1) ) + G*C20_7*( y7(i-1) - y20(i-1) ) + G*C20_8*( y8(i-1) - y20(i-1) ) + G*C20_9*( y9(i-1) - y20(i-1) ) + G*C20_10*( y10(i-1) - y20(i-1) ) + G*C20_11*( y11(i-1) - y20(i-1) ) + G*C20_12*( y12(i-1) - y20(i-1) ) + G*C20_13*( y13(i-1) - y20(i-1) ) + G*C20_14*( y14(i-1) - y20(i-1) ) + G*C20_15*( y15(i-1) - y20(i-1) ) + G*C20_16*( y16(i-1) - y20(i-1) ) + G*C20_17*( y17(i-1) - y20(i-1) ) + G*C20_18*( y18(i-1) - y20(i-1) ) + G*C20_19*( y19(i-1) - y20(i-1) ) + G*C20_21*( y21(i-1) - y20(i-1) ) + G*C20_22*( y22(i-1) - y20(i-1) ) + G*C20_23*( y23(i-1) - y20(i-1) ) )*dt;
    %Oscillator 21
    x21(i) = x21(i-1) + ( ( a21 - x21(i-1)^2 - y21(i-1)^2 )*x21(i-1) - omega21*y21(i-1) + G*C21_1*( x1(i-1) - x21(i-1) ) + G*C21_2*( x2(i-1) - x21(i-1) ) + G*C21_3*( x3(i-1) - x21(i-1) ) + G*C21_4*( x4(i-1) - x21(i-1) ) + G*C21_5*( x5(i-1) - x21(i-1) ) + G*C21_6*( x6(i-1) - x21(i-1) ) + G*C21_7*( x7(i-1) - x21(i-1) ) + G*C21_8*( x8(i-1) - x21(i-1) ) + G*C21_9*( x9(i-1) - x21(i-1) ) + G*C21_10*( x10(i-1) - x21(i-1) ) + G*C21_11*( x11(i-1) - x21(i-1) ) + G*C21_12*( x12(i-1) - x21(i-1) ) + G*C21_13*( x13(i-1) - x21(i-1) ) + G*C21_14*( x14(i-1) - x21(i-1) ) + G*C21_15*( x15(i-1) - x21(i-1) ) + G*C21_16*( x16(i-1) - x21(i-1) ) + G*C21_17*( x17(i-1) - x21(i-1) ) + G*C21_18*( x18(i-1) - x21(i-1) ) + G*C21_19*( x19(i-1) - x21(i-1) ) + G*C21_20*( x20(i-1) - x21(i-1) ) + G*C21_22*( x22(i-1) - x21(i-1) ) + G*C21_23*( x23(i-1) - x21(i-1) ) )*dt;
    y21(i) = y21(i-1) + ( ( a21 - x21(i-1)^2 - y21(i-1)^2 )*y21(i-1) + omega21*x21(i-1) + G*C21_1*( y1(i-1) - y21(i-1) ) + G*C21_2*( y2(i-1) - y21(i-1) ) + G*C21_3*( y3(i-1) - y21(i-1) ) + G*C21_4*( y4(i-1) - y21(i-1) ) + G*C21_5*( y5(i-1) - y21(i-1) ) + G*C21_6*( y6(i-1) - y21(i-1) ) + G*C21_7*( y7(i-1) - y21(i-1) ) + G*C21_8*( y8(i-1) - y21(i-1) ) + G*C21_9*( y9(i-1) - y21(i-1) ) + G*C21_10*( y10(i-1) - y21(i-1) ) + G*C21_11*( y11(i-1) - y21(i-1) ) + G*C21_12*( y12(i-1) - y21(i-1) ) + G*C21_13*( y13(i-1) - y21(i-1) ) + G*C21_14*( y14(i-1) - y21(i-1) ) + G*C21_15*( y15(i-1) - y21(i-1) ) + G*C21_16*( y16(i-1) - y21(i-1) ) + G*C21_17*( y17(i-1) - y21(i-1) ) + G*C21_18*( y18(i-1) - y21(i-1) ) + G*C21_19*( y19(i-1) - y21(i-1) ) + G*C21_20*( y20(i-1) - y21(i-1) ) + G*C21_22*( y22(i-1) - y21(i-1) ) + G*C21_23*( y23(i-1) - y21(i-1) ) )*dt;
    %Oscillator 22
    x22(i) = x22(i-1) + ( ( a22 - x22(i-1)^2 - y22(i-1)^2 )*x22(i-1) - omega22*y22(i-1) + G*C22_1*( x1(i-1) - x22(i-1) ) + G*C22_2*( x2(i-1) - x22(i-1) ) + G*C22_3*( x3(i-1) - x22(i-1) ) + G*C22_4*( x4(i-1) - x22(i-1) ) + G*C22_5*( x5(i-1) - x22(i-1) ) + G*C22_6*( x6(i-1) - x22(i-1) ) + G*C22_7*( x7(i-1) - x22(i-1) ) + G*C22_8*( x8(i-1) - x22(i-1) ) + G*C22_9*( x9(i-1) - x22(i-1) ) + G*C22_10*( x10(i-1) - x22(i-1) ) + G*C22_11*( x11(i-1) - x22(i-1) ) + G*C22_12*( x12(i-1) - x22(i-1) ) + G*C22_13*( x13(i-1) - x22(i-1) ) + G*C22_14*( x14(i-1) - x22(i-1) ) + G*C22_15*( x15(i-1) - x22(i-1) ) + G*C22_16*( x16(i-1) - x22(i-1) ) + G*C22_17*( x17(i-1) - x22(i-1) ) + G*C22_18*( x18(i-1) - x22(i-1) ) + G*C22_19*( x19(i-1) - x22(i-1) ) + G*C22_20*( x20(i-1) - x22(i-1) ) + G*C22_21*( x21(i-1) - x22(i-1) ) + G*C22_23*( x23(i-1) - x22(i-1) ) )*dt;
    y22(i) = y22(i-1) + ( ( a22 - x22(i-1)^2 - y22(i-1)^2 )*y22(i-1) + omega22*x22(i-1) + G*C22_1*( y1(i-1) - y22(i-1) ) + G*C22_2*( y2(i-1) - y22(i-1) ) + G*C22_3*( y3(i-1) - y22(i-1) ) + G*C22_4*( y4(i-1) - y22(i-1) ) + G*C22_5*( y5(i-1) - y22(i-1) ) + G*C22_6*( y6(i-1) - y22(i-1) ) + G*C22_7*( y7(i-1) - y22(i-1) ) + G*C22_8*( y8(i-1) - y22(i-1) ) + G*C22_9*( y9(i-1) - y22(i-1) ) + G*C22_10*( y10(i-1) - y22(i-1) ) + G*C22_11*( y11(i-1) - y22(i-1) ) + G*C22_12*( y12(i-1) - y22(i-1) ) + G*C22_13*( y13(i-1) - y22(i-1) ) + G*C22_14*( y14(i-1) - y22(i-1) ) + G*C22_15*( y15(i-1) - y22(i-1) ) + G*C22_16*( y16(i-1) - y22(i-1) ) + G*C22_17*( y17(i-1) - y22(i-1) ) + G*C22_18*( y18(i-1) - y22(i-1) ) + G*C22_19*( y19(i-1) - y22(i-1) ) + G*C22_20*( y20(i-1) - y22(i-1) ) + G*C22_21*( y21(i-1) - y22(i-1) ) + G*C22_23*( y23(i-1) - y22(i-1) ) )*dt;
    %Oscillator 23
    x23(i) = x23(i-1) + ( ( a23 - x23(i-1)^2 - y23(i-1)^2 )*x23(i-1) - omega23*y23(i-1) + G*C23_1*( x1(i-1) - x23(i-1) ) + G*C23_2*( x2(i-1) - x23(i-1) ) + G*C23_3*( x3(i-1) - x23(i-1) ) + G*C23_4*( x4(i-1) - x23(i-1) ) + G*C23_5*( x5(i-1) - x23(i-1) ) + G*C23_6*( x6(i-1) - x23(i-1) ) + G*C23_7*( x7(i-1) - x23(i-1) ) + G*C23_8*( x8(i-1) - x23(i-1) ) + G*C23_9*( x9(i-1) - x23(i-1) ) + G*C23_10*( x10(i-1) - x23(i-1) ) + G*C23_11*( x11(i-1) - x23(i-1) ) + G*C23_12*( x12(i-1) - x23(i-1) ) + G*C23_13*( x13(i-1) - x23(i-1) ) + G*C23_14*( x14(i-1) - x23(i-1) ) + G*C23_15*( x15(i-1) - x23(i-1) ) + G*C23_16*( x16(i-1) - x23(i-1) ) + G*C23_17*( x17(i-1) - x23(i-1) ) + G*C23_18*( x18(i-1) - x23(i-1) ) + G*C23_19*( x19(i-1) - x23(i-1) ) + G*C23_20*( x20(i-1) - x23(i-1) ) + G*C23_21*( x21(i-1) - x23(i-1) ) + G*C23_22*( x22(i-1) - x23(i-1) ) )*dt;
    y23(i) = y23(i-1) + ( ( a23 - x23(i-1)^2 - y23(i-1)^2 )*y23(i-1) + omega23*x23(i-1) + G*C23_1*( y1(i-1) - y23(i-1) ) + G*C23_2*( y2(i-1) - y23(i-1) ) + G*C23_3*( y3(i-1) - y23(i-1) ) + G*C23_4*( y4(i-1) - y23(i-1) ) + G*C23_5*( y5(i-1) - y23(i-1) ) + G*C23_6*( y6(i-1) - y23(i-1) ) + G*C23_7*( y7(i-1) - y23(i-1) ) + G*C23_8*( y8(i-1) - y23(i-1) ) + G*C23_9*( y9(i-1) - y23(i-1) ) + G*C23_10*( y10(i-1) - y23(i-1) ) + G*C23_11*( y11(i-1) - y23(i-1) ) + G*C23_12*( y12(i-1) - y23(i-1) ) + G*C23_13*( y13(i-1) - y23(i-1) ) + G*C23_14*( y14(i-1) - y23(i-1) ) + G*C23_15*( y15(i-1) - y23(i-1) ) + G*C23_16*( y16(i-1) - y23(i-1) ) + G*C23_17*( y17(i-1) - y23(i-1) ) + G*C23_18*( y18(i-1) - y23(i-1) ) + G*C23_19*( y19(i-1) - y23(i-1) ) + G*C23_20*( y20(i-1) - y23(i-1) ) + G*C23_21*( y21(i-1) - y23(i-1) ) + G*C23_22*( y22(i-1) - y23(i-1) ) )*dt;
end
plot (x1)