%% F5_16:  Call for Figure 5_16
% The initial conditions
V=-36;
n=0.8;
C=0.16;
h=0.8;
Ct=1.6;
sigma=.185;
[t,x]=ode45('lrpmlfun',[0 30],[V;n;C;h;Ct]);

figure(1);clf;subplot(411);plot(t,x(:,1));
grid on 
ylabel('V', 'FontSize', 20);
subplot(412);plot(t,x(:,3 ));
grid on 
ylabel('Ca inside', 'FontSize', 20);
Ce=(x(:,5)-x(:,3))/sigma;
subplot(413);plot(t,Ce);
grid on 
ylabel('Ca ER', 'FontSize', 20);
subplot(414);plot(t,x(:,5));
grid on 
xlabel('time (s)', 'FontSize', 20);ylabel('Ca total', 'FontSize', 20);
grid on 

%%%

figure(2)
clf
open PlasmaMembraneBifDiag_color.fig 

hold on;
plot(x(:,3 ), x(:,1), 'r:')