XPP model

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

# Theta neuron network with 50e and 20i cells with all-to-all
# connectivity and tonic input to both the e- and i-cells
#
# Parameters
par I_e=1, I_i=0
par gei=.4,gii=0.15,gee=0.1, gie=2
par tauz=50,gz=.2
par sige=0.5, sigi=.02
par taue=1, taui=3
#
#ODEs for e-cells, i-cells and adaptation
x[0..49]'=1-cos(x[j])+(1+cos(x[j]))*(I_e-gz*z[j]-gie*si+gee*se+sige*we[j])
z[0..49]'=sd(x[j])-z[j]/tauz
y[0..19]'=1-cos(y[j])+(1+cos(y[j]))*(I_i-gii*si+gei*se+sigi*wi[j])
# Synapses
se'=sum(0,49)of(sd(shift(x0,i')))/50-se/taue
si'=sum(0,19)of(sd(shift(y0,i')))/20-si/taui
global 1 x[0..49]-pi {x[j]=-pi}
global 1 y[0..19]-pi {y[j]=-pi}
#
# Initial Conditions
x[0..49](0)=ran(1)*2*pi-pi
y[0..19](0)=ran(1)*2*pi-pi
wiener we[0..49]
wiener wi[0..19]
sd(x)=exp(-b*(1-cos(x-2.5)))*b
par b=100
aux swgt=.8*se+.2*si
@ meth=euler,
@ total=1400,trans=400,maxstor=1000000
@ dt=0.02,nOut=10
@ xp=t,yp=se,xlo=600,xhi=1200,ylo=0,yhi=.5
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 [5] 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
sd x b b 1 x 2.5 b

Parameters [26] name constant value unit derived unit sbo cvterm
i_e i_e = 1 1.0 None
i_i i_i = 0 0.0 None
gei gei = .4 0.4 None
gii gii = 0.15 0.15 None
gee gee = 0.1 0.1 None
gie gie = 2 2.0 None
tauz tauz = 50 50.0 None
gz gz = .2 0.2 None
sige sige = 0.5 0.5 None
sigi sigi = .02 0.02 None
taue taue = 1 1.0 None
taui taui = 3 3.0 None
x[0..49] 0.0 None
y[0..19] 0.0 None
0.0 None
0.0 None
b b = 100 100.0 None
x[0..49] = ran(1)*2*pi-pi 0.0 dimensionless None
y[0..19] = ran(1)*2*pi-pi 0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
se 0.0 dimensionless None
si 0.0 dimensionless None
swgt 0.0 dimensionless None
t model time 0.0 dimensionless None

InitialAssignments [2] name assignment derived units sbo cvterm
= ran 1 2 None
= ran 1 2 None

Rules [7]   assignment name derived units sbo cvterm
= None None
= None None
= None None
d se/dt = None None
d si/dt = None None
swgt = 0.8 se 0.2 si None
t = time None

Events [2] name trigger priority delay assignments sbo cvterm
e0 None
initialValue = False
persistent = True
=
e1 None
initialValue = False
persistent = True
=