XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# from Rho and Prescott, PLoS Comput Biol 2012
# to be run in XPP
# standard 2-D model (see Fig 2)
# gna corresponds to gfast; gk corresponds to gslow; gl corresponds to gleak
# DIFFERENTIAL EQUATIONS
dv/dt = (Istim-gna*minf(V)*(V-Vna)-gk*w*(V-VK)-gl*(V-Vl))/cap
dw/dt = phi_w*(winf(V)-w)/tauw(V)
# FUNCTIONS AND PARAMETERS
minf(v)=.5*(1+tanh((v-beta_m)/gamma_m))
winf(v)=.5*(1+tanh((v-beta_w)/gamma_w))
tauw(v)=1/cosh((v-beta_w)/(2*gamma_w))
param Istim=0
param vna=50,vk=-100,vl=-70
param gk=20,gl=2,gna=20
param beta_m=-1.2,gamma_m=18
param beta_w=-13,gamma_w=10
# beta_w = -21 for transient spiking; beta_w = -13 for repetitive spiking
param phi_w=.15
param cap=2
# INITIAL CONDITIONS
V(0)=-70
w(0)=0.000025
# ALWAYS USE EULER! - Actually this is only true for noise
@ total=10000,dt=.05,xlo=-100,xhi=60,ylo=-.125,yhi=.6,xp=v,yp=w
@ meth=euler
@ MAXSTOR=1000000,bounds=10000
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 [7]
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
minf
v
beta_m
gamma_m
0.5
1
v
beta_m
gamma_m
FunctionDefinition
winf
v
beta_w
gamma_w
0.5
1
v
beta_w
gamma_w
FunctionDefinition
tauw
v
beta_w
gamma_w
1
v
beta_w
2
gamma_w
type
Parameters [16]
name
constant
value
unit
derived unit
sbo
cvterm
Parameter
istim
istim = 0
F
0.0
None
Parameter
vna
vna = 50
F
50.0
None
Parameter
vk
vk = -100
F
-100.0
None
Parameter
vl
vl = -70
F
-70.0
None
Parameter
gk
gk = 20
F
20.0
None
Parameter
gl
gl = 2
F
2.0
None
Parameter
gna
gna = 20
F
20.0
None
Parameter
beta_m
beta_m = -1.2
F
-1.2
None
Parameter
gamma_m
gamma_m = 18
F
18.0
None
Parameter
beta_w
beta_w = -13
F
-13.0
None
Parameter
gamma_w
gamma_w = 10
F
10.0
None
Parameter
phi_w
phi_w = .15
F
0.15
None
Parameter
cap
cap = 2
F
2.0
None
Parameter
v
v = -70
F
-70.0
None
Parameter
w
w = 0.000025
F
2.5e-05
None
Parameter
t
model time
F
0.0
dimensionless
None
type
Rules [3]
assignment
name
derived units
sbo
cvterm
Rule
d v/dt
=
istim
gna
minf
v
beta_m
gamma_m
v
vna
gk
w
v
vk
gl
v
vl
cap
None
Rule
d w/dt
=
phi_w
winf
v
beta_w
gamma_w
w
tauw
v
beta_w
gamma_w
None
Rule
t
=
time
None