clear
 clc
 load nodev.mat
 load mh.mat
 load Ibi.mat
 
 t = 0:0.005:100;
 %% geneate S4A Fig
 figure;
for i=1:21
    m(1:20001)=i-1;
    hold on
    plot3(m,t,nodev(:,i),'k');
end
zlabel('V (mV)');zlim([-110 100]);ylabel('t (ms)');xlabel('node #');xlim([-0.5 20.5]);view(35,80)

 %% geneate S4B Fig
 figure
 subplot(3,1,1)
 plot(t,nodev(:,11),'k');ylabel('V (mV)');
 subplot(3,1,2);
 plot(t,mh(:,1),'color',[0.85,0.33,0.10]);hold on; plot(t,mh(:,2),'color',[0.00,0.45,0.74]);ylabel('gating variables');ylim([-0.05 1.05]);
 subplot(3,1,3)
 plot(t,Ibi,'k');ylabel('I_B_I (mA)');xlabel('t (ms)');