XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# Ermentrout 1998 type I oscillator p gl=.2,gkdr=80,gna=80 p iapp=0.67 p C=1. p amphi=.32,amhalf=-54,amwidth=4 p bmphi=.28,bmhalf=-27,bmwidth=5 p ahphi=.128,ahhalf=-50,ahwidth=18 p bhphi=4,bhhalf=-27,bhwidth=5 p anphi=.032,anhalf=-52,anwidth=5 p bnphi=.5,bnhalf=-57,bnwidth=40 # functions am(v)=amphi*(v-amhalf)/(1-exp(-(v-amhalf)/amwidth)) bm(v)=bmphi*(v-bmhalf)/(exp((v-bmhalf)/bmwidth)-1) ah(v)=ahphi*exp(-(v-ahhalf)/ahwidth) bh(v)=bhphi/(1+exp(-(v-bhhalf)/bhwidth)) an(v)=anphi*(v-anhalf)/(1-exp(-(v-anhalf)/anwidth)) bn(v)=bnphi*exp(-(v-bnhalf)/bnwidth) Isyn(y,V)=gsyn*y*(V-Esyn) #currents ina(v,m,h)=gna*m^3*h*(v-Vna) ikdr(v,n)=gkdr*n^4*(v-Vk) il(v)=gl*(v-Vl) #diff. equ. v1'=(iapp-(ina(v1,m1,ha)+ikdr(v1,n1)+il(v1))+eps*(ua1-v1)/dx-p0*Isyn(y2,v1))/C v2'=(iapp-(ina(v2,m2,hb)+ikdr(v2,n2)+il(v2))+eps*(ub1-v2)/dx-p0*Isyn(y1,v2))/C m1'=am(v1)*(1-m1)-bm(v1)*m1 m2'=am(v2)*(1-m2)-bm(v2)*m2 n1'=an(v1)*(1-n1)-bn(v1)*n1 n2'=an(v2)*(1-n2)-bn(v2)*n2 ha'=ah(v1)*(1-ha)-bh(v1)*ha hb'=ah(v2)*(1-hb)-bh(v2)*hb ## synapse par taur=1,taud=3,thresh=-30, gsyn=1 Esyn=0 x1'=(-x1+.5*(1+tanh((v1-thresh)/3.0)))/taur x2'=(-x2+.5*(1+tanh((v2-thresh)/3.0)))/taur y1'=(-y1+x1)/taud y2'=(-y2+x2)/taud init x1=0,y1=0,x2=0,y2=0 par Vna=50,Vk=-100,Vl=-67,Vh=-65,Vrh=-40,Vsh=6,gld=0.1,gh=0.02,tauh=400 hinf(V)=1/(1+exp((V-Vh)/Vsh)) Ild(V)=V-Vl Ih(V,y)=gh*y*(V-Vrh)/gld ha[1..50]'=(hinf(ua[j])-ha[j])/tauh hb[1..50]'=(hinf(ub[j])-hb[j])/tauh # cable equation p p[0..50]=0 ua1'=((lambda/dx)^2*(ua2-2*ua1+v1)-Ild(ua1)-Ih(ua1,ha1)-p1*Isyn(y2,ua1)/gld)/tau ua[2..50]'= ((lambda/dx)^2*(ua[j+1]-2*ua[j]+ua[j-1])-Ild(ua[j])-Ih(ua[j],ha[j])-p[j]*Isyn(y2,ua[j])/gld)/tau ua51=(c1+b1*ua50/dx)/(a1+b1/dx) ub1'=((lambda/dx)^2*(ub2-2*ub1+v2)-Ild(ub1)-Ih(ub1,hb1)-p1*Isyn(y1,ub1)/gld)/tau ub[2..50]'= ((lambda/dx)^2*(ub[j+1]-2*ub[j]+ub[j-1])-Ild(ub[j])-Ih(ub[j],hb[j])-p[j]*Isyn(y1,ub[j])/gld)/tau ub51=(c1+b1*ub50/dx)/(a1+b1/dx) par lambda=1,tau=10,dx=.1,c1=0,a1=0,b1=1,c0=0,a0=0,b0=1,eps=.025 @ total=1000,dt=.05,xlo=0,xhi=1000,ylo=-100,yhi=50,bounds=1e6 i ua[1..50]=-60 i ha[1..50]=.1 i ub[1..50]=-60 i hb[1..50]=.1 d
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:
- 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.
Model :
|
![]() |
type | FunctionDefinitions [17] | name | math | sbo | cvterm |
---|---|---|---|---|---|
FunctionDefinition | max
|
minimum | |||
FunctionDefinition | min
|
maximum | |||
FunctionDefinition | heav
|
heavyside | |||
FunctionDefinition | mod
|
modulo | |||
FunctionDefinition | am
|
||||
FunctionDefinition | bm
|
||||
FunctionDefinition | ah
|
||||
FunctionDefinition | bh
|
||||
FunctionDefinition | an
|
||||
FunctionDefinition | bn
|
||||
FunctionDefinition | isyn
|
||||
FunctionDefinition | ina
|
||||
FunctionDefinition | ikdr
|
||||
FunctionDefinition | il
|
||||
FunctionDefinition | hinf
|
||||
FunctionDefinition | ild
|
||||
FunctionDefinition | ih
|
type | Parameters [73] | name | constant | value | unit | derived unit | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Parameter | gl
|
gl = .2 | F | 0.2 | None | |||
Parameter | gkdr
|
gkdr = 80 | F | 80.0 | None | |||
Parameter | gna
|
gna = 80 | F | 80.0 | None | |||
Parameter | iapp
|
iapp = 0.67 | F | 0.67 | None | |||
Parameter | c
|
c = 1. | F | 1.0 | None | |||
Parameter | amphi
|
amphi = .32 | F | 0.32 | None | |||
Parameter | amhalf
|
amhalf = -54 | F | -54.0 | None | |||
Parameter | amwidth
|
amwidth = 4 | F | 4.0 | None | |||
Parameter | bmphi
|
bmphi = .28 | F | 0.28 | None | |||
Parameter | bmhalf
|
bmhalf = -27 | F | -27.0 | None | |||
Parameter | bmwidth
|
bmwidth = 5 | F | 5.0 | None | |||
Parameter | ahphi
|
ahphi = .128 | F | 0.128 | None | |||
Parameter | ahhalf
|
ahhalf = -50 | F | -50.0 | None | |||
Parameter | ahwidth
|
ahwidth = 18 | F | 18.0 | None | |||
Parameter | bhphi
|
bhphi = 4 | F | 4.0 | None | |||
Parameter | bhhalf
|
bhhalf = -27 | F | -27.0 | None | |||
Parameter | bhwidth
|
bhwidth = 5 | F | 5.0 | None | |||
Parameter | anphi
|
anphi = .032 | F | 0.032 | None | |||
Parameter | anhalf
|
anhalf = -52 | F | -52.0 | None | |||
Parameter | anwidth
|
anwidth = 5 | F | 5.0 | None | |||
Parameter | bnphi
|
bnphi = .5 | F | 0.5 | None | |||
Parameter | bnhalf
|
bnhalf = -57 | F | -57.0 | None | |||
Parameter | bnwidth
|
bnwidth = 40 | F | 40.0 | None | |||
Parameter | taur
|
taur = 1 | F | 1.0 | None | |||
Parameter | taud
|
taud = 3 | F | 3.0 | None | |||
Parameter | thresh
|
thresh = -30 | F | -30.0 | None | |||
Parameter | gsyn
|
gsyn = 1 | F | 1.0 | None | |||
Parameter | esyn
|
esyn = 0 | F | 0.0 | None | |||
Parameter | x1
|
x1 = 0 | F | 0.0 | None | |||
Parameter | y1
|
y1 = 0 | F | 0.0 | None | |||
Parameter | x2
|
x2 = 0 | F | 0.0 | None | |||
Parameter | y2
|
y2 = 0 | F | 0.0 | None | |||
Parameter | vna
|
vna = 50 | F | 50.0 | None | |||
Parameter | vk
|
vk = -100 | F | -100.0 | None | |||
Parameter | vl
|
vl = -67 | F | -67.0 | None | |||
Parameter | vh
|
vh = -65 | F | -65.0 | None | |||
Parameter | vrh
|
vrh = -40 | F | -40.0 | None | |||
Parameter | vsh
|
vsh = 6 | F | 6.0 | None | |||
Parameter | gld
|
gld = 0.1 | F | 0.1 | None | |||
Parameter | gh
|
gh = 0.02 | F | 0.02 | None | |||
Parameter | tauh
|
tauh = 400 | F | 400.0 | None | |||
Parameter | p[0..50] = 0 | F | 0.0 | None | ||||
Parameter | lambda
|
lambda = 1 | F | 1.0 | None | |||
Parameter | tau
|
tau = 10 | F | 10.0 | None | |||
Parameter | dx
|
dx = .1 | F | 0.1 | None | |||
Parameter | c1
|
c1 = 0 | F | 0.0 | None | |||
Parameter | a1
|
a1 = 0 | F | 0.0 | None | |||
Parameter | b1
|
b1 = 1 | F | 1.0 | None | |||
Parameter | c0
|
c0 = 0 | F | 0.0 | None | |||
Parameter | a0
|
a0 = 0 | F | 0.0 | None | |||
Parameter | b0
|
b0 = 1 | F | 1.0 | None | |||
Parameter | eps
|
eps = .025 | F | 0.025 | None | |||
Parameter | ua[1..50] = -60 | F | -60.0 | None | ||||
Parameter | ha[1..50] = .1 | F | 0.1 | None | ||||
Parameter | ub[1..50] = -60 | F | -60.0 | None | ||||
Parameter | hb[1..50] = .1 | F | 0.1 | None | ||||
Parameter | v1
|
F | 0.0 | dimensionless | None | |||
Parameter | v2
|
F | 0.0 | dimensionless | None | |||
Parameter | m1
|
F | 0.0 | dimensionless | None | |||
Parameter | m2
|
F | 0.0 | dimensionless | None | |||
Parameter | n1
|
F | 0.0 | dimensionless | None | |||
Parameter | n2
|
F | 0.0 | dimensionless | None | |||
Parameter | ha
|
F | 0.0 | dimensionless | None | |||
Parameter | hb
|
F | 0.0 | dimensionless | None | |||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | ua1
|
F | 0.0 | dimensionless | None | |||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | ub1
|
F | 0.0 | dimensionless | None | |||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | ua51
|
F | 0.0 | dimensionless | None | |||
Parameter | ub51
|
F | 0.0 | dimensionless | None | |||
Parameter | t
|
model time | F | 0.0 | dimensionless | None |
type | Rules [21] | assignment | name | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
Rule | d v1/dt
|
= | None | ||||
Rule | d v2/dt
|
= | None | ||||
Rule | d m1/dt
|
= | None | ||||
Rule | d m2/dt
|
= | None | ||||
Rule | d n1/dt
|
= | None | ||||
Rule | d n2/dt
|
= | None | ||||
Rule | d ha/dt
|
= | None | ||||
Rule | d hb/dt
|
= | None | ||||
Rule | d x1/dt
|
= | None | ||||
Rule | d x2/dt
|
= | None | ||||
Rule | d y1/dt
|
= | None | ||||
Rule | d y2/dt
|
= | None | ||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | d ua1/dt
|
= | None | ||||
Rule | = | None | |||||
Rule | d ub1/dt
|
= | None | ||||
Rule | = | None | |||||
Rule | ua51
|
= | None | ||||
Rule | ub51
|
= | None | ||||
Rule | t
|
= | None |