Help in the value of a variable.

2 Ansichten (letzte 30 Tage)
Haya Ali
Haya Ali am 15 Feb. 2023
Kommentiert: Haya Ali am 15 Feb. 2023
I am advised to use the for loop for std and measure the value of P1 for each value of std. The value of P1 should be like this
P1 =
Columns 1 through 13
0 0.1429 0.2857 0.4286 0.5714 0.7143 0.8571 1.0000 1.1429 1.2857 1.4286 1.5714 1.7143
Columns 14 through 21
1.8571 2.0000 2.1429 2.2857 2.4286 2.5714 2.7143 2.8571
but when i run my code I get only a single value for P1. Please help me to resolve this. Below is my code.
%
clear all; close all; clc;
% value of constants%
G = 10;
dt = 0.01; % step size
beta=0.0001;
%Oscillator 1%
a1=0.7; omega1=20; C1_2 =0; C1_3 =0; C1_4 =0.00101; C1_5 =0; C1_6 =0; C1_7 =0.00102; C1_8 =0; C1_9 =0; C1_10 =0; C1_11 =0; C1_12 =0.00103; C1_13 =0; C1_14 =0; C1_15 =0; C1_16 =0; C1_17 =0.00104; C1_18 =0; C1_19 =0; C1_20 =0; C1_21 =0; C1_22 =0.00105; C1_23 =0.00106;
%Oscillator 2%
a2 = 0.1; omega2 = 40; C2_1 =0.003; C2_3 =0; C2_4 =0.0031; C2_5 =0; C2_6 =0; C2_7 =0; C2_8 =0.0032; C2_9 =0; C2_10 =0; C2_11 =0; C2_12 =0.0033; C2_13 =0; C2_14 =0; C2_15 =0; C2_16 =0.0034; C2_17 =0; C2_18 =0; C2_19 =0; C2_20 =0.0035; C2_21 =0; C2_22 =0; C2_23 =0;
%Oscillator 3%
a3 = 0.8; omega3 = 30; C3_1 =0; C3_2 =0; C3_4 =0.0086; C3_5 =0.008; C3_6 =0; C3_7 =0; C3_8 =0.0081; C3_9 =0; C3_10 =0; C3_11 =0; C3_12 =0; C3_13 =0; C3_14 =0.0083; C3_15 =0; C3_16 =0; C3_17 =0.0084; C3_18 =0; C3_19 =0; C3_20 =0; C3_21 =0; C3_22 =0; C3_23 =0.0085;
%Oscillator 4%
a4 = 0.5; omega4 = 60; C4_1 =0.005; C4_2 =0; C4_3 =0.0055; C4_5 =0; C4_6 =0; C4_7 =0; C4_8 =0; C4_9 =0.0051; C4_10 =0; C4_11 =0; C4_12 =0; C4_13 =0.0052; C4_14 =0; C4_15 =0; C4_16 =0; C4_17 =0.0053; C4_18 =0; C4_19 =0; C4_20 =0; C4_21 =0.0054; C4_22 =0; C4_23 =0;
%Oscillator 5%
a5 = 0.6; omega5 = 80; C5_1 =0.006; C5_2 =0; C5_3 =0.0065; C5_4 =0; C5_6 =0; C5_7 =0; C5_8 =0; C5_9 =0; C5_10 =0; C5_11 =0.0062; C5_12 =0; C5_13 =0; C5_14 =0; C5_15 =0.0063; C5_16 =0; C5_17 =0; C5_18 =0; C5_19 =0.0064; C5_20 =0; C5_21 =0; C5_22 =0; C5_23 =0;
%Oscillator 6%
a6 = 0.1; omega6 = 50; C6_1 =0.01; C6_2 =0; C6_3 =0.015; C6_4 =0; C6_5 =0; C6_7 =0; C6_8 =0; C6_9 =0.012; C6_10 =0; C6_11 =0; C6_12 =0; C6_13 =0; C6_14 =0; C6_15 =0; C6_16 =0; C6_17 =0.013; C6_18 =0; C6_19 =0; C6_20 =0; C6_21 =0; C6_22 =0; C6_23 =0.014;
%Oscillator 7%
a7 = 0.5; omega7 = 20; C7_1 =0; C7_2 =0; C7_3 =0.0045; C7_4 =0.004; C7_5 =0; C7_6 =0; C7_8 =0; C7_9 =0.0042; C7_10 =0; C7_11 =0; C7_12 =0; C7_13 =0; C7_14 =0; C7_15 =0.0043; C7_16 =0; C7_17 =0; C7_18 =0; C7_19 =0.0044; C7_20 =0; C7_21 =0; C7_22 =0.0041; C7_23 =0;
%Oscillator 8%
a8 = 0.2; omega8 = 70; C8_1 =0.002; C8_2 =0; C8_3 =0.0021; C8_4 =0; C8_5 =0; C8_6 =0; C8_7 =0; C8_9 =0; C8_10 =0.0021; C8_11 =0; C8_12 =0; C8_13 =0; C8_14 =0; C8_15 =0; C8_16 =0.0022; C8_17 =0; C8_18 =0; C8_19 =0; C8_20 =0.0023; C8_21 =0; C8_22 =0; C8_23 =0.0024;
%Oscillator 9%
a9 = 0.7; omega9 = 60; C9_1 =0; C9_2 =0; C9_3 =0.0075; C9_4 =0; C9_5 =0.007; C9_6 =0; C9_7 =0; C9_8 =0; C9_10 =0.0071; C9_11 =0; C9_12 =0; C9_13 =0; C9_14 =0; C9_15 =0; C9_16 =0.0072; C9_17 =0; C9_18 =0.0073; C9_19 =0; C9_20 =0; C9_21 =0.0074; C9_22 =0; C9_23 =0;
%Oscillator 10%
a10 = 0.3; omega10 = 50; C10_1 =0; C10_2 =0.009; C10_3 =0.0095; C10_4 =0; C10_5 =0; C10_6 =0; C10_7 =0; C10_8 =0.0091; C10_9 =0; C10_11 =0; C10_12 =0; C10_13 =0; C10_14 =0.0092; C10_15 =0; C10_16 =0; C10_17 =0; C10_18 =0.0093; C10_19 =0.0094; C10_20 =0; C10_21 =0; C10_22 =0; C10_23 =0;
%Oscillator 11%
a11 = 0.20; omega11 = 20; C11_1 =0.020; C11_2 =0; C11_3 =0.0205; C11_4 =0; C11_5 =0; C11_6 =0; C11_7 =0.0201; C11_8 =0; C11_9 =0; C11_10 =0; C11_12 =0; C11_13 =0; C11_14 =0; C11_15 =0; C11_16 =0; C11_17 =0.0202; C11_18 =0; C11_19 =0.0203; C11_20 =0; C11_21 =0; C11_22 =0; C11_23 =0.0204;
%Oscillator 12%
a12 = 0.13; omega12 = 90; C12_1 =0.013; C12_2 =0; C12_3 =0.01315; C12_4 =0; C12_5 =0; C12_6 =0; C12_7 =0.0131; C12_8 =0; C12_9 =0; C12_10 =0; C12_11 =0; C12_13 =0; C12_14 =0; C12_15 =0.0132; C12_16 =0; C12_17 =0; C12_18 =0.0133; C12_19 =0; C12_20 =0; C12_21 =0; C12_22 =0; C12_23 =0.0134;
%Oscillator 13%
a13 = 0.18; omega13 = 40; C13_1 =0; C13_2 =0; C13_3 =0.0185; C13_4 =0; C13_5 =0; C13_6 =0.018; C13_7 =0; C13_8 =0; C13_9 =0.0181; C13_10 =0; C13_11 =0; C13_12 =0; C13_14 =0; C13_15 =0.0182; C13_16 =0; C13_17 =0; C13_18 =0; C13_19 =0.0183; C13_20 =0; C13_21 =0; C13_22 =0.0184; C13_23 =0;
%Oscillator 14%
a14 = 0.23; omega14 = 30; C14_1 =0.023; C14_2 =0; C14_3 =0.0235; C14_4 =0; C14_5 =0; C14_6 =0; C14_7 =0.0231; C14_8 =0; C14_9 =0; C14_10 =0.0232; C14_11 =0; C14_12 =0; C14_13 =0; C14_15 =0; C14_16 =0.0233; C14_17 =0; C14_18 =0; C14_19 =0; C14_20 =0.0234; C14_21 =0; C14_22 =0; C14_23 =0;
%Oscillator 15%
a15 = 0.22; omega15 = 60; C15_1 =0; C15_2 =0.022; C15_3 =0.0225; C15_4 =0; C15_5 =0; C15_6 =0; C15_7 =0; C15_8 =0.0221; C15_9 =0; C15_10 =0; C15_11 =0; C15_12 =0; C15_13 =0.0222; C15_14 =0; C15_16 =0; C15_17 =0; C15_18 =0.0223; C15_19 =0; C15_20 =0.0224; C15_21 =0; C15_22 =0; C15_23 =0;
%Oscillator 16%
a16 = 0.11; omega16 = 50; C16_1 =0; C16_2 =0.011; C16_3 =0.01115; C16_4 =0; C16_5 =0; C16_6 =0; C16_7 =0; C16_8 =0; C16_9 =0.0111; C16_10 =0; C16_11 =0; C16_12 =0; C16_13 =0; C16_14 =0.0112; C16_15 =0; C16_17 =0; C16_18 =0; C16_19 =0.0113; C16_20 =0; C16_21 =0.01114; C16_22 =0; C16_23 =0;
%Oscillator 17%
a17 = 0.14; omega17 = 80; C17_1 =0.014; C17_2 =0; C17_3 =0.0145; C17_4 =0; C17_5 =0; C17_6 =0; C17_7 =0; C17_8 =0; C17_9=0.0141; C17_10 =0; C17_11 =0; C17_12 =0; C17_13 =0.0142; C17_14 =0; C17_15 =0; C17_16 =0; C17_18 =0; C17_19 =0.0143; C17_20 =0; C17_21 =0; C17_22 =0.0144; C17_23 =0;
%Oscillator 18%
a18 = 0.21; omega18 = 70; C18_1 =0.021; C18_2 =0.0211; C18_3 =0.0215; C18_4 =0; C18_5 =0; C18_6 =0; C18_7 =0; C18_8 =0; C18_9 =0; C18_10 =0.0212; C18_11 =0; C18_12 =0; C18_13 =0; C18_14 =0; C18_15 =0; C18_16 =0.0213; C18_17 =0; C18_19 =0; C18_20 =0.0214; C18_21 =0; C18_22 =0; C18_23 =0;
%Oscillator 19%
a19 = 0.17; omega19 = 30; C19_1 =0; C19_2 =0.017; C19_3 =0.0175; C19_4 =0; C19_5 =0; C19_6 =0; C19_7 =0; C19_8 =0.0171; C19_9 =0; C19_10 =0; C19_11 =0.0172; C19_12 =0; C19_13 =0; C19_14 =0; C19_15 =0; C19_16 =0; C19_17 =0.0173; C19_18 =0; C19_20 =0; C19_21 =0.0174; C19_22 =0; C19_23 =0;
%Oscillator 20%
a20 = 0.19; omega20 = 90; C20_1 =0; C20_2 =0; C20_3 =0.01915; C20_4 =0; C20_5 =0; C20_6 =0.0191; C20_7 =0; C20_8 =0; C20_9 =0; C20_10 =0; C20_11 =0; C20_12 =0.0192; C20_13 =0; C20_14 =0; C20_15 =0; C20_16 =0.0193; C20_17 =0; C20_18 =0.0194; C20_19 =0; C20_21 =0; C20_22 =0; C20_23 =0;
%Oscillator 21%
a21 = 0.15; omega21 = 20; C21_1 =0.015; C21_2 =0.0151; C21_3 =0.0155; C21_4 =0; C21_5 =0; C21_6 =0; C21_7 =0; C21_8 =0; C21_9 =0; C21_10 =0; C21_11 =0.0152; C21_12 =0; C21_13 =0; C21_14 =0; C21_15 =0; C21_16 =0.0153; C21_17 =0; C21_18 =0; C21_19 =0.0154; C21_20 =0; C21_22 =0; C21_23 =0;
%Oscillator 22%
a22 = 0.12; omega22 = 30; C22_1 =0; C22_2 =0; C22_3 =0.01215; C22_4 =0; C22_5 =0.0121; C22_6 =0; C22_7 =0; C22_8 =0; C22_9 =0.0122; C22_10 =0; C22_11 =0; C22_12 =0; C22_13 =0; C22_14 =0.0123; C22_15 =0; C22_16 =0; C22_17 =0; C22_18 =0; C22_19 =0; C22_20 =0; C22_21 =0; C22_23 =0.0124;
%Oscillator 23%
a23 = 0.16; omega23 = 40; C23_1 =0.016; C23_2 =0; C23_3 =0.0165; C23_4 =0; C23_5 =0; C23_6 =0.0161; C23_7 =0; C23_8 =0; C23_9 =0; C23_10 =0.0162; C23_11 =0; C23_12 =0; C23_13 =0; C23_14 =0; C23_15 =0.0163; C23_16 =0; C23_17 =0; C23_18 =0; C23_19 =0; C23_20 =0; C23_21 =0.0164; C23_22 =0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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 std=0:0.1:1
eta=normrnd(0,std,[1 10]); % Gaussian Noise
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) ) +beta*normrnd(0,std,[1 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) ) +beta*normrnd(0,std,[1 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
end
P1= std/a1;

Akzeptierte Antwort

Shushant
Shushant am 15 Feb. 2023
Here is the modified code, the only mistake you were doing is calculating the value of P1 after the loop has ended. The value of P1 should be calculated inside the loop and P1 should be a vector. hence I have declared a variable called idx which is initailized at 1 and increments its value after every iteration.
%
clear all; close all; clc;
% value of constants%
G = 10;
dt = 0.01; % step size
beta=0.0001;
idx = 1; % index for the P1
%Oscillator 1%
a1=0.7; omega1=20; C1_2 =0; C1_3 =0; C1_4 =0.00101; C1_5 =0; C1_6 =0; C1_7 =0.00102; C1_8 =0; C1_9 =0; C1_10 =0; C1_11 =0; C1_12 =0.00103; C1_13 =0; C1_14 =0; C1_15 =0; C1_16 =0; C1_17 =0.00104; C1_18 =0; C1_19 =0; C1_20 =0; C1_21 =0; C1_22 =0.00105; C1_23 =0.00106;
%Oscillator 2%
a2 = 0.1; omega2 = 40; C2_1 =0.003; C2_3 =0; C2_4 =0.0031; C2_5 =0; C2_6 =0; C2_7 =0; C2_8 =0.0032; C2_9 =0; C2_10 =0; C2_11 =0; C2_12 =0.0033; C2_13 =0; C2_14 =0; C2_15 =0; C2_16 =0.0034; C2_17 =0; C2_18 =0; C2_19 =0; C2_20 =0.0035; C2_21 =0; C2_22 =0; C2_23 =0;
%Oscillator 3%
a3 = 0.8; omega3 = 30; C3_1 =0; C3_2 =0; C3_4 =0.0086; C3_5 =0.008; C3_6 =0; C3_7 =0; C3_8 =0.0081; C3_9 =0; C3_10 =0; C3_11 =0; C3_12 =0; C3_13 =0; C3_14 =0.0083; C3_15 =0; C3_16 =0; C3_17 =0.0084; C3_18 =0; C3_19 =0; C3_20 =0; C3_21 =0; C3_22 =0; C3_23 =0.0085;
%Oscillator 4%
a4 = 0.5; omega4 = 60; C4_1 =0.005; C4_2 =0; C4_3 =0.0055; C4_5 =0; C4_6 =0; C4_7 =0; C4_8 =0; C4_9 =0.0051; C4_10 =0; C4_11 =0; C4_12 =0; C4_13 =0.0052; C4_14 =0; C4_15 =0; C4_16 =0; C4_17 =0.0053; C4_18 =0; C4_19 =0; C4_20 =0; C4_21 =0.0054; C4_22 =0; C4_23 =0;
%Oscillator 5%
a5 = 0.6; omega5 = 80; C5_1 =0.006; C5_2 =0; C5_3 =0.0065; C5_4 =0; C5_6 =0; C5_7 =0; C5_8 =0; C5_9 =0; C5_10 =0; C5_11 =0.0062; C5_12 =0; C5_13 =0; C5_14 =0; C5_15 =0.0063; C5_16 =0; C5_17 =0; C5_18 =0; C5_19 =0.0064; C5_20 =0; C5_21 =0; C5_22 =0; C5_23 =0;
%Oscillator 6%
a6 = 0.1; omega6 = 50; C6_1 =0.01; C6_2 =0; C6_3 =0.015; C6_4 =0; C6_5 =0; C6_7 =0; C6_8 =0; C6_9 =0.012; C6_10 =0; C6_11 =0; C6_12 =0; C6_13 =0; C6_14 =0; C6_15 =0; C6_16 =0; C6_17 =0.013; C6_18 =0; C6_19 =0; C6_20 =0; C6_21 =0; C6_22 =0; C6_23 =0.014;
%Oscillator 7%
a7 = 0.5; omega7 = 20; C7_1 =0; C7_2 =0; C7_3 =0.0045; C7_4 =0.004; C7_5 =0; C7_6 =0; C7_8 =0; C7_9 =0.0042; C7_10 =0; C7_11 =0; C7_12 =0; C7_13 =0; C7_14 =0; C7_15 =0.0043; C7_16 =0; C7_17 =0; C7_18 =0; C7_19 =0.0044; C7_20 =0; C7_21 =0; C7_22 =0.0041; C7_23 =0;
%Oscillator 8%
a8 = 0.2; omega8 = 70; C8_1 =0.002; C8_2 =0; C8_3 =0.0021; C8_4 =0; C8_5 =0; C8_6 =0; C8_7 =0; C8_9 =0; C8_10 =0.0021; C8_11 =0; C8_12 =0; C8_13 =0; C8_14 =0; C8_15 =0; C8_16 =0.0022; C8_17 =0; C8_18 =0; C8_19 =0; C8_20 =0.0023; C8_21 =0; C8_22 =0; C8_23 =0.0024;
%Oscillator 9%
a9 = 0.7; omega9 = 60; C9_1 =0; C9_2 =0; C9_3 =0.0075; C9_4 =0; C9_5 =0.007; C9_6 =0; C9_7 =0; C9_8 =0; C9_10 =0.0071; C9_11 =0; C9_12 =0; C9_13 =0; C9_14 =0; C9_15 =0; C9_16 =0.0072; C9_17 =0; C9_18 =0.0073; C9_19 =0; C9_20 =0; C9_21 =0.0074; C9_22 =0; C9_23 =0;
%Oscillator 10%
a10 = 0.3; omega10 = 50; C10_1 =0; C10_2 =0.009; C10_3 =0.0095; C10_4 =0; C10_5 =0; C10_6 =0; C10_7 =0; C10_8 =0.0091; C10_9 =0; C10_11 =0; C10_12 =0; C10_13 =0; C10_14 =0.0092; C10_15 =0; C10_16 =0; C10_17 =0; C10_18 =0.0093; C10_19 =0.0094; C10_20 =0; C10_21 =0; C10_22 =0; C10_23 =0;
%Oscillator 11%
a11 = 0.20; omega11 = 20; C11_1 =0.020; C11_2 =0; C11_3 =0.0205; C11_4 =0; C11_5 =0; C11_6 =0; C11_7 =0.0201; C11_8 =0; C11_9 =0; C11_10 =0; C11_12 =0; C11_13 =0; C11_14 =0; C11_15 =0; C11_16 =0; C11_17 =0.0202; C11_18 =0; C11_19 =0.0203; C11_20 =0; C11_21 =0; C11_22 =0; C11_23 =0.0204;
%Oscillator 12%
a12 = 0.13; omega12 = 90; C12_1 =0.013; C12_2 =0; C12_3 =0.01315; C12_4 =0; C12_5 =0; C12_6 =0; C12_7 =0.0131; C12_8 =0; C12_9 =0; C12_10 =0; C12_11 =0; C12_13 =0; C12_14 =0; C12_15 =0.0132; C12_16 =0; C12_17 =0; C12_18 =0.0133; C12_19 =0; C12_20 =0; C12_21 =0; C12_22 =0; C12_23 =0.0134;
%Oscillator 13%
a13 = 0.18; omega13 = 40; C13_1 =0; C13_2 =0; C13_3 =0.0185; C13_4 =0; C13_5 =0; C13_6 =0.018; C13_7 =0; C13_8 =0; C13_9 =0.0181; C13_10 =0; C13_11 =0; C13_12 =0; C13_14 =0; C13_15 =0.0182; C13_16 =0; C13_17 =0; C13_18 =0; C13_19 =0.0183; C13_20 =0; C13_21 =0; C13_22 =0.0184; C13_23 =0;
%Oscillator 14%
a14 = 0.23; omega14 = 30; C14_1 =0.023; C14_2 =0; C14_3 =0.0235; C14_4 =0; C14_5 =0; C14_6 =0; C14_7 =0.0231; C14_8 =0; C14_9 =0; C14_10 =0.0232; C14_11 =0; C14_12 =0; C14_13 =0; C14_15 =0; C14_16 =0.0233; C14_17 =0; C14_18 =0; C14_19 =0; C14_20 =0.0234; C14_21 =0; C14_22 =0; C14_23 =0;
%Oscillator 15%
a15 = 0.22; omega15 = 60; C15_1 =0; C15_2 =0.022; C15_3 =0.0225; C15_4 =0; C15_5 =0; C15_6 =0; C15_7 =0; C15_8 =0.0221; C15_9 =0; C15_10 =0; C15_11 =0; C15_12 =0; C15_13 =0.0222; C15_14 =0; C15_16 =0; C15_17 =0; C15_18 =0.0223; C15_19 =0; C15_20 =0.0224; C15_21 =0; C15_22 =0; C15_23 =0;
%Oscillator 16%
a16 = 0.11; omega16 = 50; C16_1 =0; C16_2 =0.011; C16_3 =0.01115; C16_4 =0; C16_5 =0; C16_6 =0; C16_7 =0; C16_8 =0; C16_9 =0.0111; C16_10 =0; C16_11 =0; C16_12 =0; C16_13 =0; C16_14 =0.0112; C16_15 =0; C16_17 =0; C16_18 =0; C16_19 =0.0113; C16_20 =0; C16_21 =0.01114; C16_22 =0; C16_23 =0;
%Oscillator 17%
a17 = 0.14; omega17 = 80; C17_1 =0.014; C17_2 =0; C17_3 =0.0145; C17_4 =0; C17_5 =0; C17_6 =0; C17_7 =0; C17_8 =0; C17_9=0.0141; C17_10 =0; C17_11 =0; C17_12 =0; C17_13 =0.0142; C17_14 =0; C17_15 =0; C17_16 =0; C17_18 =0; C17_19 =0.0143; C17_20 =0; C17_21 =0; C17_22 =0.0144; C17_23 =0;
%Oscillator 18%
a18 = 0.21; omega18 = 70; C18_1 =0.021; C18_2 =0.0211; C18_3 =0.0215; C18_4 =0; C18_5 =0; C18_6 =0; C18_7 =0; C18_8 =0; C18_9 =0; C18_10 =0.0212; C18_11 =0; C18_12 =0; C18_13 =0; C18_14 =0; C18_15 =0; C18_16 =0.0213; C18_17 =0; C18_19 =0; C18_20 =0.0214; C18_21 =0; C18_22 =0; C18_23 =0;
%Oscillator 19%
a19 = 0.17; omega19 = 30; C19_1 =0; C19_2 =0.017; C19_3 =0.0175; C19_4 =0; C19_5 =0; C19_6 =0; C19_7 =0; C19_8 =0.0171; C19_9 =0; C19_10 =0; C19_11 =0.0172; C19_12 =0; C19_13 =0; C19_14 =0; C19_15 =0; C19_16 =0; C19_17 =0.0173; C19_18 =0; C19_20 =0; C19_21 =0.0174; C19_22 =0; C19_23 =0;
%Oscillator 20%
a20 = 0.19; omega20 = 90; C20_1 =0; C20_2 =0; C20_3 =0.01915; C20_4 =0; C20_5 =0; C20_6 =0.0191; C20_7 =0; C20_8 =0; C20_9 =0; C20_10 =0; C20_11 =0; C20_12 =0.0192; C20_13 =0; C20_14 =0; C20_15 =0; C20_16 =0.0193; C20_17 =0; C20_18 =0.0194; C20_19 =0; C20_21 =0; C20_22 =0; C20_23 =0;
%Oscillator 21%
a21 = 0.15; omega21 = 20; C21_1 =0.015; C21_2 =0.0151; C21_3 =0.0155; C21_4 =0; C21_5 =0; C21_6 =0; C21_7 =0; C21_8 =0; C21_9 =0; C21_10 =0; C21_11 =0.0152; C21_12 =0; C21_13 =0; C21_14 =0; C21_15 =0; C21_16 =0.0153; C21_17 =0; C21_18 =0; C21_19 =0.0154; C21_20 =0; C21_22 =0; C21_23 =0;
%Oscillator 22%
a22 = 0.12; omega22 = 30; C22_1 =0; C22_2 =0; C22_3 =0.01215; C22_4 =0; C22_5 =0.0121; C22_6 =0; C22_7 =0; C22_8 =0; C22_9 =0.0122; C22_10 =0; C22_11 =0; C22_12 =0; C22_13 =0; C22_14 =0.0123; C22_15 =0; C22_16 =0; C22_17 =0; C22_18 =0; C22_19 =0; C22_20 =0; C22_21 =0; C22_23 =0.0124;
%Oscillator 23%
a23 = 0.16; omega23 = 40; C23_1 =0.016; C23_2 =0; C23_3 =0.0165; C23_4 =0; C23_5 =0; C23_6 =0.0161; C23_7 =0; C23_8 =0; C23_9 =0; C23_10 =0.0162; C23_11 =0; C23_12 =0; C23_13 =0; C23_14 =0; C23_15 =0.0163; C23_16 =0; C23_17 =0; C23_18 =0; C23_19 =0; C23_20 =0; C23_21 =0.0164; C23_22 =0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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 std=0:0.1:2 % MODIFIED to 2 to get 21 columns instead of 11 columns
eta=normrnd(0,std,[1 10]); % Gaussian Noise
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) ) +beta*normrnd(0,std,[1 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) ) +beta*normrnd(0,std,[1 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
P1(idx) = std/a1; %calculating P1 after every iteration
idx = idx +1; %increasing the value of idx
end
P1
P1 = 1×21
0 0.1429 0.2857 0.4286 0.5714 0.7143 0.8571 1.0000 1.1429 1.2857 1.4286 1.5714 1.7143 1.8571 2.0000 2.1429 2.2857 2.4286 2.5714 2.7143 2.8571

Weitere Antworten (0)

Community Treasure Hunt

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

Start Hunting!

Translated by