FindIth_Iext_axon.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=50;
delay=20;
Z=0;
Iexstim([0 0 Z],-1,pw,00) %(xyz(um),amp(uA),dur(us),delay(us))