#nyman3



#reactions
v1a:
 {1.0}IR = {1.0}IRins
  k1a * ins * IR + k1aBasic * IR
v1b:
 {1.0}IRins = {1.0}IR
  k1b * IRins
v1c:
 {1.0}IRins = {1.0}IRp
  k1c * IRins
v1d:
 {1.0}IRp = {1.0}IRiP
  k1d * IRp
v1e:
 {1.0}IRiP = {1.0}IRi
  IRiP * (k1e + k1f * Xp / (1 + Xp))
v1g:
 {1.0}IRp = {1.0}IR
  k1g * IRp
v1r:
 {1.0}IRi = {1.0}IR
  k1r * IRi
v2:
 {1.0}IRS = {1.0}IRSiP
  k21 * IRS * (IRp + k22 * IRiP) / (1 + km23 * Xp)
v3:
 {1.0}X = {1.0}Xp
  k3 * X * IRSiP
vm2:
 {1.0}IRSiP = {1.0}IRS
  km2 * IRSiP
vm3:
 {1.0}Xp = {1.0}X
  km3 * Xp


#parameters
ins = 100.0
k1a = 0.153418
k1aBasic = 0.0383389
k1b = 3.4699e-06
k1c = 0.574266
k1d = 4.78844
k1e = 5.25027e-05
k1f = 119.353
k1g = 4.14899
k1r = 37954.7
k21 = 538004.0
k22 = 1.7252e-06
k3 = 8.62917e-05
km2 = 262759.0
km23 = 88.9096
km3 = 0.132671
default = 1.0


#initial values
IR = 8.94067597532632
IRS = 9.43998194225544
IRSiP = 0.560018057744573
IRi = 4.83863890758515e-06
IRiP = 0.424076631823384
IRins = 0.59688996214639
IRp = 0.0383525925240207
X = 9.99635886407151
Xp = 0.00364113592848386


#assignment rules
!F IRmem = (IRp + IRins + IR) * 10.0
!F measIRp = (IRp + IRiP) * 31.8
!F measIRS1 = IRSiP * 21.2