comCoEx3.m

clear;
global TayOS; TayOS = 20 ;
format short e;
format compact;
global index;
index=0; 

StTqlo = [ 1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;   1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;   1300 ;1300 ;1300;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;   1300 ;1300 ; 1300; 1300; 1300; 1300; 1300; 1300; 1300; 1300;   1300 ;1300 ;1300 ;1300 ;1300 ;1300 ; 1300; 1300; 1300; 1300;   1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;1300 ;     1300;];


[t,y]=ode15s (@CoEx3,[0 600000],[StTqlo;]); 

long_m = length(t)
    YT1_1 = y(1:long_m, 1); 
    YT1_2 = y(1:long_m, 2); 
    YT1_3 = y(1:long_m, 3); 
    YT1_4 = y(1:long_m, 4); 
    YT1_5 = y(1:long_m, 5); 
    YT1_6 = y(1:long_m, 6); 
    YT1_7 = y(1:long_m, 7); 
    YT1_8 = y(1:long_m, 8); 
    YT1_9 = y(1:long_m, 9); 
    YT1_10 = y(1:long_m, 10); 
    YT2_1 = y(1:long_m, 11); 
    YT2_2 = y(1:long_m, 12); 
    YT2_3 = y(1:long_m, 13); 
    YT2_4 = y(1:long_m, 14); 
    YT2_5 = y(1:long_m, 15); 
    YT2_6 = y(1:long_m, 16); 
    YT2_7 = y(1:long_m, 17); 
    YT2_8 = y(1:long_m, 18); 
    YT2_9 = y(1:long_m, 19); 
    YT2_10 = y(1:long_m, 20); 
    YT3_1 = y(1:long_m, 21); 
    YT3_2 = y(1:long_m, 22); 
    YT3_3 = y(1:long_m, 23); 
    YT3_4 = y(1:long_m, 24); 
    YT3_5 = y(1:long_m, 25); 
    YT3_6 = y(1:long_m, 26); 
    YT3_7 = y(1:long_m, 27); 
    YT3_8 = y(1:long_m, 28); 
    YT3_9 = y(1:long_m, 29); 
    YT3_10 = y(1:long_m, 30); 
    YT4_1 = y(1:long_m, 31); 
    YT4_2 = y(1:long_m, 32); 
    YT4_3 = y(1:long_m, 33); 
    YT4_4 = y(1:long_m, 34); 
    YT4_5 = y(1:long_m, 35); 
    YT4_6 = y(1:long_m, 36); 
    YT4_7 = y(1:long_m, 37); 
    YT4_8 = y(1:long_m, 38); 
    YT4_9 = y(1:long_m, 39); 
    YT4_10 = y(1:long_m, 40); 
    YT5_1 = y(1:long_m, 41); 
    YT5_2 = y(1:long_m, 42); 
    YT5_3 = y(1:long_m, 43); 
    YT5_4 = y(1:long_m, 44); 
    YT5_5 = y(1:long_m, 45); 
    YT5_6 = y(1:long_m, 46); 
    YT5_7 = y(1:long_m, 47); 
    YT5_8 = y(1:long_m, 48); 
    YT5_9 = y(1:long_m, 49); 
    YT5_10 = y(1:long_m, 50); 
    YT6_1 = y(1:long_m, 51); 
    YT6_2 = y(1:long_m, 52); 
    YT6_3 = y(1:long_m, 53); 
    YT6_4 = y(1:long_m, 54); 
    YT6_5 = y(1:long_m, 55); 
    YT6_6 = y(1:long_m, 56); 
    YT6_7 = y(1:long_m, 57); 
    YT6_8 = y(1:long_m, 58); 
    YT6_9 = y(1:long_m, 59); 
    YT6_10 = y(1:long_m, 60); 
	 YTC = y(1:long_m, 61); 

 save my_file_ode

plot(t, y)
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' )
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' )
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' )
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' )
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' )
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' )
grid ; 
pause ;