XPP model

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

#this is the stimulation protocol with three electrodes attached to 5th, 6th and 7th STN units 
#with a local (w1=0.3,w2=0) or non-local (w1=0.3,w2=0.1) spread of stimulation signals.
#Used to generate figures 5 and 7 in PLoS One paper.

#stimulation parameters
par w1=0.3,w2=0,Kn=60,gsyn=1.4,iapp=8,fr=0.01,t0=1000,Cn=0.04

# stn parameters
p vl=-60.,vna=55.,vk=-80.,thetam=30,sm=15
p gl=2.25,gna=50.,gk=40.,tn=1.,th=0.05
p gahp=9.,gca=.5,vca=140.,k1=15.,eps=5e-05
p kca=22.5,thetas=39.,ss=2.,xp=1.,i=0.,ssca=8.
p thetah=-39.,sh=3.1,thetan=-32.,sn=-8.
p taun0=1.,taun1=100.,thn=80.,sigman=26.
p tauh0=1.,tauh1=500,thh=57.,sigmah=3.,phi=5.
p thetat=-63.,kt=-7.8,gt=.5,phir=2.
p thetar=-67.,kr=2.,taur0=7.1,taur1=17.5,thr=-68,sigmar=2.2
p alpha=5,beta=1.,ab=-30.,vsyn=-100.
p rth=.25,rsig=-.07,i1=0.,rho1=.5,a1=.9

# gp parameters
p gnag=120.,gkg=30.,gahpg=30.,gtg=.5,gcag=.1,glg=.1
p vnag=55.,vkg=-80.,vcag=120.,vlg=-55.
p thetasg=-57.,ksg=2.,thetas1g=-35.,ks1g=2.
p thetarg=-70.,krg=-2.,taurg=30.
p thetamg=-37.,sigmamg=10.
p thetang=-50.,sigmang=14.
p taun0g=.05,taun1g=.27,thng=-40.,sng=-12.
p thetahg=-58.,sigmahg=-12.
p tauh0g=.05,tauh1g=.27,thhg=-40.,shg=-12.
p k1g=30.,kcag=3.,epsg=0.0055
p phig=1.,deltang=.3,deltahg=.1
p gsyngg=.0,vsyngg=-80.
p gsyng=0.4,vsyng=0,alphag=2.,betag=.14,abg=-20
p ka=.33

# stn functions
sinf(v)=1./(1.+exp(-(v+thetas)/ss))
sinfca(v)=1./(1.+exp(-(v+thetas)/ssca))
minf(v)=1./(1.+exp(-(v+thetam)/sm))
hinf(v)=1./(1.+exp((v-thetah)/sh))
ninf(v)=1./(1.+exp((v-thetan)/sn))
taun(v)=taun0+taun1/(1+exp((v+thn)/sigman))
tauh(v)=tauh0+tauh1/(1+exp((v+thh)/sigmah))
rinf(v)=1/(1+exp((v-thetar)/kr))
taur(v)=taur0+taur1/(1+exp((v+thr)/sigmar))
tinf(v)=1/(1+exp((v-thetat)/kt))
rnew(r)=1/(1+exp((r-rth)/rsig))-1/(1+exp(-rth/rsig))

# gp functions
sinfg(vg)=1/(1+exp(-(vg-thetasg)/ksg))
sinf1g(vg)=1/(1+exp(-(vg-thetas1g)/ks1g))
rinfg(vg)=1/(1+exp(-(vg-thetarg)/krg))
minfg(vg)=1./(1.+exp(-(vg-thetamg)/sigmamg))
ninfg(vg)=1./(1.+exp(-(vg-thetang)/sigmang))
taung(vg)=taun0g+taun1g/(1+exp(-(vg-thng)/sng))
hinfg(vg)=1./(1.+exp(-(vg-thetahg)/sigmahg))
tauhg(vg)=tauh0g+tauh1g/(1+exp(-(vg-thhg)/shg))
hv(v)=1./(1.+exp(-v/.001))

