clear ;
global TayOS; TayOS = 20 ;
format short e;
format compact;
global index;
global CN1; global CN2; global CN3; global CN4; global CN5; global CN6;
global q1; global q2; global q3; global q4; global q5; global q6;
global Nag_Ct1; global Nag_Ct2; global Nag_Ct3; global Nag_Ct4; global Nag_Ct5; global Nag_Ct6;
CN1=50; CN2=50; CN3=50; CN4=50; CN5=50; CN6=50;
q1=3000; q2=3000; q3=3000; q4=3000; q5=3000; q6=3000;
Nag_Ct1=60; Nag_Ct2=60; Nag_Ct3=60; Nag_Ct4=60; Nag_Ct5=60; Nag_Ct6=60;
index=0;
TnSteni = 20 ;
TnTqlo = 20 ;
TnNagr = 20 ;
TnT2 = 20 ;
StSteni = [ TnSteni;TnSteni;TnSteni;TnSteni; TnSteni;TnSteni;TnSteni;TnSteni; TnSteni;TnSteni;TnSteni;TnSteni; TnSteni;TnSteni;TnSteni;TnSteni; TnSteni;TnSteni;TnSteni;TnSteni; TnSteni;TnSteni;TnSteni;TnSteni; ];
StTqlo = [ TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo;TnTqlo; TnTqlo;];
StNagr = [ TnNagr;TnNagr;TnNagr;TnNagr;TnNagr;TnNagr; ];
StT2 = [ TnT2;TnT2;TnT2;TnT2;TnT2;TnT2; ];
[t,y]=ode15s (@ExRK3,[0 2000000],[StSteni; StTqlo; StNagr; StT2;]);
long_m = length(t)
YS1_1 = y(1:long_m, 1);
YS1_2 = y(1:long_m, 2);
YS1_3 = y(1:long_m, 3);
YS1_4 = y(1:long_m, 4);
YS2_1 = y(1:long_m, 5);
YS2_2 = y(1:long_m, 6);
YS2_3 = y(1:long_m, 7);
YS2_4 = y(1:long_m, 8);
YS3_1 = y(1:long_m, 9);
YS3_2 = y(1:long_m, 10);
YS3_3 = y(1:long_m, 11);
YS3_4 = y(1:long_m, 12);
YS4_1 = y(1:long_m, 13);
YS4_2 = y(1:long_m, 14);
YS4_3 = y(1:long_m, 15);
YS4_4 = y(1:long_m, 16);
YS5_1 = y(1:long_m, 17);
YS5_2 = y(1:long_m, 18);
YS5_3 = y(1:long_m, 19);
YS5_4 = y(1:long_m, 20);
YS6_1 = y(1:long_m, 21);
YS6_2 = y(1:long_m, 22);
YS6_3 = y(1:long_m, 23);
YS6_4 = y(1:long_m, 24);
YT1_1 = y(1:long_m, 25);
YT1_2 = y(1:long_m, 26);
YT1_3 = y(1:long_m, 27);
YT1_4 = y(1:long_m, 28);
YT1_5 = y(1:long_m, 29);
YT1_6 = y(1:long_m, 30);
YT1_7 = y(1:long_m, 31);
YT1_8 = y(1:long_m, 32);
YT1_9 = y(1:long_m, 33);
YT1_10 = y(1:long_m, 34);
YT2_1 = y(1:long_m, 35);
YT2_2 = y(1:long_m, 36);
YT2_3 = y(1:long_m, 37);
YT2_4 = y(1:long_m, 38);
YT2_5 = y(1:long_m, 39);
YT2_6 = y(1:long_m, 40);
YT2_7 = y(1:long_m, 41);
YT2_8 = y(1:long_m, 42);
YT2_9 = y(1:long_m, 43);
YT2_10 = y(1:long_m, 44);
YT3_1 = y(1:long_m, 45);
YT3_2 = y(1:long_m, 46);
YT3_3 = y(1:long_m, 47);
YT3_4 = y(1:long_m, 48);
YT3_5 = y(1:long_m, 49);
YT3_6 = y(1:long_m, 50);
YT3_7 = y(1:long_m, 51);
YT3_8 = y(1:long_m, 52);
YT3_9 = y(1:long_m, 53);
YT3_10 = y(1:long_m, 54);
YT4_1 = y(1:long_m, 55);
YT4_2 = y(1:long_m, 56);
YT4_3 = y(1:long_m, 57);
YT4_4 = y(1:long_m, 58);
YT4_5 = y(1:long_m, 59);
YT4_6 = y(1:long_m, 60);
YT4_7 = y(1:long_m, 61);
YT4_8 = y(1:long_m, 62);
YT4_9 = y(1:long_m, 63);
YT4_10 = y(1:long_m, 64);
YT5_1 = y(1:long_m, 65);
YT5_2 = y(1:long_m, 66);
YT5_3 = y(1:long_m, 67);
YT5_4 = y(1:long_m, 68);
YT5_5 = y(1:long_m, 69);
YT5_6 = y(1:long_m, 70);
YT5_7 = y(1:long_m, 71);
YT5_8 = y(1:long_m, 72);
YT5_9 = y(1:long_m, 73);
YT5_10 = y(1:long_m, 74);
YT6_1 = y(1:long_m, 75);
YT6_2 = y(1:long_m, 76);
YT6_3 = y(1:long_m, 77);
YT6_4 = y(1:long_m, 78);
YT6_5 = y(1:long_m, 79);
YT6_6 = y(1:long_m, 80);
YT6_7 = y(1:long_m, 81);
YT6_8 = y(1:long_m, 82);
YT6_9 = y(1:long_m, 83);
YT6_10 = y(1:long_m, 84);
YTC = y(1:long_m, 85);
QT1 = y(1:long_m, 86);
QT2 = y(1:long_m, 87);
QT3 = y(1:long_m, 88);
QT4 = y(1:long_m, 89);
QT5 = y(1:long_m, 90);
QT6 = y(1:long_m, 91);
Tr2_1 = y(1:long_m, 92);
Tr2_2 = y(1:long_m, 93);
Tr2_3 = y(1:long_m, 94);
Tr2_4 = y(1:long_m, 95);
Tr2_5 = y(1:long_m, 96);
Tr2_6 = y(1:long_m, 97);
save my_file_ode
plot(t, y)
grid
pause
plot( t,YS1_1,'-r', t,YS1_2,'-r', t,YS1_3,'-r', t,YS1_4,'-r', t,QT1,'--r', t,Tr2_1,'--b')
grid ;
pause ;
plot( t,YS2_1,'-r', t,YS2_2,'-r', t,YS2_3,'-r', t,YS2_4,'-r', t,QT2,'--r', t,Tr2_2,'--b')
grid ;
pause ;
plot( t,YS3_1,'-r', t,YS3_2,'-r', t,YS3_3,'-r', t,YS3_4,'-r', t,QT3,'--r', t,Tr2_3,'--b')
grid ;
pause ;
plot( t,YS4_1,'-r', t,YS4_2,'-r', t,YS4_3,'-r', t,YS4_4,'-r', t,QT4,'--r', t,Tr2_4,'--b')
grid ;
pause ;
plot( t,YS5_1,'-r', t,YS5_2,'-r', t,YS5_3,'-r', t,YS5_4,'-r', t,QT5,'--r', t,Tr2_5,'--b')
grid ;
pause ;
plot( t,YS6_1,'-r', t,YS6_2,'-r', t,YS6_3,'-r', t,YS6_4,'-r', t,QT6,'--r', t,Tr2_6,'--b')
grid ;
pause ;
plot( t,YT1_1,'-b', t,YT1_2,'-b', t,YT1_3,'-b', t,YT1_4,'-b', t,YT1_5,'-b', t,YT1_6,'-b', t,YT1_7,'-b', t,YT1_8,'-b', t,YT1_9,'-b', t,YT1_10,'-b', t,YTC,'--g', t,QT1,'--r', t,Tr2_1,'--b')
grid ;
pause ;
plot( t,YT2_1,'-b', t,YT2_2,'-b', t,YT2_3,'-b', t,YT2_4,'-b', t,YT2_5,'-b', t,YT2_6,'-b', t,YT2_7,'-b', t,YT2_8,'-b', t,YT2_9,'-b', t,YT2_10,'-b', t,YTC,'--g', t,QT2,'--r', t,Tr2_2,'--b')
grid ;
pause ;
plot( t,YT3_1,'-b', t,YT3_2,'-b', t,YT3_3,'-b', t,YT3_4,'-b', t,YT3_5,'-b', t,YT3_6,'-b', t,YT3_7,'-b', t,YT3_8,'-b', t,YT3_9,'-b', t,YT3_10,'-b', t,YTC,'--g', t,QT3,'--r', t,Tr2_3,'--b')
grid ;
pause ;
plot( t,YT4_1,'-b', t,YT4_2,'-b', t,YT4_3,'-b', t,YT4_4,'-b', t,YT4_5,'-b', t,YT4_6,'-b', t,YT4_7,'-b', t,YT4_8,'-b', t,YT4_9,'-b', t,YT4_10,'-b', t,YTC,'--g', t,QT4,'--r', t,Tr2_4,'--b')
grid ;
pause ;
plot( t,YT5_1,'-b', t,YT5_2,'-b', t,YT5_3,'-b', t,YT5_4,'-b', t,YT5_5,'-b', t,YT5_6,'-b', t,YT5_7,'-b', t,YT5_8,'-b', t,YT5_9,'-b', t,YT5_10,'-b', t,YTC,'--g', t,QT5,'--r', t,Tr2_5,'--b')
grid ;
pause ;
plot( t,YT6_1,'-b', t,YT6_2,'-b', t,YT6_3,'-b', t,YT6_4,'-b', t,YT6_5,'-b', t,YT6_6,'-b', t,YT6_7,'-b', t,YT6_8,'-b', t,YT6_9,'-b', t,YT6_10,'-b', t,YTC,'--g', t,QT6,'--r', t,Tr2_6,'--b')
grid ;
pause ;
plot( t,QT1,'--r', t,Tr2_1,'--b', t,QT2,'--r', t,Tr2_2,'--b', t,QT3,'--r', t,Tr2_3,'--b', t,QT4,'--r', t,Tr2_4,'--b', t,QT5,'--r', t,Tr2_5,'--b', t,QT6,'--r', t,Tr2_6,'--b')
grid ;
pause ;
