#karin2



#reactions
v1:
 $pool = {1.0}z
  z * lambdaPlus
v2:
 {1.0}z = $pool
  z * lambdaMinus
v3:
 $pool = {1.0}y
  mu * (M - (z + zmut) * y)
v4:
 $pool = {1.0}zmut
  zmut * (k * y / 10 - 0.5)


#parameters
M = 25.0
k = 4.0
lambdaMinus = 0.5
mu = 0.25
default = 1.0


#initial values
y = 5.0
z = 5.0
zmut = 0.0


#assignment rules
!F lambdaPlus = y / 10.0