# stn currents
il(v)=gl*(v-vl)
ina(v,h)=gna*(minf(v))^3*h*(v-vna)
ik(v,n)=gk*n^4*(v-vk)
iahp(v,ca)=gahp*(v-vk)*ca/(ca+k1)
ica(v)=gca*((sinfca(v))^2)*(v-vca)
it(v,r)=gt*(tinf(v)**3)*(rnew(r)^2)*(v-vca)

###################################################################
circular structure
###################################################################

isyn1=gsyn*(sg10+sg1+sg2)*(v1-vsyn)
isyn[2..9]=gsyn*(sg[j-1]+sg[j]+sg[j+1])*(v[j]-vsyn)
isyn10=gsyn*(sg9+sg10+sg1)*(v10-vsyn)

###################################################################

# gp currents
itg(vg,rg)=gtg*(sinfg(vg)^3)*rg*(vg-vcag)
inag(vg,hg)=gnag*(minfg(vg)^3)*hg*(vg-vnag)
ikg(vg,ng)=gkg*(ng^4)*(vg-vkg)
iahpg(vg,cag)=gahpg*(vg-vkg)*cag/(cag+k1g)
icag(vg)=gcag*((sinf1g(vg))^2)*(vg-vcag)
ilg(vg)=glg*(vg-vlg)

###########################################
isyng[1..10]=gsyng*s[j]*(vg[j]-vsyng)

###########################################

#stimulation current

C1=Cn*Heav(t-t0)
tau=0.5/fr
tau0=0.5/fr

Istim1=0.
Istim2=0.
Istim3=w2*Istim5
Istim4=w1*Istim5+w2*Istim6
Istim5=2.*Kn*(1./(1.+exp(-(C1*(delay(x1,tau0+tau)-delay(x1,tau0)))/14.))-0.5)+w1*Istim6+w2*Istim7
Istim6=2.*Kn*(1./(1.+exp(-(C1*(delay(x2,tau0+tau)-delay(x2,tau0)))/14.))-0.5)+w1*(Istim5+Istim7)
Istim7=2.*Kn*(1./(1.+exp(-(C1*(delay(x3,tau0+tau)-delay(x3,tau0)))/14.))-0.5)+w1*Istim6+w2*Istim5
Istim8=w1*Istim7+w2*Istim6
istim9=w2*Istim7
Istim10=0.

aux Stim5=Istim5
aux Stim6=Istim6
aux Stim7=Istim7

# LFP of the neuronal population,assumes recording with electrode at 5th STN neuron
Lfp5=0.1*(-isyn5+Istim5-w1*(isyn4+isyn6)-w2*(isyn3+isyn7))
Lfp6=0.1*(-isyn6+Istim6-w1*(isyn5+isyn7)-w2*(isyn4+isyn8))
Lfp7=0.1*(-isyn7+Istim7-w1*(isyn6+isyn8)-w2*(isyn5+isyn9))

aux Muf5=Lfp5
aux Muf6=Lfp6
aux Muf7=Lfp7

# stn initial conditions
i V1=-71.2,V2=-66.2,V3=-55.8,V4=-64.6,V5=-61.2,V6=-69.0,V7=-71.4,V8=-69.3,V9=-68.6,V10=-72.0
i H1=0.18,H2=0.35,H3=0.26,H4=0.28,H5=0.19,H6=0.31,H7=0.15,H8=0.23,H9=0.22,H10=0.34
i N1=0.33,N2=0.02,N3=0.1,N4=0.05,N5=0.26,N6=0.05,N7=0.44,N8=0.11,N9=0.11,N10=0.03
i R1=0.51,R2=0.57,R3=0.3,R4=0.63,R5=0.48,R6=0.54,R7=0.42,R8=0.62,R9=0.5,R10=0.61
i CA1=0.53,CA2=0.52,CA3=0.51,CA4=0.51,CA5=0.51,CA6=0.51,CA7=0.5,CA8=0.50,CA9=0.45,CA10=0.45
i S1=0.015,S2=0.,S3=0.,S4=0.,S5=0.003,S6=0.,S7=0.22,S8=0.,S9=0.,S10=0.

