XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# HH-Vext.ode: an XppAut file for the Hodgkin-Huxley equations
# Ref: Chen et al., Chinese J Physiol 2010, 54 (6), 423-429.
# + strength of external electrical field
# + trp channel
# time in msec, voltage in mV and frequency in Hz
# initial conditions
init v=-71, m=0.038, h=0.69, n=0.159
# the parameters
param vext=8, iapp=0, freq=0.3
param gna=120
param gk=36, gl=0.3, gtrp=0.03, cm=1
number vna=50, vk=-80, vtrp=0, vl=-49
# time constants and steady state functions for gating variables
am(v)=0.1*(v+40)/(1-exp(-(v+40)/10))
bm(v)=4*exp(-(v+65)/18)
ah(v)=0.07*exp(-(v+65)/20)
bh(v)=1/( 1+ exp(-(v+35)/10) )
an(v)=0.01*(v+55)/(1-exp(-(V + 55)/10))
bn(v)=0.125*exp(-(V+65)/80)
vex(t)=vext*(sin(2*pi*(freq/1000)*t))
ik=gk*n^4*(v+vex(t)-vk)
ina=gna*m^3*h*(v+vex(t)-vna)
il=gl*(v+vex(t)-vl)
itrp=gtrp*(v+vex(t)-vtrp)
# the equations
v'=(-ik-ina-il-itrp+iapp)/cm
m'=(am(v)*(1-m)-bm(v)*m)
n'=(an(v)*(1-n)-bn(v)*n)
h'=(ah(v)*(1-h)-bh(v)*h)
aux vex=vex(t)
aux itrp=itrp
# set xpp parameters
@ total=10000, xp=t, yp=v, xlo=0, xhi=10000, ylo=-90, yhi=50, bounds=10000000, dt=0.01
@ maxstor=10000000
done
Terms of use
Copyright © 2017 Matthias Koenig
Redistribution and use of any part of this model, with or without modification, are permitted provided that
the following conditions are met:
Redistributions of this SBML file must retain the above copyright notice, this list of conditions
and the following disclaimer.
Redistributions in a different form must reproduce the above copyright notice, this list of
conditions and the following disclaimer in the documentation and/or other materials provided
with the distribution.
This model is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
Model :
id
name
time
substance
extent
volume
area
length
L3V1
type
FunctionDefinitions [11]
name
math
sbo
cvterm
FunctionDefinition
max
minimum
x
y
x
x
y
y
FunctionDefinition
min
maximum
x
y
x
x
y
y
FunctionDefinition
heav
heavyside
x
0
x
0
0.5
x
0
1
x
0
0
FunctionDefinition
mod
modulo
x
y
x
y
x
y
x
0
y
0
x
y
x
y
FunctionDefinition
am
v
0.1
v
40
1
v
40
10
FunctionDefinition
bm
v
4
v
65
18
FunctionDefinition
ah
v
0.07
v
65
20
FunctionDefinition
bh
v
1
1
v
35
10
FunctionDefinition
an
v
0.01
v
55
1
v
55
10
FunctionDefinition
bn
v
0.125
v
65
80
FunctionDefinition
vex
t
freq
vext
vext
2
freq
1000
t
type
Parameters [22]
name
constant
value
unit
derived unit
sbo
cvterm
Parameter
v
v = -71
F
-71.0
None
Parameter
m
m = 0.038
F
0.038
None
Parameter
h
h = 0.69
F
0.69
None
Parameter
n
n = 0.159
F
0.159
None
Parameter
vext
vext = 8
F
8.0
None
Parameter
iapp
iapp = 0
F
0.0
None
Parameter
freq
freq = 0.3
F
0.3
None
Parameter
gna
gna = 120
F
120.0
None
Parameter
gk
gk = 36
F
36.0
None
Parameter
gl
gl = 0.3
F
0.3
None
Parameter
gtrp
gtrp = 0.03
F
0.03
None
Parameter
cm
cm = 1
F
1.0
None
Parameter
vna
vna = 50
F
50.0
None
Parameter
vk
vk = -80
F
-80.0
None
Parameter
vtrp
vtrp = 0
F
0.0
None
Parameter
vl
vl = -49
F
-49.0
None
Parameter
ik
F
0.0
dimensionless
None
Parameter
ina
F
0.0
dimensionless
None
Parameter
il
F
0.0
dimensionless
None
Parameter
itrp
F
0.0
dimensionless
None
Parameter
vex
F
0.0
dimensionless
None
Parameter
t
model time
F
0.0
dimensionless
None
type
Rules [10]
assignment
name
derived units
sbo
cvterm
Rule
d v/dt
=
ik
ina
il
itrp
iapp
cm
None
Rule
d m/dt
=
am
v
1
m
bm
v
m
None
Rule
d n/dt
=
an
v
1
n
bn
v
n
None
Rule
d h/dt
=
ah
v
1
h
bh
v
h
None
Rule
ik
=
gk
n
4
v
vex
t
freq
vext
vk
None
Rule
ina
=
gna
m
3
h
v
vex
t
freq
vext
vna
None
Rule
il
=
gl
v
vex
t
freq
vext
vl
None
Rule
itrp
=
gtrp
v
vex
t
freq
vext
vtrp
None
Rule
vex
=
vex
t
freq
vext
None
Rule
t
=
time
None