correcting color mismatch in subfig.3 and initialstate/initial_state name error in fig.5 script

reduced_model_figure5.m

@@ -17,12 +17,12 @@ x_limit=50;  y_limit_tnf=6000; y_limit_il1b=70; y_lmit_il10=130;
 color_IL10_0h=[240 230 140]./255;    color_control=[9 121 105]./255; 
 color_IL10_4h=[252 141 89]./255;     color_LPS_block=[25 110 170]./255;
 cytokine_num=4; % number of the measured cytokines
-initialstate = zeros(1,10); simtime=0:0.1:50; 
+initial_state = zeros(1,10); simtime=0:0.1:50; 
 %% Uncertainty of predictions by the reduced model
 Xt=selectedX; 
 max=zeros(cytokine_num,length(simtime));  min=inf*ones(cytokine_num,length(simtime));
 for h=1:size(Xt,1)
-    IC =ALLFUNCS.SteadyState(model,Xt(h,:),initialstate); 
+    IC =ALLFUNCS.SteadyState(model,Xt(h,:),initial_state); 
     out1= ALLFUNCS.case2sim(IC,model,simtime,Xt(h,:));
     k=0;
     for j=[5,4,6,7]
@@ -177,29 +177,29 @@ for h=1:size(Xt,1)
 end
 %% TNF uncertainty - "IL-10-Block" - reduced model
 figure(3); subplot(2,2,1);
-plot (simtime,min(1,:),'linewidth',2,'color',color_control); hold on;
-plot (simtime,max(1,:),'linewidth',2,'color',color_control); 
+plot (simtime,min(1,:),'linewidth',2,'color',color_LPS_block); hold on;
+plot (simtime,max(1,:),'linewidth',2,'color',color_LPS_block); 
 p1=patch([simtime fliplr(simtime)], [min(1,:) fliplr(max(1,:))],'y','FaceAlpha',0.5)
 p1.EdgeColor=color_LPS_block; p1.FaceColor=color_LPS_block; hold on;
 ax=gca; ax.YAxis.Exponent = 3; xticks([0,10,20,30,40,50]);
 %% IL-10 uncertainty - "IL-10-Block" - reduced model
 figure(3); subplot(2,2,2);
-plot (simtime,min(2,:),'linewidth',2,'color',color_control); hold on;
-plot (simtime,max(2,:),'linewidth',2,'color',color_control); 
+plot (simtime,min(2,:),'linewidth',2,'color',color_LPS_block); hold on;
+plot (simtime,max(2,:),'linewidth',2,'color',color_LPS_block); 
 p2=patch([simtime fliplr(simtime)], [min(2,:) fliplr(max(2,:))],'y','FaceAlpha',0.5)
 p2.EdgeColor=color_LPS_block; p2.FaceColor=color_LPS_block; hold on;
 ax=gca; ax.YAxis.Exponent = 1; xticks([0,10,20,30,40,50]);
 %% IL-1b uncertainty - "IL-10-Block" - reduced model
 figure(3); subplot(2,2,3);
-plot (simtime,min(3,:),'linewidth',2,'color',color_control); hold on;
-plot (simtime,max(3,:),'linewidth',2,'color',color_control);
+plot (simtime,min(3,:),'linewidth',2,'color',color_LPS_block); hold on;
+plot (simtime,max(3,:),'linewidth',2,'color',color_LPS_block);
 p3=patch([simtime fliplr(simtime)], [min(3,:) fliplr(max(3,:))],'y','FaceAlpha',0.5)
 p3.EdgeColor=color_LPS_block; p3.FaceColor=color_LPS_block; hold on;
 ax=gca; ax.YAxis.Exponent = 1; xticks([0,10,20,30,40,50]);
 %% IL-1Ra uncertainty - "IL-10-Block" - reduced model
 figure(3); subplot(2,2,4);
-plot (simtime,min(4,:),'linewidth',2,'color',color_control); hold on;
-plot (simtime,max(4,:),'linewidth',2,'color',color_control); 
+plot (simtime,min(4,:),'linewidth',2,'color',color_LPS_block); hold on;
+plot (simtime,max(4,:),'linewidth',2,'color',color_LPS_block); 
 p4=patch([simtime fliplr(simtime)], [min(4,:) fliplr(max(4,:))],'y','FaceAlpha',0.5)
 p4.EdgeColor=color_LPS_block; p4.FaceColor=color_LPS_block; hold on;
 ax=gca; ax.YAxis.Exponent = 3; xticks([0,10,20,30,40,50]);