XPP model

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

# doublebvp.ode
# used for Bose and Rubin, IJBC, Vol. 25, June 30, 2015

# integrate 4 trajectories simultaneously

# (v,g) is the suppressed cell on its last cycle before it becomes active
v'=P*(I-v-g*(v-E)) 
g'=P*(-beta*g)

# (vhat,ghat) is the active cell on its last cycle of being active; ghat \equiv 0
vhat'=P*(I-vhat)

# (vtilde,gtilde) is the formerly suppressed cell on the cycle after its first reset
vtilde'=Q*(I-vtilde-gtilde*(vtilde-E))
gtilde'=Q*(-beta*gtilde)

# (vstar,gstar) is the formerly active cell on the cycle after it is first inhibited
vstar'=Q*(I-vstar-gstar*(vstar-E))
gstar'=Q*(-beta*gstar)

# two trivial equations are added, so that number of eqns = number of BCs
P'=0
Q'=0

# BCs
b v-1
b v'-1
b g-g0
b vhat
b vtilde
b vtilde'-1
b g'-gtilde
b vhat'-vstar
b vstar'-1

init P=0.65,Q=1.02,v=1,g=1.1,vstar=.96,gstar=1.24
par beta=0.2,I=2,E=-0.1,g0=1.1
@ total=1,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 [14] name constant value unit derived unit sbo cvterm
p p = 0.65 0.65 None
q q = 1.02 1.02 None
v v = 1 1.0 None
g g = 1.1 1.1 None
vstar vstar = .96 0.96 None
gstar gstar = 1.24 1.24 None
beta beta = 0.2 0.2 None
i i = 2 2.0 None
e e = -0.1 -0.1 None
g0 g0 = 1.1 1.1 None
vhat 0.0 dimensionless None
vtilde 0.0 dimensionless None
gtilde 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [10]   assignment name derived units sbo cvterm
d v/dt = p i v g v e None
d g/dt = p beta g None
d vhat/dt = p i vhat None
d vtilde/dt = q i vtilde gtilde vtilde e None
d gtilde/dt = q beta gtilde None
d vstar/dt = q i vstar gstar vstar e None
d gstar/dt = q beta gstar None
d p/dt = 0 None
d q/dt = 0 None
t = time None