FindIth_Iext_axon_C2A3.m

create('axon');
%-------------------------------------------------------------------------------------------
global axon
%-------------------------------------------------------------------------------------------
axon.nseg		=73;                           %       change the number of segments
axon.dia		=10;                            % um    diameter
axon.Lmode		=2;
axon.Linter_D	=100;
axon.da_D 		=0.6;
axon.lnodal		=1;
axon.dn_D		=0.33;
axon.xyzi   	=[1 0 0];                                   % unitless
axon.xyzc       =[0 0 0];                                   % um: the coordinate of the node at the center
insert(axon,'Schwarz')
%-------------------------------------------------------------------------------------------
pw=10;
delay=00;
Z=200;
Iexstim([-1500 0 Z], 1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us)) 
Iexstim([ -750 0 Z],-1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us)) 
Iexstim([    0 0 Z], 1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us)) 
Iexstim([  750 0 Z],-1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us)) 
Iexstim([ 1500 0 Z], 1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us))