clear
 clc
 load vstep.mat
 load inav11_step.mat
 load inav16_step.mat
 load nodev_step.mat
 load vpso.mat
 load inav11_pso.mat
 load inav16_pso.mat
 load nodev_pso.mat
 
 t = 0:0.005:100;
%% geneate Fig 2A
 figure;
 subplot(3,1,1);plot(t,vstep,'k');ylabel('V_s_t_e_p (mV)')
 subplot(3,1,2);plot(t,inav11_step(:,2),'k');ylabel('I_N_a_v_1_1 (mA/cm^2)')
 subplot(3,1,3);plot(t,inav16_step(:,2),'k');ylabel('I_N_a_v_1_6 (mA/cm^2)');xlabel('t (ms)')
%% geneate Fig 2B
 figure;
 subplot(2,1,1);plot(t,inav11_step(:,3),'k');hold on;plot(t,inav11_step(:,5),'b');plot(t,inav11_step(:,6),'r');plot(t,inav11_step(:,7),'g'); ylabel('fraction')
 subplot(2,1,2);plot(t,inav16_step(:,3),'k');hold on;plot(t,inav16_step(:,5),'b');plot(t,inav16_step(:,6),'r');plot(t,inav16_step(:,7),'g'); ylabel('fraction');xlabel('t (ms)')
%% geneate Fig 2C
 figure;
 t=0:0.005:100;
for i=1:21
    m(1:20001)=i-1;
    hold on
    plot3(m,t,nodev_step(:,i),'k');
end
zlabel('V (mV)');ylabel('t (ms)');xlabel('node #');xlim([-0.5 20.5]);view(35,26)
%% geneate Fig 2D
 figure;
 subplot(3,1,1);plot(t,vpso,'k');ylabel('V_P_S_O (mV)')
 subplot(3,1,2);plot(t,inav11_pso(:,2),'k');ylabel('I_N_a_v_1_1 (mA/cm^2)')
 subplot(3,1,3);plot(t,inav16_pso(:,2),'k');ylabel('I_N_a_v_1_6 (mA/cm^2)');xlabel('t (ms)')
  %% geneate Fig 2E
 figure;
 subplot(2,1,1);plot(t,inav11_pso(:,3),'k');hold on;plot(t,inav11_pso(:,5),'b');plot(t,inav11_pso(:,6),'r');plot(t,inav11_pso(:,7),'g'); ylabel('fraction')
 subplot(2,1,2);plot(t,inav16_pso(:,3),'k');hold on;plot(t,inav16_pso(:,5),'b');plot(t,inav16_pso(:,6),'r');plot(t,inav16_pso(:,7),'g'); ylabel('fraction');xlabel('t (ms)')
  %% geneate Fig 2F
 figure;
 t=0:0.005:100;
for i=1:21
    m(1:20001)=i-1;
    hold on
    plot3(m,t,nodev_pso(:,i),'k');
end
zlabel('V (mV)');ylabel('t (ms)');xlabel('node #');xlim([-0.5 20.5]);view(35,26)