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');