%% Figure 6_10
%% Calls ode function 'mlds2cellfun'
%% Carol Lucas

clear tdelray;
clear vmax1;
clear vmax2;
figure(1);clf;

v1init=-60.9;w1init=.0149;seinit=0.;siinit=0.;
v2init=-60.9;w2init=.0149;s21init=0.;tdelinit=0;se2init=0.;

init=[v1init;w1init;v2init;w2init;seinit;siinit;se2init;s21init];

tdelray=[-40:5:-15 -14:14 15:5:40];

for k=1:length(tdelray)
    tdel=tdelray(k);
[t,x]=ode23('mlds2cellfun',[0 200],init,[],tdel);
figure(1);plot(t,x(:,1),t,x(:,3));
xlabel('time (ms)');ylabel(' Voltages, V1 and V2');title(['tdel = ' int2str(tdel)]);pause(1);

vmax1(k)=max(x(:,1));
vmax2(k)=max(x(:,3));

end

figure(2);clf;plot(tdelray,vmax1,'bo-',tdelray,vmax2,'g+:');
xlabel('time delay');ylabel('V maxes of V2');