C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 33 entries in each of the rate and state variable arrays.
C There are a total of 80 entries in the constant variable array.
C
C
C VOI is time in component environment (s).
C CONSTS(1) is kf1 in component PI3K (second_order_rate_constant).
C CONSTS(2) is kb1 in component PI3K (first_order_rate_constant).
C CONSTS(3) is kf2 in component PI3K (second_order_rate_constant).
C CONSTS(4) is kb2 in component PI3K (first_order_rate_constant).
C CONSTS(5) is kf3 in component PI3K (first_order_rate_constant).
C CONSTS(6) is kb3 in component PI3K (first_order_rate_constant).
C CONSTS(7) is kf34 in component PI3K (first_order_rate_constant).
C CONSTS(8) is kb34 in component PI3K (first_order_rate_constant).
C CONSTS(9) is V4 in component PI3K (flux).
C CONSTS(10) is k4 in component PI3K (nm).
C CONSTS(11) is kf5 in component PI3K (second_order_rate_constant).
C CONSTS(12) is kb5 in component PI3K (first_order_rate_constant).
C CONSTS(13) is kf6 in component PI3K (first_order_rate_constant).
C CONSTS(14) is kb6 in component PI3K (first_order_rate_constant).
C CONSTS(15) is kf7 in component PI3K (first_order_rate_constant).
C CONSTS(16) is kb7 in component PI3K (first_order_rate_constant).
C CONSTS(17) is kf8 in component PI3K (first_order_rate_constant).
C CONSTS(18) is kb8 in component PI3K (second_order_rate_constant).
C CONSTS(19) is kf9 in component PI3K (first_order_rate_constant).
C CONSTS(20) is kb9 in component PI3K (second_order_rate_constant).
C CONSTS(21) is V10 in component PI3K (flux).
C CONSTS(22) is k10 in component PI3K (nm).
C CONSTS(23) is kf23 in component PI3K (second_order_rate_constant).
C CONSTS(24) is kb23 in component PI3K (first_order_rate_constant).
C CONSTS(25) is kf24 in component PI3K (first_order_rate_constant).
C CONSTS(26) is kb24 in component PI3K (first_order_rate_constant).
C CONSTS(27) is kf25 in component PI3K (first_order_rate_constant).
C CONSTS(28) is kb25 in component PI3K (second_order_rate_constant).
C CONSTS(29) is V26 in component PI3K (flux).
C CONSTS(30) is k26 in component PI3K (nm).
C STATES(1) is R in component PI3K (nm).
C STATES(2) is Shc in component PI3K (nm).
C STATES(3) is PI3K in component PI3K (nm).
C STATES(4) is HRG in component PI3K (nm).
C STATES(5) is R_HRG in component PI3K (nm).
C STATES(6) is R_HRG2 in component PI3K (nm).
C STATES(7) is Internalisation in component PI3K (nm).
C STATES(8) is RP in component PI3K (nm).
C STATES(9) is R_Shc in component PI3K (nm).
C STATES(10) is R_ShP in component PI3K (nm).
C STATES(11) is ShP in component PI3K (nm).
C STATES(12) is R_ShGS in component PI3K (nm).
C STATES(13) is ShGS in component PI3K (nm).
C STATES(14) is GS in component PI3K (nm).
C STATES(15) is R_PI3K in component PI3K (nm).
C STATES(16) is R_PI3Kstar in component PI3K (nm).
C STATES(17) is PI3Kstar in component PI3K (nm).
C CONSTS(31) is two in component PI3K (dimensionless).
C STATES(18) is RasGTP in component RasGDPtoRasGTP (nm).
C CONSTS(32) is kf11 in component RasGDPtoRasGTP (first_order_rate_constant).
C CONSTS(33) is k11 in component RasGDPtoRasGTP (nm).
C CONSTS(34) is V12 in component RasGDPtoRasGTP (flux).
C CONSTS(35) is k12 in component RasGDPtoRasGTP (nm).
C STATES(19) is RasGDP in component RasGDPtoRasGTP (nm).
C STATES(20) is Akt_PIPP in component Akt (nm).
C STATES(21) is RAF_star in component RAF (nm).
C CONSTS(36) is kf13 in component RAF (first_order_rate_constant).
C CONSTS(37) is k13 in component RAF (nm).
C CONSTS(38) is kf14 in component RAF (first_order_rate_constant).
C CONSTS(39) is k14 in component RAF (nm).
C CONSTS(40) is E in component RAF (nm).
C STATES(22) is RAF in component RAF (nm).
C STATES(23) is MEKP in component MEK (nm).
C STATES(24) is MEKPP in component MEK (nm).
C CONSTS(41) is kf27 in component Akt (first_order_rate_constant).
C CONSTS(42) is k27 in component Akt (nm).
C CONSTS(43) is V28 in component Akt (flux).
C CONSTS(44) is k28 in component Akt (nm).
C CONSTS(45) is kf29 in component Akt (second_order_rate_constant).
C CONSTS(46) is kb29 in component Akt (first_order_rate_constant).
C CONSTS(47) is V30 in component Akt (flux).
C CONSTS(48) is k30 in component Akt (nm).
C CONSTS(49) is kf31 in component Akt (first_order_rate_constant).
C CONSTS(50) is k31 in component Akt (nm).
C CONSTS(51) is V32 in component Akt (flux).
C CONSTS(52) is k32 in component Akt (nm).
C CONSTS(53) is kf33 in component Akt (first_order_rate_constant).
C CONSTS(54) is k33 in component Akt (nm).
C CONSTS(55) is k16 in component Akt (nm).
C CONSTS(56) is k18 in component Akt (nm).
C STATES(25) is P in component Akt (nm).
C STATES(26) is PIP3 in component Akt (nm).
C STATES(27) is Akt in component Akt (nm).
C STATES(28) is Akt_PIP3 in component Akt (nm).
C STATES(29) is Akt_PIP in component Akt (nm).
C CONSTS(57) is PP2A in component Akt (nm).
C CONSTS(58) is one in component Akt (dimensionless).
C CONSTS(59) is PP2A in component MEK (nm).
C STATES(30) is MEK in component MEK (nm).
C CONSTS(60) is kf15 in component MEK (first_order_rate_constant).
C CONSTS(61) is k15 in component MEK (nm).
C CONSTS(62) is kf16 in component MEK (first_order_rate_constant).
C CONSTS(63) is k16 in component MEK (nm).
C CONSTS(64) is kf17 in component MEK (first_order_rate_constant).
C CONSTS(65) is k17 in component MEK (nm).
C CONSTS(66) is kf18 in component MEK (first_order_rate_constant).
C CONSTS(67) is k18 in component MEK (nm).
C CONSTS(68) is k31 in component MEK (nm).
C CONSTS(69) is k33 in component MEK (nm).
C CONSTS(70) is one in component MEK (dimensionless).
C CONSTS(71) is MKP3 in component ERK (nm).
C STATES(31) is ERK in component ERK (nm).
C STATES(32) is ERKP in component ERK (nm).
C STATES(33) is ERKPP in component ERK (nm).
C CONSTS(72) is kf19 in component ERK (first_order_rate_constant).
C CONSTS(73) is k19 in component ERK (nm).
C CONSTS(74) is kf20 in component ERK (first_order_rate_constant).
C CONSTS(75) is k20 in component ERK (nm).
C CONSTS(76) is kf21 in component ERK (first_order_rate_constant).
C CONSTS(77) is k21 in component ERK (nm).
C CONSTS(78) is kf22 in component ERK (first_order_rate_constant).
C CONSTS(79) is k22 in component ERK (nm).
C CONSTS(80) is one in component ERK (dimensionless).
C RATES(1) is d/dt R in component PI3K (nm).
C RATES(4) is d/dt HRG in component PI3K (nm).
C RATES(5) is d/dt R_HRG in component PI3K (nm).
C RATES(6) is d/dt R_HRG2 in component PI3K (nm).
C RATES(8) is d/dt RP in component PI3K (nm).
C RATES(7) is d/dt Internalisation in component PI3K (nm).
C RATES(9) is d/dt R_Shc in component PI3K (nm).
C RATES(2) is d/dt Shc in component PI3K (nm).
C RATES(10) is d/dt R_ShP in component PI3K (nm).
C RATES(14) is d/dt GS in component PI3K (nm).
C RATES(11) is d/dt ShP in component PI3K (nm).
C RATES(12) is d/dt R_ShGS in component PI3K (nm).
C RATES(13) is d/dt ShGS in component PI3K (nm).
C RATES(15) is d/dt R_PI3K in component PI3K (nm).
C RATES(3) is d/dt PI3K in component PI3K (nm).
C RATES(16) is d/dt R_PI3Kstar in component PI3K (nm).
C RATES(17) is d/dt PI3Kstar in component PI3K (nm).
C RATES(19) is d/dt RasGDP in component RasGDPtoRasGTP (nm).
C RATES(18) is d/dt RasGTP in component RasGDPtoRasGTP (nm).
C RATES(22) is d/dt RAF in component RAF (nm).
C RATES(21) is d/dt RAF_star in component RAF (nm).
C RATES(25) is d/dt P in component Akt (nm).
C RATES(26) is d/dt PIP3 in component Akt (nm).
C RATES(27) is d/dt Akt in component Akt (nm).
C RATES(28) is d/dt Akt_PIP3 in component Akt (nm).
C RATES(29) is d/dt Akt_PIP in component Akt (nm).
C RATES(20) is d/dt Akt_PIPP in component Akt (nm).
C RATES(30) is d/dt MEK in component MEK (nm).
C RATES(23) is d/dt MEKP in component MEK (nm).
C RATES(24) is d/dt MEKPP in component MEK (nm).
C RATES(31) is d/dt ERK in component ERK (nm).
C RATES(33) is d/dt ERKPP in component ERK (nm).
C RATES(32) is d/dt ERKP in component ERK (nm).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.0012
      CONSTS(2) = 0.00076
      CONSTS(3) = 0.01
      CONSTS(4) = 0.1
      CONSTS(5) = 1
      CONSTS(6) = 0.01
      CONSTS(7) = 0.001
      CONSTS(8) = 0
      CONSTS(9) = 62.5
      CONSTS(10) = 50
      CONSTS(11) = 0.1
      CONSTS(12) = 1
      CONSTS(13) = 20
      CONSTS(14) = 5
      CONSTS(15) = 60
      CONSTS(16) = 546
      CONSTS(17) = 2040
      CONSTS(18) = 15700
      CONSTS(19) = 40.8
      CONSTS(20) = 0
      CONSTS(21) = 0.0154
      CONSTS(22) = 340
      CONSTS(23) = 0.1
      CONSTS(24) = 2
      CONSTS(25) = 9.85
      CONSTS(26) = 0.0985
      CONSTS(27) = 45.8
      CONSTS(28) = 0.047
      CONSTS(29) = 2620
      CONSTS(30) = 3680
      STATES(1) = 80
      STATES(2) = 1000
      STATES(3) = 10
      STATES(4) = 10
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 0
      STATES(11) = 0
      STATES(12) = 0
      STATES(13) = 0
      STATES(14) = 10
      STATES(15) = 0
      STATES(16) = 0
      STATES(17) = 0
      CONSTS(31) = 2
      STATES(18) = 0
      CONSTS(32) = 0.222
      CONSTS(33) = 0.181
      CONSTS(34) = 0.289
      CONSTS(35) = 0.0571
      STATES(19) = 120
      STATES(20) = 0.0
      STATES(21) = 100
      CONSTS(36) = 1.53
      CONSTS(37) = 11.7
      CONSTS(38) = 0.00673
      CONSTS(39) = 8.07
      CONSTS(40) = 7
      STATES(22) = 0
      STATES(23) = 0
      STATES(24) = 0
      CONSTS(41) = 16.9
      CONSTS(42) = 39.1
      CONSTS(43) = 17000
      CONSTS(44) = 9.02
      CONSTS(45) = 507
      CONSTS(46) = 234
      CONSTS(47) = 20000
      CONSTS(48) = 80000
      CONSTS(49) = 0.107
      CONSTS(50) = 4.35
      CONSTS(51) = 20000
      CONSTS(52) = 80000
      CONSTS(53) = 0.211
      CONSTS(54) = 12
      CONSTS(55) = 2200
      CONSTS(56) = 60
      STATES(25) = 800
      STATES(26) = 0
      STATES(27) = 10
      STATES(28) = 0
      STATES(29) = 0
      CONSTS(57) = 11.4
      CONSTS(58) = 1
      CONSTS(59) = 11.4
      STATES(30) = 120
      CONSTS(60) = 3.5
      CONSTS(61) = 317
      CONSTS(62) = 0.058
      CONSTS(63) = 2200
      CONSTS(64) = 2.9
      CONSTS(65) = 317
      CONSTS(66) = 0.058
      CONSTS(67) = 60
      CONSTS(68) = 4.35
      CONSTS(69) = 12
      CONSTS(70) = 1
      CONSTS(71) = 2.4
      STATES(31) = 1000
      STATES(32) = 0
      STATES(33) = 0
      CONSTS(72) = 9.5
      CONSTS(73) = 146000
      CONSTS(74) = 0.3
      CONSTS(75) = 160
      CONSTS(76) = 16
      CONSTS(77) = 146000
      CONSTS(78) = 0.27
      CONSTS(79) = 60
      CONSTS(80) = 1
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) =  - CONSTS(1)*STATES(1)*STATES(4)+ CONSTS(2)*STATES(5)
      RATES(4) =  - CONSTS(1)*STATES(1)*STATES(4)+ CONSTS(2)*STATES(5)
      RATES(5) = ( CONSTS(1)*STATES(1)*STATES(4) -  CONSTS(2)*STATES(5)) -  CONSTS(31)*( CONSTS(3)*STATES(5)*STATES(5) -  CONSTS(4)*STATES(6))
      RATES(6) = (( CONSTS(3)*STATES(5)*STATES(5) -  CONSTS(4)*STATES(6)) - ( CONSTS(5)*STATES(6) -  CONSTS(6)*STATES(8)))+( CONSTS(9)*STATES(8))/(CONSTS(10)+STATES(8))
      RATES(8) = (((((( CONSTS(5)*STATES(6) -  CONSTS(6)*STATES(8)) - ( CONSTS(9)*STATES(8))/(CONSTS(10)+STATES(8))) - ( CONSTS(11)*STATES(8)*STATES(2) -  CONSTS(12)*STATES(9)))+( CONSTS(17)*STATES(10) -  CONSTS(18)*STATES(13)*STATES(8))) - ( CONSTS(23)*STATES(8)*STATES(3) -  CONSTS(24)*STATES(15)))+( CONSTS(27)*STATES(16) -  CONSTS(28)*STATES(8)*STATES(17))) - ( CONSTS(7)*STATES(8) -  CONSTS(8)*STATES(7))
      RATES(7) =  CONSTS(7)*STATES(8) -  CONSTS(8)*STATES(7)
      RATES(9) = ( CONSTS(11)*STATES(8)*STATES(2) -  CONSTS(12)*STATES(9)) - ( CONSTS(13)*STATES(9) -  CONSTS(14)*STATES(10))
      RATES(2) = - ( CONSTS(11)*STATES(8)*STATES(2) -  CONSTS(12)*STATES(9))+( CONSTS(21)*STATES(11))/(CONSTS(22)+STATES(11))
      RATES(10) = ( CONSTS(13)*STATES(9) -  CONSTS(14)*STATES(10)) - ( CONSTS(15)*STATES(10) -  CONSTS(16)*STATES(12))
      RATES(14) = - ( CONSTS(15)*STATES(10) -  CONSTS(16)*STATES(12))+( CONSTS(19)*STATES(13) -  CONSTS(20)*STATES(14)*STATES(11))
      RATES(11) = ( CONSTS(19)*STATES(13) -  CONSTS(20)*STATES(14)*STATES(11)) - ( CONSTS(21)*STATES(11))/(CONSTS(22)+STATES(11))
      RATES(12) = ( CONSTS(15)*STATES(10) -  CONSTS(16)*STATES(12)) - ( CONSTS(17)*STATES(10) -  CONSTS(18)*STATES(13)*STATES(8))
      RATES(13) = ( CONSTS(17)*STATES(10) -  CONSTS(18)*STATES(13)*STATES(8)) - ( CONSTS(19)*STATES(13) -  CONSTS(20)*STATES(14)*STATES(11))
      RATES(15) = ( CONSTS(23)*STATES(8)*STATES(3) -  CONSTS(24)*STATES(15)) - ( CONSTS(25)*STATES(15) -  CONSTS(26)*STATES(16))
      RATES(3) = - ( CONSTS(23)*STATES(8)*STATES(3) -  CONSTS(24)*STATES(15))+( CONSTS(29)*STATES(17))/(CONSTS(30)+STATES(17))
      RATES(16) = ( CONSTS(25)*STATES(15) -  CONSTS(26)*STATES(16)) - ( CONSTS(27)*STATES(16) -  CONSTS(28)*STATES(8)*STATES(17))
      RATES(17) = ( CONSTS(27)*STATES(16) -  CONSTS(28)*STATES(8)*STATES(17)) - ( CONSTS(29)*STATES(17))/(CONSTS(30)+STATES(17))
      RATES(19) = - (( CONSTS(32)*STATES(13)*STATES(19))/(CONSTS(33)+STATES(19)))+( CONSTS(34)*STATES(18))/(CONSTS(35)+STATES(18))
      RATES(18) = ( CONSTS(32)*STATES(13)*STATES(19))/(CONSTS(33)+STATES(19)) - ( CONSTS(34)*STATES(18))/(CONSTS(35)+STATES(18))
      RATES(22) = ( CONSTS(38)*(STATES(20)+CONSTS(40))*STATES(21))/(CONSTS(39)+STATES(21)) - ( CONSTS(36)*STATES(18)*STATES(22))/(CONSTS(37)+STATES(22))
      RATES(21) = - (( CONSTS(38)*(STATES(20)+CONSTS(40))*STATES(21))/(CONSTS(39)+STATES(21)))+( CONSTS(36)*STATES(18)*STATES(22))/(CONSTS(37)+STATES(22))
      RATES(25) = ( CONSTS(43)*STATES(26))/(CONSTS(44)+STATES(26)) - ( CONSTS(41)*STATES(17)*STATES(25))/(CONSTS(42)+STATES(25))
      RATES(26) = (- (( CONSTS(43)*STATES(26))/(CONSTS(44)+STATES(26)))+( CONSTS(41)*STATES(17)*STATES(25))/(CONSTS(42)+STATES(25))) - ( CONSTS(45)*STATES(26)*STATES(27) -  CONSTS(46)*STATES(28))
      RATES(27) = - ( CONSTS(45)*STATES(26)*STATES(27) -  CONSTS(46)*STATES(28))
      RATES(28) = (( CONSTS(45)*STATES(26)*STATES(27) -  CONSTS(46)*STATES(28)) - ( CONSTS(47)*STATES(28))/( CONSTS(48)*(CONSTS(58)+STATES(29)/CONSTS(52))+STATES(28)))+( CONSTS(49)*CONSTS(57)*STATES(29))/( CONSTS(50)*(CONSTS(58)+STATES(23)/CONSTS(55)+STATES(24)/CONSTS(56)+STATES(20)/CONSTS(54))+STATES(29))
      RATES(29) = ((( CONSTS(47)*STATES(28))/( CONSTS(48)*(CONSTS(58)+STATES(29)/CONSTS(52))+STATES(28)) - ( CONSTS(49)*CONSTS(57)*STATES(29))/( CONSTS(50)*(CONSTS(58)+STATES(23)/CONSTS(55)+STATES(24)/CONSTS(56)+STATES(20)/CONSTS(54))+STATES(29))) - ( CONSTS(51)*STATES(29))/( CONSTS(52)*(CONSTS(58)+STATES(28)/CONSTS(48))+STATES(29)))+( CONSTS(53)*CONSTS(57)*STATES(20))/( CONSTS(54)*(CONSTS(58)+STATES(23)/CONSTS(55)+STATES(24)/CONSTS(56)+STATES(29)/CONSTS(50))+STATES(20))
      RATES(20) = ( CONSTS(51)*STATES(29))/( CONSTS(52)*(CONSTS(58)+STATES(28)/CONSTS(48))+STATES(29)) - ( CONSTS(53)*CONSTS(57)*STATES(20))/( CONSTS(54)*(CONSTS(58)+STATES(23)/CONSTS(55)+STATES(24)/CONSTS(56)+STATES(29)/CONSTS(50))+STATES(20))
      RATES(30) = - (( CONSTS(60)*STATES(21)*STATES(30))/( CONSTS(61)*(CONSTS(70)+STATES(23)/CONSTS(65))+STATES(30)))+( CONSTS(62)*CONSTS(59)*STATES(23))/( CONSTS(63)*(CONSTS(70)+STATES(24)/CONSTS(67)+STATES(29)/CONSTS(68)+STATES(20)/CONSTS(69))+STATES(23))
      RATES(23) = ((( CONSTS(60)*STATES(21)*STATES(30))/( CONSTS(61)*(CONSTS(70)+STATES(23)/CONSTS(65))+STATES(30)) - ( CONSTS(62)*CONSTS(59)*STATES(23))/( CONSTS(63)*(CONSTS(70)+STATES(24)/CONSTS(67)+STATES(29)/CONSTS(68)+STATES(20)/CONSTS(69))+STATES(23))) - ( CONSTS(64)*STATES(21)*STATES(23))/( CONSTS(65)*(CONSTS(70)+STATES(30)/CONSTS(61))+STATES(23)))+( CONSTS(66)*CONSTS(59)*STATES(24))/( CONSTS(67)*(CONSTS(70)+STATES(23)/CONSTS(63)+STATES(29)/CONSTS(68)+STATES(20)/CONSTS(69))+STATES(24))
      RATES(24) = ( CONSTS(64)*STATES(21)*STATES(23))/( CONSTS(65)*(CONSTS(70)+STATES(30)/CONSTS(61))+STATES(23)) - ( CONSTS(66)*CONSTS(59)*STATES(24))/( CONSTS(67)*(CONSTS(70)+STATES(23)/CONSTS(63)+STATES(29)/CONSTS(68)+STATES(20)/CONSTS(69))+STATES(24))
      RATES(31) = - (( CONSTS(72)*STATES(24)*STATES(31))/( CONSTS(73)*(CONSTS(80)+STATES(32)/CONSTS(77))+STATES(31)))+( CONSTS(74)*CONSTS(71)*STATES(32))/( CONSTS(75)*(CONSTS(80)+STATES(33)/CONSTS(79))+STATES(32))
      RATES(33) = ( CONSTS(76)*STATES(24)*STATES(32))/( CONSTS(77)*(CONSTS(80)+STATES(31)/CONSTS(73))+STATES(32)) - ( CONSTS(78)*CONSTS(71)*STATES(33))/( CONSTS(79)*(CONSTS(80)+STATES(32)/CONSTS(75))+STATES(33))
      RATES(32) = ((( CONSTS(72)*STATES(24)*STATES(31))/( CONSTS(73)*(CONSTS(80)+STATES(32)/CONSTS(77))+STATES(31)) - ( CONSTS(74)*CONSTS(71)*STATES(32))/( CONSTS(75)*(CONSTS(80)+STATES(33)/CONSTS(79))+STATES(32))) - ( CONSTS(76)*STATES(24)*STATES(32))/( CONSTS(77)*(CONSTS(80)+STATES(31)/CONSTS(73))+STATES(32)))+( CONSTS(78)*CONSTS(71)*STATES(33))/( CONSTS(79)*(CONSTS(80)+STATES(32)/CONSTS(75))+STATES(33))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END