%% F5_19   Generates figures as in 5_19 and 5_20
%% Calls mlbetafun for various values of VLPM
%% VLPMRAY=[.10 .15 .18 .20 .50];
%% Carol Lucas

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

vinit=-65.0; ninit=0.00016; cinit=0.2;
VLPMRAY=[.10 .15 .18 .20 .50];

for k=1:length(VLPMRAY)
    VLPM=VLPMRAY(k);
    [t,x]=ode45('mlbetafun',[0:30000], [vinit;ninit;cinit],[],VLPM);
    figure(1);clf;subplot(211);plot(t,x(:,1));ylabel('V');title(['Vlpm =' num2str(VLPM)]);
    subplot(212);plot(t,x(:,3));xlabel('time (ms)');ylabel('[Ca2+]i');
    pause;
    figure(2);clf;plot(x(:,3),x(:,1));
    ylabel('[Ca2+]i');xlabel('V');title(['Vlpm =' num2str(VLPM)]);axis([0 .4 -70 -10]);pause;
    
end