%%  F5_14 Simulates 5_14 and Quick and Dirty Versions of 5_15
%%  Simulates Figure 5.14 and 5.15 A and B
%%  Carol Lucas

mess=' Program pauses after each figure.  Hit return to continue.'

[t,x]=ode45('lropenfun',[0 300],[.2;.8;4]);
figure(1);clf;plot(t,x);
xlabel('time (s)'); pause;
figure(2);clf;plot(x(:,3),x(:,1));
xlabel('[Ca2+]T');ylabel('[Ca2+]i');title('Phase Plane - F5.15A');pause;
[beat,tp,xbeat]=fanal(t,x(:,1));
figure(3);clf;plot(x(:,3),beat);
xlabel('[Ca2+]T');ylabel('Periods(s)');
axis([1.5 4 0 30]);title('Period vs [CA2+]T - F5.15B');pause;