# gp initial conditions
i VG1=-66.42,VG2=-59.05,VG3=32.92,VG4=-59.91,VG5=-62.2,VG6=-59.7,VG7=-77.47,VG8=-61.83,VG9=-42.9,VG10=-74.73
i NG1=0.24,NG2=0.34,NG3=0.46,NG4=0.33,NG5=0.28,NG6=0.33,NG7=0.62,NG8=0.29,NG9=0.38,NG10=0.34
i HG1=0.55,HG2=0.52,HG3=0.51,HG4=0.55,HG5=0.58,HG6=0.52,HG7=0.29,HG8=0.59,HG9=0.54,HG10=0.42
i RG1=0.22,RG2=0.1,RG3=0.16,RG4=0.12,RG5=0.19,RG6=0.12,RG7=0.15,RG8=0.18,RG9=0.16,RG10=0.13
i CAG1=0.43,CAG2=0.28,CAG3=0.35,CAG4=0.31,CAG5=0.39,CAG6=0.31,CAG7=0.40,CAG8=0.39,CAG9=0.34,CAG10=0.33
i SG1=0.48,SG2=0.02,SG3=0.33,SG4=0.04,SG5=0.37,SG6=0.04,SG7=0.67,SG8=0.13,SG9=0.29,SG10=0.54

#delay initial conditions
i x1=0.,x2=0.,x3=0.

p Istn=10.,Igpe=0.

aux Mf5=0.001*delay(x1,tau0)
aux Mf6=0.001*delay(x2,tau0)
aux Mf7=0.001*delay(x3,tau0)

om=2.*pi*fr

# harmonic oscillators' equations

x1'=y1
y1'=-1.*om*y1-(om^2)*x1+Lfp5

x2'=y2
y2'=-1.*om*y2-(om^2)*x2+Lfp6

x3'=y3
y3'=-1.*om*y3-(om^2)*x3+Lfp7

# stn equations
v[1..10]'=-(il(v[j])+ina(v[j],h[j])+ik(v[j],n[j])+iahp(v[j],ca[j])+ica(v[j])+it(v[j],r[j]))-isyn[j]+Istn+Istim[j]
r[1..10]'=phir*(rinf(v[j])-r[j])/taur(v[j])
h[1..10]'=phi*( hinf(v[j])-h[j] )/tauh(v[j])
n[1..10]'=phi*( ninf(v[j])-n[j] )/taun(v[j])
ca[1..10]'=phi*eps*(-ica(v[j])-it(v[j],r[j]) - kca*ca[j])
s[1..10]'=alpha*(1-s[j])*sinf(v[j]+ab)-beta*s[j]

# gp equations
vg[1..10]'= -(itg(vg[j],rg[j])+inag(vg[j],hg[j])+ikg(vg[j],ng[j])+iahpg(vg[j],cag[j])+icag(vg[j])+ilg(vg[j]))+iapp-isyng[j]+Igpe
ng[1..10]'= deltang*(ninfg(vg[j])-ng[j])/taung(vg[j])
hg[1..10]'= deltahg*(hinfg(vg[j])-hg[j])/tauhg(vg[j])
rg[1..10]'= phig*(rinfg(vg[j])-rg[j])/taurg
cag[1..10]'=epsg*(-icag(vg[j])-itg(vg[j],rg[j]) - kcag*cag[j])
sg[1..10]'=alphag*(1-sg[j])*sinfg(vg[j]+abg)-betag*sg[j]


