XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# Theta neuron model 50e and 20i cells with all-to-all # connectivity and periodic input to both the e- and i-cells # # Parameters par I_e=0.4,I_i=0 par tauz=50,gz=.5 par sige=0.5,sigi=.01 par gie=3,gei=.4,gii=0.15,gee=2, par taue=1,taui=3 # # Periodic stimuli description with parameters E(t)=heav(t-t_on)*heav(stim+t_on-t) PO(t)=heav(d-t) P(t)=PO(mod(t,per)) z'=(-z+p(t))/tau I_x(t)=Ampx*z*E(t) I_y(t)=Ampy*z*E(t) par Ampx=70,Ampy=10,d=1,per=25,tau=20 par t_on=0,stim=1000 # # ODEs for e-cells(x), i-cells(y) and adaptation(z) x[0..49]'=1-cos(x[j])+(1+cos(x[j]))*(I_e+I_x(t)-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+I_y(t)-gii*si+gei*se+sigi*wi[j]) # global 1 x[0..49]-pi {x[j]=-pi;out_put=1} global 1 y[0..19]-pi {y[j]=-pi;out_put=1} # 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 # #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 aux per_input=I_x(t) @ dt=0.02,nOut=10 @ meth=euler,total=1400,trans=400,maxstor=1000000 @ yp=se,xlo=400,xhi=1400,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:
- 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 [10] | name | math | sbo | cvterm |
---|---|---|---|---|---|
FunctionDefinition | max
|
minimum | |||
FunctionDefinition | min
|
maximum | |||
FunctionDefinition | heav
|
heavyside | |||
FunctionDefinition | mod
|
modulo | |||
FunctionDefinition | e
|
||||
FunctionDefinition | po
|
||||
FunctionDefinition | p
|
||||
FunctionDefinition | i_x
|
||||
FunctionDefinition | i_y
|
||||
FunctionDefinition | sd
|
type | Parameters [35] | name | constant | value | unit | derived unit | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Parameter | i_e
|
i_e = 0.4 | F | 0.4 | None | |||
Parameter | i_i
|
i_i = 0 | F | 0.0 | None | |||
Parameter | tauz
|
tauz = 50 | F | 50.0 | None | |||
Parameter | gz
|
gz = .5 | F | 0.5 | None | |||
Parameter | sige
|
sige = 0.5 | F | 0.5 | None | |||
Parameter | sigi
|
sigi = .01 | F | 0.01 | None | |||
Parameter | gie
|
gie = 3 | F | 3.0 | None | |||
Parameter | gei
|
gei = .4 | F | 0.4 | None | |||
Parameter | gii
|
gii = 0.15 | F | 0.15 | None | |||
Parameter | gee
|
gee = 2 | F | 2.0 | None | |||
Parameter | taue
|
taue = 1 | F | 1.0 | None | |||
Parameter | taui
|
taui = 3 | F | 3.0 | None | |||
Parameter | ampx
|
ampx = 70 | F | 70.0 | None | |||
Parameter | ampy
|
ampy = 10 | F | 10.0 | None | |||
Parameter | d
|
d = 1 | F | 1.0 | None | |||
Parameter | per
|
per = 25 | F | 25.0 | None | |||
Parameter | tau
|
tau = 20 | F | 20.0 | None | |||
Parameter | t_on
|
t_on = 0 | F | 0.0 | None | |||
Parameter | stim
|
stim = 1000 | F | 1000.0 | None | |||
Parameter | x[0..49] | F | 0.0 | None | ||||
Parameter | y[0..19] | F | 0.0 | None | ||||
Parameter | T | 0.0 | None | |||||
Parameter | T | 0.0 | None | |||||
Parameter | b
|
b = 100 | F | 100.0 | None | |||
Parameter | x[0..49] = ran(1)*2*pi-pi | T | 0.0 | dimensionless | None | |||
Parameter | y[0..19] = ran(1)*2*pi-pi | T | 0.0 | dimensionless | None | |||
Parameter | z
|
F | 0.0 | dimensionless | None | |||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | se
|
F | 0.0 | dimensionless | None | |||
Parameter | si
|
F | 0.0 | dimensionless | None | |||
Parameter | swgt
|
F | 0.0 | dimensionless | None | |||
Parameter | per_input
|
F | 0.0 | dimensionless | None | |||
Parameter | t
|
model time | F | 0.0 | dimensionless | None |
type | InitialAssignments [2] | name | assignment | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
InitialAssignment | = | None | |||||
InitialAssignment | = | None |
type | Rules [9] | assignment | name | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
Rule | d z/dt
|
= | None | ||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | d se/dt
|
= | None | ||||
Rule | d si/dt
|
= | None | ||||
Rule | swgt
|
= | None | ||||
Rule | per_input
|
= | None | ||||
Rule | t
|
= | None |
type | Events [2] | name | trigger | priority | delay | assignments | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Event | e0
|
None initialValue = False persistent = True |
=
out_put = |
|||||
Event | e1
|
None initialValue = False persistent = True |
=
out_put = |