clear
 clc
 load nodev_single_47.mat
 load nodev_single_119.mat
 load nodev_200Hz_47.mat
 load nodev_200Hz_119.mat
 load ina_119_single.mat
 load ina_47_single.mat
 load ina_119_200Hz.mat
 load ina_47_200Hz.mat
 
 t = 0:0.005:200;
 %% geneate S1A Fig, left panels
 figure;title('single test pulse')
for i=1:21
    m(1:40001)=i-1;
    hold on
    plot3(m,t,nodev_single_47(:,i),'color',[0.00,0.45,0.74]);
end
zlabel('V (mV)');ylabel('t (ms)');ylim([99 105]);xlabel('node #');xlim([-0.5 20.5]);view(12,18)
figure
subplot(2,2,1)
plot(t,nodev_single_47(:,4),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,2)
plot(t,nodev_single_47(:,11),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,3)
plot(t,ina_47_single(:,1),'r');ylabel('I_N_a_f (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,4)
plot(t,ina_47_single(:,2),'r');ylabel('I_N_a_v_1_1 (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);

 %% geneate S1A Fig, right panels
 figure;title('200 Hz test pulses')
for i=1:21
    m(1:40001)=i-1;
    hold on
    plot3(m,t,nodev_200Hz_47(:,i),'color',[0.00,0.45,0.74]);
end
zlabel('V (mV)');ylabel('t (ms)');ylim([99 105]);xlabel('node #');xlim([-0.5 20.5]);view(12,18)
figure
subplot(2,2,1)
plot(t,nodev_200Hz_47(:,4),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,2)
plot(t,nodev_200Hz_47(:,11),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,3)
plot(t,ina_47_200Hz(:,1),'r');ylabel('I_N_a_f (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,4)
plot(t,ina_47_200Hz(:,2),'r');ylabel('I_N_a_v_1_1 (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);


 %% geneate S1B Fig, left panels
 figure;title('single test pulse')
for i=1:21
    m(1:40001)=i-1;
    hold on
    plot3(m,t,nodev_single_119(:,i),'color',[0.00,0.45,0.74]);
end
zlabel('V (mV)');ylabel('t (ms)');ylim([99 105]);xlabel('node #');xlim([-0.5 20.5]);view(12,18)
figure
subplot(2,2,1)
plot(t,nodev_single_119(:,4),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,2)
plot(t,nodev_single_119(:,11),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,3)
plot(t,ina_119_single(:,1),'r');ylabel('I_N_a_f (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,4)
plot(t,ina_119_single(:,2),'r');ylabel('I_N_a_v_1_1 (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);

 %% geneate S1B Fig, right panels
 figure;title('200 Hz test pulses')
for i=1:21
    m(1:40001)=i-1;
    hold on
    plot3(m,t,nodev_200Hz_119(:,i),'color',[0.00,0.45,0.74]);
end
zlabel('V (mV)');ylabel('t (ms)');ylim([99 105]);xlabel('node #');xlim([-0.5 20.5]);view(12,18)
figure
subplot(2,2,1)
plot(t,nodev_200Hz_119(:,4),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,2)
plot(t,nodev_200Hz_119(:,11),'k');ylabel('V (mV)');ylim([-99 50]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,3)
plot(t,ina_119_200Hz(:,1),'r');ylabel('I_N_a_f (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);
subplot(2,2,4)
plot(t,ina_119_200Hz(:,2),'r');ylabel('I_N_a_v_1_1 (mA/cm^2)');ylim([-60 0]);xlabel('t (ms)');xlim([99.8 102]);