Ex_07_Iext_axon_with_myelin.m

create('axon');

%-------------------------------------------------------------------------------------------
global axon
%-------------------------------------------------------------------------------------------

axon.nseg		=55;                                    %       change the number of segments
axon.dia		=10;                                    % um    diameter
axon.Lmode		=2;
axon.Linter_D	=100;
axon.L   		=axon.Linter_D*axon.dia*(axon.nseg);    %um     change the length of axon
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')
Iexstim([ 10000 0 500],-1200,500,10)                         % [xyz(um),amp(uA),dur(us),delay(us)]
%instim(axon,1,0.001,500,20000)                    %intracellular stimulation     [variable positon(0-1) amp(uA) pw(us)  delay(us)]
%Iexstim([-5000 0 500],'Isin')                         % [xyz(um),amp(uA),dur(us),delay(us)]
%Multiple stimulators supported
if 0
Iexstim([ 10000 0 500],-1200,500,10)                         % [xyz(um),amp(uA),dur(us),delay(us)]
Iexstim([-5000 0 500],'Isin')                         % [xyz(um),amp(uA),dur(us),delay(us)]
Iexstim([ 5000 0 500],'SquarePulse')                         % [xyz(um),amp(uA),dur(us),delay(us)]
end

%The parameters that are used in Neurocal are
%.name			%  	    internal use.
%.varnum        %  	    internal use.
%.nseg			% 	    number of segments for the variable.
%.dia       	%(um)	the outer diameter of the variable.
%.da_D  		%		axon/Fiber diameter ratio for calculating the axonal resistance 
                %         - for myelinated axon. if this value is empty(default), 
%                         - .dia will be used instead.
%.dn_D      	%	    node/Fiber diameter ratio for node of Ranvior for calculating the membrane resistance
                %         - and membrane capacitance for myelinated axon. If this value is empty(default), 
                %         - .dia will be used instead.
%.Linter_D      %       Internodal distance / Diameter                
%.L  		    %(um)	the total length of the variable.
%.Lmode         %       1=fix=L, 2 = Linter_D x D x nseg
%.lnodal      	%(um)	Lengh of node of Ranvior for calculating the membrane resistance
                %         - and membrane capacitance for myelinated axon. If this value is empty(default), 
                %         - L/nseg will be used instead.
%.cm	       	%(uf/cm2)for calculating membrane capacitance
%.ra           	%(ohm-cm)for calculating axonal resistance
%.rm           	%(kohm-cm2)for calculating membrane resistance.
%.vini	    	%(mV)
%.model         %
%.vestim		%(mV)   Extracellular voltage (when Ve is not calculated,
%                         - but assigned manually
%.xyz			%(um)   Coordinate of the nodes used for calculating the Ve
%.xyzi    	    % unitless: Vector of the variable (axon) for calculating the coordinate for all nodes
%.xyzc          % um:       The coordinate of the node at the center