No BSD License  

Highlights from
Neurocal

image thumbnail
from Neurocal by Zeng Lertmanorat
Simulation describing the electrical activity of nerve cell (neuron) by solving cable equation

Ex_08_Simple_Neuron.m
%-------------------------------------------------------------------------------------------
create('axon');
create('soma');
create('dendrite',3)            %create dendrite(1),dendrite(2),dendrite(3)

global axon soma dendrite

axon.L=2000;
axon.nseg=21;

dendrite{1}.L=2000;
dendrite{1}.nseg=20;
dendrite{2}.L=4000;
dendrite{2}.nseg=20;
dendrite{3}.L=2000;
dendrite{3}.nseg=20;

connect(soma,0,dendrite{1},0); %[main, position(0-1), branch, position (0 OR 1)]
connect(soma,0,dendrite{2},0); %the branch position is either 0 or 1
connect(soma,0,dendrite{3},0); %the branch position is either 0 or 1
connect(axon,0,soma,0);        %the branch position is either 0 or 1

insert(axon,'hh')
Iinstim(soma,1,0.05,500,0) %[variable positon(0-1) amp(uA) pw(us)  delay(us)]

%-------------------------------------------------------------------------------------------
%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.
%.Lmode         %   1=fix, 2 = L/D x D x nseg
%.Linter_D      %       - Internodal distance / Diameter                
%.L  		    %	um	- the total length of the variable.
%.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
%.vestim		%	mV  -extracellular stimulator
%.xyz			%	um
%.xyzvector		%	um

Contact us at files@mathworks.com