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
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 [11] 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
am v 0.1 v 40 1 v 40 10
bm v 4 v 65 18
ah v 0.07 v 65 20
bh v 1 1 v 35 10
an v 0.01 v 55 1 v 55 10
bn v 0.125 v 65 80
vex t freq vext vext 2 freq 1000 t

Parameters [22] name constant value unit derived unit sbo cvterm
v v = -71 -71.0 None
m m = 0.038 0.038 None
h h = 0.69 0.69 None
n n = 0.159 0.159 None
vext vext = 8 8.0 None
iapp iapp = 0 0.0 None
freq freq = 0.3 0.3 None
gna gna = 120 120.0 None
gk gk = 36 36.0 None
gl gl = 0.3 0.3 None
gtrp gtrp = 0.03 0.03 None
cm cm = 1 1.0 None
vna vna = 50 50.0 None
vk vk = -80 -80.0 None
vtrp vtrp = 0 0.0 None
vl vl = -49 -49.0 None
ik 0.0 dimensionless None
ina 0.0 dimensionless None
il 0.0 dimensionless None
itrp 0.0 dimensionless None
vex 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [10]   assignment name derived units sbo cvterm
d v/dt = ik ina il itrp iapp cm None
d m/dt = am v 1 m bm v m None
d n/dt = an v 1 n bn v n None
d h/dt = ah v 1 h bh v h None
ik = gk n 4 v vex t freq vext vk None
ina = gna m 3 h v vex t freq vext vna None
il = gl v vex t freq vext vl None
itrp = gtrp v vex t freq vext vtrp None
vex = vex t freq vext None
t = time None