@ dt=.05,total=6999,meth=qualrk,xp=t,yp=v1,xlo=0,xhi=7000,ylo=-80,yhi=20.,bound=500001,maxstor=500001,delay=100.

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 [36] 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
sinf v ss thetas 1 1 v thetas ss
sinfca v ssca thetas 1 1 v thetas ssca
minf v sm thetam 1 1 v thetam sm
hinf v sh thetah 1 1 v thetah sh
ninf v sn thetan 1 1 v thetan sn
taun v sigman taun0 taun1 thn taun0 taun1 1 v thn sigman
tauh v sigmah tauh0 tauh1 thh tauh0 tauh1 1 v thh sigmah
rinf v kr thetar 1 1 v thetar kr
taur v sigmar taur0 taur1 thr taur0 taur1 1 v thr sigmar
tinf v kt thetat 1 1 v thetat kt
rnew r rsig rth 1 1 r rth rsig 1 1 rth rsig
sinfg vg ksg thetasg 1 1 vg thetasg ksg
sinf1g vg ks1g thetas1g 1 1 vg thetas1g ks1g
rinfg vg krg thetarg 1 1 vg thetarg krg
minfg vg sigmamg thetamg 1 1 vg thetamg sigmamg
ninfg vg sigmang thetang 1 1 vg thetang sigmang
taung vg sng taun0g taun1g thng taun0g taun1g 1 vg thng sng
hinfg vg sigmahg thetahg 1 1 vg thetahg sigmahg
tauhg vg shg tauh0g tauh1g thhg tauh0g tauh1g 1 vg thhg shg
hv v 1 1 v 0.001
il v gl vl gl v vl
ina v h gna sm thetam vna gna minf v thetam sm 3 h v vna
ik v n gk vk gk n 4 v vk
iahp v ca gahp k1 vk gahp v vk ca ca k1
ica v gca ssca thetas vca gca sinfca v thetas ssca 2 v vca
it v r gt kt rsig rth thetat vca gt tinf v kt thetat 3 rnew r rsig rth 2 v vca
itg vg rg gtg ksg thetasg vcag gtg sinfg vg ksg thetasg 3 rg vg vcag
inag vg hg gnag sigmamg thetamg vnag gnag minfg vg sigmamg thetamg 3 hg vg vnag
ikg vg ng gkg vkg gkg ng 4 vg vkg
iahpg vg cag gahpg k1g vkg gahpg vg vkg cag cag k1g
icag vg gcag ks1g thetas1g vcag gcag sinf1g vg thetas1g ks1g 2 vg vcag
ilg vg glg vlg glg vg vlg

