XPP model

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

# 200 e and 50 I integrate & fire
# random applied current, random conductances
# 
# prevent tables from being recomputed each time params change
@ autoeval=0
table wee % 40000 0 39999 ran(1)/100
table wei % 10000 0 9999  ran(1)/100
table wie % 10000 0 9999  ran(1)/25
table wii % 2500 0 2499  ran(1)/25
# multiply by the synapses
special see=mmult(200,200,wee,se0)
special sei=mmult(200,50,wei,se0)
special sie=mmult(50,200,wie,si0)
special sii=mmult(50,50,wii,si0)
# random currents
table r_e % 200 0 199  ran(1)-.5
table r_i % 50 0 49 ran(1)-.5
# parameters
par tau_e=20,tau_i=20,ele=-65,eli=-65
par vte=-50,vti=-50
par ver=-60,vir=-70
par taue=4,taui=10
# ODEs
ve[0..199]'=(-(ve[j]-ele) + ie0+ie1*r_e([j])-gee*see([j])*(ve[j]-eex)-gie*sie([j])*(ve[j]-ein))/tau_e
vi[0..49]'=(-(vi[j]-eli) + ii0+ii1*r_i([j])-gei*sei([j])*(vi[j]-eex)-gii*sii([j])*(ve[j]-ein))/tau_i
se[0..199]'=-se[j]/taue
si[0..49]'=-si[j]/taui
# each time votlage crosses threshold, reset and set synapse to 1
global 1 ve[0..199]-vte {se[j]=1;ve[j]=ver}
global 1 vi[0..49]-vti {si[j]=1;vi[j]=vir}
# more parameters
par ie0=10,ie1=.5
par ii0=10,ii1=.5
par gee=.02,gie=.05,gii=.02,gei=.02
par eex=0,ein=-75
# initial data
init ve[0..199]=-65
init vi[0..49]=-65
# numerical stuff
@ total=200,meth=euler,nout=10,dt=.01
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 [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
tau_e tau_e = 20 20.0 None
tau_i tau_i = 20 20.0 None
ele ele = -65 -65.0 None
eli eli = -65 -65.0 None
vte vte = -50 -50.0 None
vti vti = -50 -50.0 None
ver ver = -60 -60.0 None
vir vir = -70 -70.0 None
taue taue = 4 4.0 None
taui taui = 10 10.0 None
ie0 ie0 = 10 10.0 None
ie1 ie1 = .5 0.5 None
ii0 ii0 = 10 10.0 None
ii1 ii1 = .5 0.5 None
gee gee = .02 0.02 None
gie gie = .05 0.05 None
gii gii = .02 0.02 None
gei gei = .02 0.02 None
eex eex = 0 0.0 None
ein ein = -75 -75.0 None
ve[0..199] = -65 -65.0 None
vi[0..49] = -65 -65.0 None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [5]   assignment name derived units sbo cvterm
= None None
= None None
= None None
= None None
t = time None

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