function Figure5_4

close all;
boxes=40;

load dataFigure5_4.dat;
dataSTO=dataFigure5_4;

morblue=zeros(boxes,1);
morred=zeros(boxes,1);
morblue2=zeros(boxes,1);
morred2=zeros(boxes,1);

for ii=1:boxes
    if (dataSTO(ii,2)>200)
       morblue(ii)=dataSTO(ii,2);
    else
       morred(ii)=dataSTO(ii,2);
    end;
    if (dataSTO(ii,3)>75)
       morblue2(ii)=dataSTO(ii,3);
    else
       morred2(ii)=dataSTO(ii,3);
    end;
end;

figure(1);
set(gca,'Fontsize',20);
bar(dataSTO(:,1),morblue(:,1),'b');
hold on;
bar(dataSTO(:,1),morred(:,1),'r');
xlabel('x [mm]');
ylabel('number of A molecules');
text(0.027,476,'time=30 min','Fontsize',20);
h=line([0 1],[200 200]);
set(h,'Color','k','LineStyle','--');
axis([0 1 0 510]);

figure(2);
set(gca,'Fontsize',20);
bar(dataSTO(:,1),morblue2(:,1),'b');
hold on;
bar(dataSTO(:,1),morred2(:,1),'r');
xlabel('x [mm]');
ylabel('number of B molecules');
text(0.027,112,'time=30 min','Fontsize',20);
h=line([0 1],[75 75]);
set(h,'Color','k','LineStyle','--');
axis([0 1 0 120]);
