m1=0.54; m2=0.24; m3=0.2; m4=0.2; m5=0.3; m6=0.3; m7=0.7; m8=0.4; m9=1.1; m10=1; m11=1; m12=1; m13=0.32; m14=0.4; m15=0.7; 
m16=0.5; m17=0.5; m18=3.4; m19=0.2; m20=0.6; m21=0.08; m22=0.1; m23=1.8; m24=0.1; m25=1.8; m26=0.5; m27=0.1; m28=20; m29=5; 
m30=3; m31=0.3; m32=0.2; m33=13; m34=0.6; m35=0.3; m36=0.1; m37=0.8; m38=0.5; m39=0.3; 

n1=2.6; n2=0.64; n3=0.29; n4=0.07; n5=0.23; n6=20; n7=0.2; n8=0.5; n9=0.2; n10=0.4; n11=0.6; n12=12.5; n13=1.3; n14=0.1; n15=0.4; n16=0.23; n17=0.1; 

p1=0.13; p2=0.27; p3=0.1; p4=0.56; p5=4; p6=0.6; p7=0.3; p8=0.6; p9=0.8; p10=0.54; p11=0.51; p12=3.4; p13=0.1; p14=0.14; p15=3;
p16=0.62; p17=4.8; p18=4; p19=1; p20=0.1; p21=1; p22=0.5; p23=0.37; p24=10; p25=8; p26=0.3; p27=0.8; p28=2; p29=0.1; p30=0.9; p31=0.1;

g1=0.1; g2=0.01; g3=0.6; g4=0.01; g5=0.15; g6=0.3; g7=0.6; g8=0.01; g9=0.3; g10=0.5; g11=0.7; g12=0.2; g13=1; g14=0.004; g15=0.4;
g16=0.3; g17=0.03; g18=0.3; g19=0.5; g20=1.1; g21=0.5; g22=1.5; g23=0.01; g24=3; g25=0.03; g26=0.3; g27=0.1; g28=0;

a=2; b=2; c=2; d=2; e=2; f=2; g=0; h=0; 

q1=1.2; q2=1.56; q3=2.8; q4=0; 

period=24;
dusk=12;
day_numb=5;

t = [0 day_numb*period];

%Solving the ODEs and draw traectories
options=odeset('stat','off');

y0=[1.0151 0.956 0.0755 0.0041 0.506 0.0977 0.0238 0.0731 0.0697 0.0196 0.0435 0.2505 0.0709 0.1017 0.0658 0.4016 0.1167 0.1012 0.207 0.0788 0.3102 0.0553 0.2991 0.1503 0.0286 0.65 0.2566 0.1012 0.576 0.3269]; %12:12 wt

[T, Y] = ode23s(@Model,t,y0,options,q1,q2,q3,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16,p17,p18,p19,p20,p21,p22,p23,p24,p25,p26,p27,p28,p29,p30,p31,m1,m2,m3,m4,m5,m6,m7,m8,m9,m10,m11,m12,m13,m14,m15,m16,m17,m18,m19,m20,m21,m22,m23,m24,m25,m26,m27,m28,m29,m30,m31,m32,m33,m34,m35,m36,m37,m38,m39,n1,n2,n3,n4,n5,n6,n7,n8,n9,n10,n11,n12,n13,n14,g1,g2,g3,g4,g5,g6,g7,g8,g9,g10,g11,g12,g13,g14,g15,g16,a,b,c,d,e,f,dusk,period);

%T=T-24*3;

GN=p28*Y(:,10)./(p29+m19+p17*Y(:,25));
EGN=(p18*Y(:,4)+p17*Y(:,25).*GN)./(m38*Y(:,26)+m39*Y(:,27)+p31);
E34=p25*Y(:,19).*Y(:,25)./(p26*Y(:,29)+p21+m36*Y(:,26)+m37*Y(:,27));

figure (1)
plot(T,Y(:,6),'r');
hold on;
plot(T,Y(:,8),'b');
hold on;
plot(T,Y(:,16),'g');
hold on;
plot(T,Y(:,15),'m');
hold on;
plot(T,Y(:,1),'k');
hold on;
title('LHY mRNA-black, PRR9-magenta; PRR7-green; NI-blue; protein P-dot red');

figure (2)
plot(T,Y(:,28),'b');
hold on;
plot(T,Y(:,13),'g');
hold on;
title('LUX mRNA-blue; EC-green');