XPP model

This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6.

# numbers:
number Cm=1.0

# parameters:
p Iapp=3.35
p gA=0.39
p theta_m=-24.0
p gNa=112.5,gK=225.0,gL=0.25
p sigma_m=11.5
p theta_h=-58.3  , sigma_h=-6.7
p theta_n=-12.4  , sigma_n=6.8
p theta_t_h=-60  , sigma_t_h=-12.0
p theta_tna=-14.6, sigma_tna=-8.6
p theta_tnb=1.3  , sigma_tnb=18.7
p theta_a=-50    , sigma_a=20
p theta_b=-70    , sigma_b=-6
p tau_b=150, tau_a=2
p power_n=2.0
p V_Na=50.0,V_K=-90.0,V_L=-70.0

# auxilary functions:
GAMMAF(VV,theta,sigma)=1.0/(1.0+exp(-(VV-theta)/sigma))

# innitial conditions
V(0)=-70.038
h(0)=0.8522
n(0)=0.000208
a(0)=0.2686
b(0)=0.5016
# functions:
m_inf=GAMMAF(V,theta_m,sigma_m)
h_inf=GAMMAF(V,theta_h,sigma_h)
n_inf=GAMMAF(V,theta_n,sigma_n)
a_inf=GAMMAF(V,theta_a,sigma_a)
b_inf=GAMMAF(V,theta_b,sigma_b)
tau_h(V)=0.5+14.0*GAMMAF(V,theta_t_h,sigma_t_h)
tau_n(V)=(0.087+11.4*GAMMAF(V,theta_tna,sigma_tna))*(0.087+11.4*GAMMAF(V,theta_tnb,sigma_tnb))

# ode's
V'=-gNa*m_inf^3*h*(V-V_Na)-gK*(n^power_n)*(V-V_K)-gL*(V-V_L)-gA*a^3*b*(V-V_k)+Iapp
h'=(h_inf-h)/tau_h (V)
n'=(n_inf-n)/tau_n(V)
a'=(a_inf-a)/tau_a
b'=(b_inf-b)/tau_b

# 
@ MAXSTOR=16000000
@ BACK=Black
@ XP=T
@ YP=Vs
@ AXES=2
@ TOTAL=1000.0
@ DT=0.01
@ NJMP=1
@ T0=0.0
@ TRANS=0.0
@ NMESH=40
@ METH=rungekutta
@ DTMIN=0.001
@ DTMAX=1.0
@ TOLER=0.00001
@ BOUND=10000.0
@ DELAY=0.0
@ XLO=0.0, XHI=1000.0, YLO=-90.0, YHI=50.0
@ NTST=50,NMAX=20000,NPR=50
@ DS=0.02,DSMIN=0.001,DSMAX=0.5
@ PARMIN=-10,PARMAX=200,NORMMIN=0.0,NORMMAX=10000.0
@ AUTOVAR=V,AUTOXMIN=-10.0,AUTOXMAX=100.0,AUTOYMIN=-90.0,AUTOYMAX=90.0

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 [7] 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
gammaf vv theta sigma 1 1 vv theta sigma
tau_h v sigma_t_h theta_t_h 0.5 14 gammaf v theta_t_h sigma_t_h
tau_n v sigma_tna sigma_tnb theta_tna theta_tnb 0.087 11.4 gammaf v theta_tna sigma_tna 0.087 11.4 gammaf v theta_tnb sigma_tnb

Parameters [39] name constant value unit derived unit sbo cvterm
cm cm = 1.0 1.0 None
iapp iapp = 3.35 3.35 None
ga ga = 0.39 0.39 None
theta_m theta_m = -24.0 -24.0 None
gna gna = 112.5 112.5 None
gk gk = 225.0 225.0 None
gl gl = 0.25 0.25 None
sigma_m sigma_m = 11.5 11.5 None
theta_h theta_h = -58.3 -58.3 None
sigma_h sigma_h = -6.7 -6.7 None
theta_n theta_n = -12.4 -12.4 None
sigma_n sigma_n = 6.8 6.8 None
theta_t_h theta_t_h = -60 -60.0 None
sigma_t_h sigma_t_h = -12.0 -12.0 None
theta_tna theta_tna = -14.6 -14.6 None
sigma_tna sigma_tna = -8.6 -8.6 None
theta_tnb theta_tnb = 1.3 1.3 None
sigma_tnb sigma_tnb = 18.7 18.7 None
theta_a theta_a = -50 -50.0 None
sigma_a sigma_a = 20 20.0 None
theta_b theta_b = -70 -70.0 None
sigma_b sigma_b = -6 -6.0 None
tau_b tau_b = 150 150.0 None
tau_a tau_a = 2 2.0 None
power_n power_n = 2.0 2.0 None
v_na v_na = 50.0 50.0 None
v_k v_k = -90.0 -90.0 None
v_l v_l = -70.0 -70.0 None
v v = -70.038 -70.038 None
h h = 0.8522 0.8522 None
n n = 0.000208 0.000208 None
a a = 0.2686 0.2686 None
b b = 0.5016 0.5016 None
m_inf 0.0 dimensionless None
h_inf 0.0 dimensionless None
n_inf 0.0 dimensionless None
a_inf 0.0 dimensionless None
b_inf 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [11]   assignment name derived units sbo cvterm
d v/dt = gna m_inf 3 h v v_na gk n power_n v v_k gl v v_l ga a 3 b v v_k iapp None
d h/dt = h_inf h tau_h v None
d n/dt = n_inf n tau_n v sigma_tna sigma_tnb theta_tna theta_tnb None
d a/dt = a_inf a tau_a None
d b/dt = b_inf b tau_b None
m_inf = gammaf v theta_m sigma_m None
h_inf = gammaf v theta_h sigma_h None
n_inf = gammaf v theta_n sigma_n None
a_inf = gammaf v theta_a sigma_a None
b_inf = gammaf v theta_b sigma_b None
t = time None