clear clc load vpso_DC.mat load vnode10_KHFAC.mat load IPSO_DC.mat load IPSO_KHFAC.mat load inav11_DC.mat load inav11_KHFAC.mat load nodev_KHFAC.mat t = 0:0.005:100; %% geneate Fig 4A figure; subplot(3,1,1);plot(t,vnode10_KHFAC,'k');hold on;plot(t,vpso,'r');ylabel('V (mV)') subplot(3,1,2);plot(t,inav11_KHFAC(:,1),'k');hold on;plot(t,inav11_DC(:,1),'r');ylabel('I_N_a_v_1_1 (mA/cm^2)'); subplot(3,1,3);plot(t,1.5*IPSO_KHFAC,'k');hold on;plot(t,IPSO,'r');ylabel('stimulus (mA)');xlabel('t (ms)') %% geneate Fig 4B figure; subplot(3,1,1);plot(t,inav11_KHFAC(:,2),'k');hold on;plot(t,inav11_DC(:,2),'r');ylabel('C11') subplot(3,1,2);plot(t,inav11_KHFAC(:,3),'k');hold on;plot(t,inav11_DC(:,3),'r');ylabel('I12') subplot(3,1,3);plot(t,inav11_KHFAC(:,4),'k');hold on;plot(t,inav11_DC(:,4),'r');ylabel('O11');xlabel('t (ms)') %% geneate Fig 4C figure; t=0:0.005:100; for i=1:21 m(1:20001)=i-1; hold on plot3(m,t,nodev(:,i),'k'); end zlabel('V (mV)');ylabel('t (ms)');xlabel('node #');xlim([-0.5 20.5]);view(7,47)