XPP model

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

# Simplified 2D HH model
# formatted for use with XPPAUT

# Parameters as in the model of Figure 8, Lundstrom et al, J Comp Neurosci
# Note: the z variable introduces stimulus noise and must be removed for bifurcation analysis

wiener nz
v' = -(Gl*(v - El) + Gk*n^4*(v - Ek) + Gna*minf^3*h*(v - Ena))/cm + i/cm + z
n' = (ninf - n) / ntau  
z'=-z/tau + nsig*sqrt(2/tau)*nz
par tau=1

par i=0, mVh=-40,mk=7,nVh=-45,nk=15, Gna=50,Gk=36,Gl=5,Ena=50,Ek=-77,El=-54,cm=1
par nsig=0

par Cbase=3,Camp=50,Vmax=-50,sig=1

ntau = Cbase + Camp*exp(-(Vmax - V)^2/sig^2)
ninf = 1.0 / (1+exp( (nVh - v)/nk) )
minf = 1.0 / (1+exp( (mVh - v)/mk) )
h = .89 - 1.1*n

init v=-65,n=0.3

@ meth=euler,total=20ee0,dt=.005,nout=50
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 [27] name constant value unit derived unit sbo cvterm
nz 0.0 None
tau tau = 1 1.0 None
i i = 0 0.0 None
mvh mvh = -40 -40.0 None
mk mk = 7 7.0 None
nvh nvh = -45 -45.0 None
nk nk = 15 15.0 None
gna gna = 50 50.0 None
gk gk = 36 36.0 None
gl gl = 5 5.0 None
ena ena = 50 50.0 None
ek ek = -77 -77.0 None
el el = -54 -54.0 None
cm cm = 1 1.0 None
nsig nsig = 0 0.0 None
cbase cbase = 3 3.0 None
camp camp = 50 50.0 None
vmax vmax = -50 -50.0 None
sig sig = 1 1.0 None
v v = -65 -65.0 None
n n = 0.3 0.3 None
z 0.0 dimensionless None
ntau 0.0 dimensionless None
ninf 0.0 dimensionless None
minf 0.0 dimensionless None
h 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [8]   assignment name derived units sbo cvterm
d v/dt = gl v el gk n 4 v ek gna minf 3 h v ena cm i cm z None
d n/dt = ninf n ntau None
d z/dt = z tau nsig 2 2 tau nz None
ntau = cbase camp vmax v 2 sig 2 None
ninf = 1 1 nvh v nk None
minf = 1 1 mvh v mk None
h = 0.89 1.1 n None
t = time None