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
This file has been produced by sbmlutils.

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:

  1. Redistributions of this SBML file must retain the above copyright notice, this list of conditions and the following disclaimer.
  2. 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
Access SBML model  L3V1

FunctionDefinitions [4] name math sbo cvterm
max minimum x y x x y y
min maximum x y x x y y
heav heavyside x 0 x 0 0.5 x 0 1 x 0 0
mod modulo x y x y x y x 0 y 0 x y x y

Parameters [45] name constant value unit derived unit sbo cvterm
vm vm = -95 -95.0 None
vt vt = -95 -95.0 None
m m = 0.0 0.0 None
h h = 1.0 1.0 None
n n = 0.0 0.0 None
period period = 200 200.0 None
istim_mag istim_mag = 3 3.0 None
istim_beg istim_beg = 1 1.0 None
istim_dur istim_dur = 0.5 0.5 None
gna_max gna_max = 1.0 1.0 None
gk_max gk_max = 0.415 0.415 None
gl_max gl_max = 0.0024 0.0024 None
ena ena = 50.0 50.0 None
ek ek = -70.0 -70.0 None
el el = -95.0 -95.0 None
cm cm = 0.0090 0.009 None
ct ct = 0.04 0.04 None
rs rs = 15.0 15.0 None
alpha_m_max alpha_m_max = 0.208 0.208 None
beta_m_max beta_m_max = 2.081 2.081 None
alpha_n_max alpha_n_max = 0.0229 0.0229 None
beta_n_max beta_n_max = 0.09616 0.09616 None
alpha_h_max alpha_h_max = 0.0156 0.0156 None
beta_h_max beta_h_max = 3.382 3.382 None
em em = -42.0 -42.0 None
en en = -40.0 -40.0 None
eh eh = -41.0 -41.0 None
v_alpha_m v_alpha_m = 10.0 10.0 None
v_alpha_h v_alpha_h = 14.7 14.7 None
v_alpha_n v_alpha_n = 7.0 7.0 None
v_beta_m v_beta_m = 18.0 18.0 None
v_beta_h v_beta_h = 7.6 7.6 None
v_beta_n v_beta_n = 40.0 40.0 None
istim 0.0 dimensionless None
ina 0.0 dimensionless None
beta_n 0.0 dimensionless None
beta_m 0.0 dimensionless None
beta_h 0.0 dimensionless None
it 0.0 dimensionless None
il 0.0 dimensionless None
ik 0.0 dimensionless None
alpha_n 0.0 dimensionless None
alpha_m 0.0 dimensionless None
alpha_h 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [17]   assignment name derived units sbo cvterm
d vm/dt = istim ina ik il it cm None
d m/dt = alpha_m 1 m beta_m m None
d h/dt = alpha_h 1 h beta_h h None
d n/dt = alpha_n 1 n beta_n n None
d vt/dt = vm vt rs ct None
istim = istim_mag heav mod t period istim_beg heav istim_beg istim_dur mod t period None
ina = gna_max m m m h vm ena None
beta_n = beta_n_max en vm v_beta_n None
beta_m = beta_m_max em vm v_beta_m None
beta_h = beta_h_max 1 eh vm v_beta_h None
it = vm vt rs None
il = gl_max vm el None
ik = gk_max n n n n vm ek None
alpha_n = alpha_n_max vm en 1 en vm v_alpha_n None
alpha_m = alpha_m_max vm em 1 em vm v_alpha_m None
alpha_h = alpha_h_max eh vm v_alpha_h None
t = time None