%% Carol Lucas
%% Demonstrates Figure 5.6 a for full compared to reduced receptor model 
%% Part of Problem 3C
%% Recompiled with different tjump values
%% Calls klfun and klfunPo

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

%%Jump at t=2
[ttotnj,xnj]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,22);
Pototnj=xnj(:,2)+xnj(:,3);
[ttot,x]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,2);
Potot=x(:,2)+x(:,3);
[tponj,Ponnj]=ode45('klfunPo',[0 20],poinit,[],cinit,cfinal,22);
[tpo,Po]=ode45('klfunPo',[0 20],poinit,[],cinit,cfinal,2);
figure(1);clf;plot(ttotnj,Pototnj,tponj,Ponnj,ttot,Potot,tpo,Po,'m+:');axis([0 20 0 .5]);
xlabel('time(sec)');ylabel('Po');title('Full vs. Reduced, Jump at 2');
mess= ' Jump at 2 seconds, hit return to continue', pause;

%% Jump at t=8
[ttot,x]=ode45('klfun',[0 20],fullinit,[],cinit,cfinal,8);
Potot=x(:,2)+x(:,3);
[tpo,Po]=ode45('klfunPo',[0 20],poinit,[],cinit,cfinal,8);
figure(2);clf;plot(ttotnj,Pototnj,tponj,Ponnj,ttot,Potot,tpo,Po,'m+:');axis([0 20 0 .5]);
xlabel('time(sec)');ylabel('Po');title('Full vs. Reduced, Jump at 8');
mess='Jump at 8 seconds, hit return to finish', pause;