XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
% Model for action potential of skeletal muscle fiber
% Adrian RH, Chandler WK, Hodgkin AL. Voltage clamp experiments in striated muscle
% fibres. J Physiol 208:607-44, 1970.
% implemented by Dr. Sheng-Nan Wu
% Initial values
Initial Vm=-95, Vt=-95, m=0.0, h=1.0, n=0.0
% Stimulus
Param period=200, iStim_mag=3, iStim_beg=1, iStim_dur=0.5
iStim= iStim_mag * heav(mod(t,period)-iStim_beg) * heav(iStim_beg+iStim_dur-mod(t,period))
% Values of the model parameters
Param gNa_max=1.0, gK_max=0.415, gL_max=0.0024
Param ENa=50.0, EK=-70.0, EL=-95.0, Cm=0.0090, Ct=0.04, Rs=15.0
Param alpha_m_max=0.208
Param beta_m_max=2.081
Param alpha_n_max=0.0229
Param beta_n_max=0.09616
Param alpha_h_max=0.0156
Param beta_h_max=3.382
Param Em=-42.0, En=-40.0, Eh=-41.0
Param v_alpha_m=10.0, v_alpha_h=14.7, v_alpha_n=7.0
Param v_beta_m=18.0, v_beta_h=7.6, v_beta_n=40.0
% Expressions
Ina= (gNa_max * m * m * m * h * (Vm - ENa))
beta_n= (beta_n_max * exp(((En - Vm) / v_beta_n)))
beta_m= (beta_m_max * exp(((Em - Vm) / v_beta_m)))
beta_h= (beta_h_max / (1.0 + exp(((Eh - Vm) / v_beta_h))))
IT= ((Vm - Vt) / Rs)
IL= (gL_max * (Vm - EL))
IK= (gK_max * n * n * n * n * (Vm - EK))
alpha_n= (alpha_n_max * (Vm - En) / (1.0 - exp(((En - Vm) / v_alpha_n))))
alpha_m= (alpha_m_max * (Vm - Em) / (1.0 - exp(((Em - Vm) / v_alpha_m))))
alpha_h= (alpha_h_max * exp(((Eh - Vm) / v_alpha_h)))
% Differential equations
dVm/dt= ((iStim - (INa + IK + IL + IT)) / Cm)
dm/dt= ((alpha_m * (1.0 - m)) - (beta_m * m))
dh/dt= ((alpha_h * (1.0 - h)) - (beta_h * h))
dn/dt= ((alpha_n * (1.0 - n)) - (beta_n * n))
dVt/dt=((Vm - Vt) / (Rs * Ct))
% Numerical and plotting parameters for xpp
@ meth=Euler, dt=0.01, total=25
@ yp=v, yhi=50, ylo=-110, xlo=0, xhi=25, bounds=5000
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 [4]
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
type
Parameters [45]
name
constant
value
unit
derived unit
sbo
cvterm
Parameter
vm
vm = -95
F
-95.0
None
Parameter
vt
vt = -95
F
-95.0
None
Parameter
m
m = 0.0
F
0.0
None
Parameter
h
h = 1.0
F
1.0
None
Parameter
n
n = 0.0
F
0.0
None
Parameter
period
period = 200
F
200.0
None
Parameter
istim_mag
istim_mag = 3
F
3.0
None
Parameter
istim_beg
istim_beg = 1
F
1.0
None
Parameter
istim_dur
istim_dur = 0.5
F
0.5
None
Parameter
gna_max
gna_max = 1.0
F
1.0
None
Parameter
gk_max
gk_max = 0.415
F
0.415
None
Parameter
gl_max
gl_max = 0.0024
F
0.0024
None
Parameter
ena
ena = 50.0
F
50.0
None
Parameter
ek
ek = -70.0
F
-70.0
None
Parameter
el
el = -95.0
F
-95.0
None
Parameter
cm
cm = 0.0090
F
0.009
None
Parameter
ct
ct = 0.04
F
0.04
None
Parameter
rs
rs = 15.0
F
15.0
None
Parameter
alpha_m_max
alpha_m_max = 0.208
F
0.208
None
Parameter
beta_m_max
beta_m_max = 2.081
F
2.081
None
Parameter
alpha_n_max
alpha_n_max = 0.0229
F
0.0229
None
Parameter
beta_n_max
beta_n_max = 0.09616
F
0.09616
None
Parameter
alpha_h_max
alpha_h_max = 0.0156
F
0.0156
None
Parameter
beta_h_max
beta_h_max = 3.382
F
3.382
None
Parameter
em
em = -42.0
F
-42.0
None
Parameter
en
en = -40.0
F
-40.0
None
Parameter
eh
eh = -41.0
F
-41.0
None
Parameter
v_alpha_m
v_alpha_m = 10.0
F
10.0
None
Parameter
v_alpha_h
v_alpha_h = 14.7
F
14.7
None
Parameter
v_alpha_n
v_alpha_n = 7.0
F
7.0
None
Parameter
v_beta_m
v_beta_m = 18.0
F
18.0
None
Parameter
v_beta_h
v_beta_h = 7.6
F
7.6
None
Parameter
v_beta_n
v_beta_n = 40.0
F
40.0
None
Parameter
istim
F
0.0
dimensionless
None
Parameter
ina
F
0.0
dimensionless
None
Parameter
beta_n
F
0.0
dimensionless
None
Parameter
beta_m
F
0.0
dimensionless
None
Parameter
beta_h
F
0.0
dimensionless
None
Parameter
it
F
0.0
dimensionless
None
Parameter
il
F
0.0
dimensionless
None
Parameter
ik
F
0.0
dimensionless
None
Parameter
alpha_n
F
0.0
dimensionless
None
Parameter
alpha_m
F
0.0
dimensionless
None
Parameter
alpha_h
F
0.0
dimensionless
None
Parameter
t
model time
F
0.0
dimensionless
None
type
Rules [17]
assignment
name
derived units
sbo
cvterm
Rule
d vm/dt
=
istim
ina
ik
il
it
cm
None
Rule
d m/dt
=
alpha_m
1
m
beta_m
m
None
Rule
d h/dt
=
alpha_h
1
h
beta_h
h
None
Rule
d n/dt
=
alpha_n
1
n
beta_n
n
None
Rule
d vt/dt
=
vm
vt
rs
ct
None
Rule
istim
=
istim_mag
heav
mod
t
period
istim_beg
heav
istim_beg
istim_dur
mod
t
period
None
Rule
ina
=
gna_max
m
m
m
h
vm
ena
None
Rule
beta_n
=
beta_n_max
en
vm
v_beta_n
None
Rule
beta_m
=
beta_m_max
em
vm
v_beta_m
None
Rule
beta_h
=
beta_h_max
1
eh
vm
v_beta_h
None
Rule
it
=
vm
vt
rs
None
Rule
il
=
gl_max
vm
el
None
Rule
ik
=
gk_max
n
n
n
n
vm
ek
None
Rule
alpha_n
=
alpha_n_max
vm
en
1
en
vm
v_alpha_n
None
Rule
alpha_m
=
alpha_m_max
vm
em
1
em
vm
v_alpha_m
None
Rule
alpha_h
=
alpha_h_max
eh
vm
v_alpha_h
None
Rule
t
=
time
None