%% Program to generate figure 5.8
%% Carol Lucas

[t,x]=ode45('klopenfun_red',[0:500],[.2;2.6],[],1.5,50,400,.1);
figure(1);clf;plot(t,x(:,1));
xlabel('time (S)');ylabel('[Ca2+]i');
title('Figure 5.8 - Reduced open-cell Keizer-Levine, jin=1.5 from 50 to 400 secs');