%% 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(311);plot(t,x(:,1));ylabel('V');
subplot(312);plot(t,x(:,3 ));ylabel('Ca inside');
Ce=(x(:,5)-x(:,3))/sigma;
subplot(313);plot(t,Ce);xlabel('time (s)');ylabel('Ca ER');