The code below is from McIntyre's model(2002) and I add some comments to make it easier to read.
Code: Select all
for i=0, paranodes1-1 {
MYSA[i]{
nseg=1
diam=fiberD
L=paralength1
Ra=rhoa*(1/(paraD1/fiberD)^2)/10000
cm=2*paraD1/fiberD
insert pas
g_pas=0.001*paraD1/fiberD
e_pas=v_init
insert extracellular xraxial=Rpn1 xg=mygm/(nl*2) xc=mycm/(nl*2) //nl:number of myelin lemella, mygm:myelin conductance
}
1.Why 'xg=mygm/(nl*2)'? As there are nl layers of transverse resistive in series. I think it should be mygm/nl.
2.What does 'xraxial' means and what's the unit of it?
Please tell me.
Thank you!