%% DEMONSTRATES Chapter 6 Figures 2 and 3
%% Calls 'mlgapfun' and fanal to get frequency
%% Carol Lucas

%% Program pauses after each figure.  Hit return to continue.'

initgap=[-20;.2;-30;0];

gc=2;   %% GAPS STRONG

Titlegc=['  gc =  ' num2str(gc)];

[t,x]= ode45('mlgapfun',[0:.1:250],initgap,[],gc);
figure(1);clf;plot(t,x(:,1),t,x(:,3)); 
title(Titlegc); xlabel('time');ylabel('Voltages');pause;
beat1V1=fanal(t,x(:,1)); 
figure(2);clf;plot(t,x(:,1),t,beat1V1);
title([Titlegc ['  Cell 1']]); xlabel('time');ylabel('Voltages & Inter Spike Intervals');pause;
beat1V2=fanal(t,x(:,3));
plot(t,x(:,3),t,beat1V2);
title([Titlegc  ['  Cell 2']]); xlabel('time');ylabel('Voltages & Inter Spike Intervals');pause;
figure(2);plot(t,beat1V1,t,beat1V2);
title(Titlegc); xlabel('time');ylabel('Inter Spike Intervals Comparisons');pause;

gc=1;     %%  GAPS WEAK

Titlegc=['  gc =  ' num2str(gc)];

[t,x]= ode45('mlgapfun',[0:.5:250],initgap,[],gc);
figure(3);plot(t,x(:,1),t,x(:,3));
title(Titlegc); xlabel('time');ylabel('Voltages');pause;
beat1V1=fanal(t,x(:,1));
beat1V2=fanal(t,x(:,3));
figure(4);plot(t,beat1V1,t,beat1V2);
title(Titlegc); xlabel('time');ylabel('Inter Spike Intervals Comparisons');pause;