%% Carol Lucas
%% Demonstrates Figure 5.6A for full model with jumps in C

fullinit=[.3;.2;.2];poinit=[.4];cinit=.35;cfinal=.7;

%% No jump late (22 secs)


[ttotnj,xnj]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,22);
Pototnj=xnj(:,2)+xnj(:,3);
figure(1);clf;plot(ttotnj,Pototnj,'b');axis([0 20 0 .5]);
xlabel('time(sec)');ylabel('Po = O1 + O2');hold on;
mess = 'No visible jump, hit return to continue',pause;

%% Jump at 2 seconds

[ttot2,x2]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,2);
Potot2=x2(:,2)+x2(:,3);
plot(ttot2,Potot2,'g');
mess= 'Jump at 2 seconds, hit return to continue',pause;

%% Jump at 8 seconds

[ttot8,x8]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,8);
Potot8=x8(:,2)+x8(:,3);
plot(ttot8,Potot8,'r');
mess='Jump at 8 seconds, hit return to finish'