Someone could help me.
I post only a part of the programs
(excuse me if i made some mistake)
Code: Select all
//definition parametre du grand axone
create axon_grand
access axon_grand
axon_grand{
nseg=1
diam=18.8
L=18.8
Ra=123
insert hh
gnabar_hh=0.25
gl_hh=0.0001666
el_hh=-60
celsius=37
istim=2
delay=1
pw=0.1
axonnodes=21
paranodes1=40
paranodes2=40
axoninter=120
axontotal=221
fiberD_grand=5.6
paralength1=3
nodelength=1.0
space_p1=0.002
space_p2=0.004
space_i=0.004
rhoa=0.7e6
mycm=0.1
mygm=0.001}
//pas d'evolution
D = 0.1
//fonction pour faire evoluer
proc setD_grand(){
fiberD1_grand = fiberD_grand + D
fiberD_grand = fiberD1_grand}
//titre du tableau
xpanel("Fibres Grand Axon", 0)
//case pour regler le retard
xvalue("Diametre Grande Fibre", "fiberD_grand", 1, "setD_grand()", 0, 0)
//position de la fenetre
xpanel(38,890)
//case selon la taille de la fibre
if (fiberD_grand==5.7){
axon_grand{
g=0.605
axonD=3.4
nodeD=1.9
paraD1=1.9
paraD2=3.4
deltax=500
paralength2=35
nl=80
}
}
if (fiberD_grand==7.3){
axon_grand{
g=0.630
axonD=4.6
nodeD=2.4
paraD1=2.4
paraD2=4.6
deltax=750
paralength2=38
nl=100
}
}
if (fiberD_grand==8.7){
axon_grand{
g=0.661
axonD=5.8
nodeD=2.8
paraD1=2.8
paraD2=5.8
deltax=1000
paralength2=40
nl=110
}
}
if (fiberD_grand==10.0){
axon_grand{
g=0.690
axonD=6.9
nodeD=3.3
paraD1=3.3
paraD2=6.9
deltax=1150
paralength2=46
nl=120
}
}
if (fiberD_grand==11.5){
axon_grand{
g=0.700
axonD=8.1
nodeD=3.7
paraD1=3.7
paraD2=8.1
deltax=1250
paralength2=50
nl=130
}
}
if (fiberD_grand==12.8){
axon_grand{
g=0.719
axonD=9.2
nodeD=4.2
paraD1=4.2
paraD2=9.2
deltax=1350
paralength2=54
nl=135
}
}
if (fiberD_grand==14.0){
axon_grand{
g=0.739
axonD=10.4
nodeD=4.7
paraD1=4.7
paraD2=10.4
deltax=1400
paralength2=56
nl=140
}
}
if (fiberD_grand==15.0){
axon_grand{
g=0.767
axonD=11.5
nodeD=5.0
paraD1=5.0
paraD2=11.5
deltax=1450
paralength2=58
nl=145
}
}
if (fiberD_grand==16.0){
axon_grand{
g=0.791
axonD=12.7
nodeD=5.5
paraD1=5.5
paraD2=12.7
deltax=1500
paralength2=60
nl=150
}
}