%% F5_7A
%% Generates nullclines as shown in Figure 5.7
%% Calls klclosejac--which does each analysis separately
%% Then all will be superimposed in Figure 4
%% Carol Lucas


CTRAY=[1.8 2.2 3.2];

for k=1:3
    
CT=CTRAY(k);
figure(k);clf;
[Cs,Ws,lams,CRAY,nullC,nullw]=klclosejac(CT);
nullCRAY(k,:)=nullC;
end

figure(4);clf;plot(CRAY,nullw,CRAY,nullCRAY);axis([0 3 0 3]);
xlabel('[Ca2+]');ylabel('w');hold on;

mess='hit return to see some simulations around these values'
pause;

klclose_anal(1,1,3.2);
klclose_anal(2.5,2,2.2)
klclose_anal(.5,.5,2.2)
klclose_anal(.1,.1,2.2)
klclose_anal(.7,.2,2.2)
klclose_anal(1,1,1.8)