Parameters [277] name constant value unit derived unit sbo cvterm
w1 w1 = 0.3 0.3 None
w2 w2 = 0 0.0 None
kn kn = 60 60.0 None
gsyn gsyn = 1.4 1.4 None
iapp iapp = 8 8.0 None
fr fr = 0.01 0.01 None
t0 t0 = 1000 1000.0 None
cn cn = 0.04 0.04 None
vl vl = -60. -60.0 None
vna vna = 55. 55.0 None
vk vk = -80. -80.0 None
thetam thetam = 30 30.0 None
sm sm = 15 15.0 None
gl gl = 2.25 2.25 None
gna gna = 50. 50.0 None
gk gk = 40. 40.0 None
tn tn = 1. 1.0 None
th th = 0.05 0.05 None
gahp gahp = 9. 9.0 None
gca gca = .5 0.5 None
vca vca = 140. 140.0 None
k1 k1 = 15. 15.0 None
eps eps = 5e-05 5e-05 None
kca kca = 22.5 22.5 None
thetas thetas = 39. 39.0 None
ss ss = 2. 2.0 None
xp xp = 1. 1.0 None
i i = 0. 0.0 None
ssca ssca = 8. 8.0 None
thetah thetah = -39. -39.0 None
sh sh = 3.1 3.1 None
thetan thetan = -32. -32.0 None
sn sn = -8. -8.0 None
taun0 taun0 = 1. 1.0 None
taun1 taun1 = 100. 100.0 None
thn thn = 80. 80.0 None
sigman sigman = 26. 26.0 None
tauh0 tauh0 = 1. 1.0 None
tauh1 tauh1 = 500 500.0 None
thh thh = 57. 57.0 None
sigmah sigmah = 3. 3.0 None
phi phi = 5. 5.0 None
thetat thetat = -63. -63.0 None
kt kt = -7.8 -7.8 None
gt gt = .5 0.5 None
phir phir = 2. 2.0 None
thetar thetar = -67. -67.0 None
kr kr = 2. 2.0 None
taur0 taur0 = 7.1 7.1 None
taur1 taur1 = 17.5 17.5 None
thr thr = -68 -68.0 None
sigmar sigmar = 2.2 2.2 None
alpha alpha = 5 5.0 None
beta beta = 1. 1.0 None
ab ab = -30. -30.0 None
vsyn vsyn = -100. -100.0 None
rth rth = .25 0.25 None
rsig rsig = -.07 -0.07 None
i1 i1 = 0. 0.0 None
rho1 rho1 = .5 0.5 None
a1 a1 = .9 0.9 None
gnag gnag = 120. 120.0 None
gkg gkg = 30. 30.0 None
gahpg gahpg = 30. 30.0 None
gtg gtg = .5 0.5 None
gcag gcag = .1 0.1 None
glg glg = .1 0.1 None
vnag vnag = 55. 55.0 None
vkg vkg = -80. -80.0 None
vcag vcag = 120. 120.0 None
vlg vlg = -55. -55.0 None
thetasg thetasg = -57. -57.0 None
ksg ksg = 2. 2.0 None
thetas1g thetas1g = -35. -35.0 None
ks1g ks1g = 2. 2.0 None
thetarg thetarg = -70. -70.0 None
krg krg = -2. -2.0 None
taurg taurg = 30. 30.0 None
thetamg thetamg = -37. -37.0 None
sigmamg sigmamg = 10. 10.0 None
thetang thetang = -50. -50.0 None
sigmang sigmang = 14. 14.0 None
taun0g taun0g = .05 0.05 None
taun1g taun1g = .27 0.27 None
thng thng = -40. -40.0 None
sng sng = -12. -12.0 None
thetahg thetahg = -58. -58.0 None
sigmahg sigmahg = -12. -12.0 None
tauh0g tauh0g = .05 0.05 None
tauh1g tauh1g = .27 0.27 None
thhg thhg = -40. -40.0 None
shg shg = -12. -12.0 None
k1g k1g = 30. 30.0 None
kcag kcag = 3. 3.0 None
epsg epsg = 0.0055 0.0055 None
phig phig = 1. 1.0 None
deltang deltang = .3 0.3 None
deltahg deltahg = .1 0.1 None
gsyngg gsyngg = .0 0.0 None
vsyngg vsyngg = -80. -80.0 None
gsyng gsyng = 0.4 0.4 None
vsyng vsyng = 0 0.0 None
alphag alphag = 2. 2.0 None
betag betag = .14 0.14 None
abg abg = -20 -20.0 None
ka ka = .33 0.33 None
v1 v1 = -71.2 -71.2 None
v2 v2 = -66.2 -66.2 None
v3 v3 = -55.8 -55.8 None
v4 v4 = -64.6 -64.6 None
v5 v5 = -61.2 -61.2 None
v6 v6 = -69.0 -69.0 None
v7 v7 = -71.4 -71.4 None
v8 v8 = -69.3 -69.3 None
v9 v9 = -68.6 -68.6 None
v10 v10 = -72.0 -72.0 None
h1 h1 = 0.18 0.18 None
h2 h2 = 0.35 0.35 None
h3 h3 = 0.26 0.26 None
h4 h4 = 0.28 0.28 None
h5 h5 = 0.19 0.19 None
h6 h6 = 0.31 0.31 None
h7 h7 = 0.15 0.15 None
h8 h8 = 0.23 0.23 None
h9 h9 = 0.22 0.22 None
h10 h10 = 0.34 0.34 None
n1 n1 = 0.33 0.33 None
n2 n2 = 0.02 0.02 None
n3 n3 = 0.1 0.1 None
n4 n4 = 0.05 0.05 None
n5 n5 = 0.26 0.26 None
n6 n6 = 0.05 0.05 None
n7 n7 = 0.44 0.44 None
n8 n8 = 0.11 0.11 None
n9 n9 = 0.11 0.11 None
n10 n10 = 0.03 0.03 None
r1 r1 = 0.51 0.51 None
r2 r2 = 0.57 0.57 None
r3 r3 = 0.3 0.3 None
r4 r4 = 0.63 0.63 None
r5 r5 = 0.48 0.48 None
r6 r6 = 0.54 0.54 None
r7 r7 = 0.42 0.42 None
r8 r8 = 0.62 0.62 None
r9 r9 = 0.5 0.5 None
r10 r10 = 0.61 0.61 None
ca1 ca1 = 0.53 0.53 None
ca2 ca2 = 0.52 0.52 None
ca3 ca3 = 0.51 0.51 None
ca4 ca4 = 0.51 0.51 None
ca5 ca5 = 0.51 0.51 None
ca6 ca6 = 0.51 0.51 None
ca7 ca7 = 0.5 0.5 None
ca8 ca8 = 0.50 0.5 None
ca9 ca9 = 0.45 0.45 None
ca10 ca10 = 0.45 0.45 None
s1 s1 = 0.015 0.015 None
s2 s2 = 0. 0.0 None
s3 s3 = 0. 0.0 None
s4 s4 = 0. 0.0 None
s5 s5 = 0.003 0.003 None
s6 s6 = 0. 0.0 None
s7 s7 = 0.22 0.22 None
s8 s8 = 0. 0.0 None
s9 s9 = 0. 0.0 None
s10 s10 = 0. 0.0 None
vg1 vg1 = -66.42 -66.42 None
vg2 vg2 = -59.05 -59.05 None
vg3 vg3 = 32.92 32.92 None
vg4 vg4 = -59.91 -59.91 None
vg5 vg5 = -62.2 -62.2 None
vg6 vg6 = -59.7 -59.7 None
vg7 vg7 = -77.47 -77.47 None
vg8 vg8 = -61.83 -61.83 None
vg9 vg9 = -42.9 -42.9 None
vg10 vg10 = -74.73 -74.73 None
ng1 ng1 = 0.24 0.24 None
ng2 ng2 = 0.34 0.34 None
ng3 ng3 = 0.46 0.46 None
ng4 ng4 = 0.33 0.33 None
ng5 ng5 = 0.28 0.28 None
ng6 ng6 = 0.33 0.33 None
ng7 ng7 = 0.62 0.62 None
ng8 ng8 = 0.29 0.29 None
ng9 ng9 = 0.38 0.38 None
ng10 ng10 = 0.34 0.34 None
hg1 hg1 = 0.55 0.55 None
hg2 hg2 = 0.52 0.52 None
hg3 hg3 = 0.51 0.51 None
hg4 hg4 = 0.55 0.55 None
hg5 hg5 = 0.58 0.58 None
hg6 hg6 = 0.52 0.52 None
hg7 hg7 = 0.29 0.29 None
hg8 hg8 = 0.59 0.59 None
hg9 hg9 = 0.54 0.54 None
hg10 hg10 = 0.42 0.42 None
rg1 rg1 = 0.22 0.22 None
rg2 rg2 = 0.1 0.1 None
rg3 rg3 = 0.16 0.16 None
rg4 rg4 = 0.12 0.12 None
rg5 rg5 = 0.19 0.19 None
rg6 rg6 = 0.12 0.12 None
rg7 rg7 = 0.15 0.15 None
rg8 rg8 = 0.18 0.18 None
rg9 rg9 = 0.16 0.16 None
rg10 rg10 = 0.13 0.13 None
cag1 cag1 = 0.43 0.43 None
cag2 cag2 = 0.28 0.28 None
cag3 cag3 = 0.35 0.35 None
cag4 cag4 = 0.31 0.31 None
cag5 cag5 = 0.39 0.39 None
cag6 cag6 = 0.31 0.31 None
cag7 cag7 = 0.40 0.4 None
cag8 cag8 = 0.39 0.39 None
cag9 cag9 = 0.34 0.34 None
cag10 cag10 = 0.33 0.33 None
sg1 sg1 = 0.48 0.48 None
sg2 sg2 = 0.02 0.02 None
sg3 sg3 = 0.33 0.33 None
sg4 sg4 = 0.04 0.04 None
sg5 sg5 = 0.37 0.37 None
sg6 sg6 = 0.04 0.04 None
sg7 sg7 = 0.67 0.67 None
sg8 sg8 = 0.13 0.13 None
sg9 sg9 = 0.29 0.29 None
sg10 sg10 = 0.54 0.54 None
x1 x1 = 0. 0.0 None
x2 x2 = 0. 0.0 None
x3 x3 = 0. 0.0 None
istn istn = 10. 10.0 None
igpe igpe = 0. 0.0 None
y1 0.0 dimensionless None
y2 0.0 dimensionless None
y3 0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
isyn1 0.0 dimensionless None
0.0 dimensionless None
isyn10 0.0 dimensionless None
0.0 dimensionless None
c1 0.0 dimensionless None
tau 0.0 dimensionless None
tau0 0.0 dimensionless None
istim1 0.0 dimensionless None
istim2 0.0 dimensionless None
istim3 0.0 dimensionless None
istim4 0.0 dimensionless None
istim5 0.0 dimensionless None
istim6 0.0 dimensionless None
istim7 0.0 dimensionless None
istim8 0.0 dimensionless None
istim9 0.0 dimensionless None
istim10 0.0 dimensionless None
stim5 0.0 dimensionless None
stim6 0.0 dimensionless None
stim7 0.0 dimensionless None
lfp5 0.0 dimensionless None
lfp6 0.0 dimensionless None
lfp7 0.0 dimensionless None
muf5 0.0 dimensionless None
muf6 0.0 dimensionless None
muf7 0.0 dimensionless None
mf5 0.0 dimensionless None
mf6 0.0 dimensionless None
mf7 0.0 dimensionless None
om 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [49]   assignment name derived units sbo cvterm
d x1/dt = y1 None
d y1/dt = 1 om y1 om 2 x1 lfp5 None
d x2/dt = y2 None
d y2/dt = 1 om y2 om 2 x2 lfp6 None
d x3/dt = y3 None
d y3/dt = 1 om y3 om 2 x3 lfp7 None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
isyn1 = gsyn sg10 sg1 sg2 v1 vsyn None
= None None
isyn10 = gsyn sg9 sg10 sg1 v10 vsyn None
= None None
c1 = cn heav t t0 None
tau = 0.5 fr None
tau0 = 0.5 fr None
istim1 = 0 None
istim2 = 0 None
istim3 = w2 istim5 None
istim4 = w1 istim5 w2 istim6 None
istim5 = 2 kn 1 1 c1 delay x1 tau0 tau delay x1 tau0 14 0.5 w1 istim6 w2 istim7 None
istim6 = 2 kn 1 1 c1 delay x2 tau0 tau delay x2 tau0 14 0.5 w1 istim5 istim7 None
istim7 = 2 kn 1 1 c1 delay x3 tau0 tau delay x3 tau0 14 0.5 w1 istim6 w2 istim5 None
istim8 = w1 istim7 w2 istim6 None
istim9 = w2 istim7 None
istim10 = 0 None
stim5 = istim5 None
stim6 = istim6 None
stim7 = istim7 None
lfp5 = 0.1 isyn5 istim5 w1 isyn4 isyn6 w2 isyn3 isyn7 None
lfp6 = 0.1 isyn6 istim6 w1 isyn5 isyn7 w2 isyn4 isyn8 None
lfp7 = 0.1 isyn7 istim7 w1 isyn6 isyn8 w2 isyn5 isyn9 None
muf5 = lfp5 None
muf6 = lfp6 None
muf7 = lfp7 None
mf5 = 0.001 delay x1 tau0 None
mf6 = 0.001 delay x2 tau0 None
mf7 = 0.001 delay x3 tau0 None
om = 2 fr None
t = time None