XPP model

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

# Rate-based model of coupled neuron populations to 
# perform discrimination using integral feedback control.
# Continuous attractor version.
par w1_2=1.0 w2_1=-1.0 i0=0.5 tau=0.05 di=0.001 
par tstim=0.5 t1=5.0 t2=10.0 tauapp=0.005 dtapp=0.001
par f1=20 f2=16 

init r1=0.0 r2=0.0 
# init sapp1=0.0 sapp2=0.0 
init s1=0.0 s2=0.0

sapp1=f1/((1+exp(-2*(t-t1)/dtapp))*(1+exp(-2*(t1+tstim-t)/dtapp)))
sapp2=f2/((1+exp(-2*(t-t2)/dtapp))*(1+exp(-2*(t2+tstim-t)/dtapp)))

ds1/dt=(sapp1-s1)/tauapp
ds2/dt=(sapp2-s2)/tauapp

dr2/dt=w1_2*r1/tau
i1=w2_1*r2+s1+s2-i0
dr1/dt=(-r1+i1/(1+exp(-2*i1/di)))/tau
# dsapp1/dt=(s1-sapp1)/tauapp
# dsapp2/dt=(s2-sapp2)/tauapp
# ds1/dt=0.0
# ds2/dt=0.0

# global 1 {t-t1} {s1=f1}
# global 2 {t-t1-tstim} {s1=0.0}
# global 3 {t-t2} {s2=f2}
# global 4 {t-t2-tstim} {s2=0}

@total=15,bound=100,dt=.001,dtmin=1e-5,dtmax=10,atoler=1e-4
@toler=1e-5,xhi=15,yhi=2,ylo=0 njmp=50,bell=0
@bell=off,nout=50                                                             
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 [20] name constant value unit derived unit sbo cvterm
w1_2 w1_2 = 1.0 1.0 None
w2_1 w2_1 = -1.0 -1.0 None
i0 i0 = 0.5 0.5 None
tau tau = 0.05 0.05 None
di di = 0.001 0.001 None
tstim tstim = 0.5 0.5 None
t1 t1 = 5.0 5.0 None
t2 t2 = 10.0 10.0 None
tauapp tauapp = 0.005 0.005 None
dtapp dtapp = 0.001 0.001 None
f1 f1 = 20 20.0 None
f2 f2 = 16 16.0 None
r1 r1 = 0.0 0.0 None
r2 r2 = 0.0 0.0 None
s1 s1 = 0.0 0.0 None
s2 s2 = 0.0 0.0 None
sapp1 0.0 dimensionless None
sapp2 0.0 dimensionless None
i1 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [8]   assignment name derived units sbo cvterm
d s1/dt = sapp1 s1 tauapp None
d s2/dt = sapp2 s2 tauapp None
d r2/dt = w1_2 r1 tau None
d r1/dt = r1 i1 1 2 i1 di tau None
sapp1 = f1 1 2 t t1 dtapp 1 2 t1 tstim t dtapp None
sapp2 = f2 1 2 t t2 dtapp 1 2 t2 tstim t dtapp None
i1 = w2_1 r2 s1 s2 i0 None
t = time None