C
C There are a total of 32 entries in the algebraic variable array.
C There are a total of 23 entries in each of the rate and state variable arrays.
C There are a total of 50 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is EGF in component EGF (nanomolar).
C ALGBRC(1) is v1 in component v1 (flux).
C STATES(2) is R in component R (nanomolar).
C STATES(3) is Ra in component Ra (nanomolar).
C ALGBRC(9) is v2 in component v2 (flux).
C STATES(4) is R2 in component R2 (nanomolar).
C ALGBRC(10) is v3 in component v3 (flux).
C ALGBRC(11) is v4 in component v4 (flux).
C STATES(5) is RP in component RP (nanomolar).
C ALGBRC(14) is v7 in component v7 (flux).
C ALGBRC(19) is v11 in component v11 (flux).
C ALGBRC(23) is v15 in component v15 (flux).
C ALGBRC(26) is v18 in component v18 (flux).
C ALGBRC(28) is v20 in component v20 (flux).
C ALGBRC(12) is v5 in component v5 (flux).
C ALGBRC(15) is v9 in component v9 (flux).
C ALGBRC(21) is v13 in component v13 (flux).
C STATES(6) is R_PL in component R_PL (nanomolar).
C ALGBRC(13) is v6 in component v6 (flux).
C STATES(7) is R_PLP in component R_PLP (nanomolar).
C STATES(8) is R_G in component R_G (nanomolar).
C ALGBRC(17) is v10 in component v10 (flux).
C STATES(9) is R_G_S in component R_G_S (nanomolar).
C STATES(10) is R_Sh in component R_Sh (nanomolar).
C ALGBRC(22) is v14 in component v14 (flux).
C STATES(11) is R_ShP in component R_ShP (nanomolar).
C ALGBRC(32) is v24 in component v24 (flux).
C ALGBRC(25) is v17 in component v17 (flux).
C STATES(12) is R_Sh_G in component R_Sh_G (nanomolar).
C ALGBRC(27) is v19 in component v19 (flux).
C STATES(13) is R_Sh_G_S in component R_Sh_G_S (nanomolar).
C STATES(14) is G_S in component G_S (nanomolar).
C ALGBRC(31) is v23 in component v23 (flux).
C ALGBRC(20) is v12 in component v12 (flux).
C STATES(15) is ShP in component ShP (nanomolar).
C ALGBRC(29) is v21 in component v21 (flux).
C ALGBRC(24) is v16 in component v16 (flux).
C STATES(16) is Sh_G in component Sh_G (nanomolar).
C ALGBRC(30) is v22 in component v22 (flux).
C STATES(17) is Sh_G_S in component Sh_G_S (nanomolar).
C STATES(18) is PLC_gamma in component PLC_gamma (nanomolar).
C ALGBRC(16) is v8 in component v8 (flux).
C STATES(19) is PLC_gamma_P in component PLC_gamma_P (nanomolar).
C ALGBRC(18) is v25 in component v25 (flux).
C STATES(20) is PLC_gamma_P_I in component PLC_gamma_P_I (nanomolar).
C STATES(21) is Grb in component Grb (nanomolar).
C STATES(22) is Shc in component Shc (nanomolar).
C STATES(23) is SOS in component SOS (nanomolar).
C CONSTS(1) is k1 in component v1 (second_order_rate_constant).
C CONSTS(2) is k1_ in component v1 (first_order_rate_constant).
C CONSTS(3) is k2 in component v2 (second_order_rate_constant).
C CONSTS(4) is k2_ in component v2 (first_order_rate_constant).
C CONSTS(5) is k3 in component v3 (first_order_rate_constant).
C CONSTS(6) is k3_ in component v3 (first_order_rate_constant).
C CONSTS(7) is K4 in component v4 (nanomolar).
C CONSTS(8) is V4 in component v4 (flux).
C CONSTS(9) is k5 in component v5 (second_order_rate_constant).
C CONSTS(10) is k5_ in component v5 (first_order_rate_constant).
C CONSTS(11) is k6 in component v6 (first_order_rate_constant).
C CONSTS(12) is k6_ in component v6 (first_order_rate_constant).
C CONSTS(13) is k7 in component v7 (first_order_rate_constant).
C CONSTS(14) is k7_ in component v7 (second_order_rate_constant).
C CONSTS(15) is K8 in component v8 (nanomolar).
C CONSTS(16) is V8 in component v8 (flux).
C CONSTS(17) is k9 in component v9 (second_order_rate_constant).
C CONSTS(18) is k9_ in component v9 (first_order_rate_constant).
C CONSTS(19) is k10 in component v10 (second_order_rate_constant).
C CONSTS(20) is k10_ in component v10 (first_order_rate_constant).
C CONSTS(21) is k11 in component v11 (first_order_rate_constant).
C CONSTS(22) is k11_ in component v11 (second_order_rate_constant).
C CONSTS(23) is k12 in component v12 (first_order_rate_constant).
C CONSTS(24) is k12_ in component v12 (second_order_rate_constant).
C CONSTS(25) is k13 in component v13 (second_order_rate_constant).
C CONSTS(26) is k13_ in component v13 (first_order_rate_constant).
C CONSTS(27) is k14 in component v14 (first_order_rate_constant).
C CONSTS(28) is k14_ in component v14 (first_order_rate_constant).
C CONSTS(29) is k15 in component v15 (first_order_rate_constant).
C CONSTS(30) is k15_ in component v15 (second_order_rate_constant).
C CONSTS(31) is K16 in component v16 (nanomolar).
C CONSTS(32) is V16 in component v16 (flux).
C CONSTS(33) is k17 in component v17 (second_order_rate_constant).
C CONSTS(34) is k17_ in component v17 (first_order_rate_constant).
C CONSTS(35) is k18 in component v18 (first_order_rate_constant).
C CONSTS(36) is k18_ in component v18 (second_order_rate_constant).
C CONSTS(37) is k19 in component v19 (second_order_rate_constant).
C CONSTS(38) is k19_ in component v19 (first_order_rate_constant).
C CONSTS(39) is k20 in component v20 (first_order_rate_constant).
C CONSTS(40) is k20_ in component v20 (second_order_rate_constant).
C CONSTS(41) is k21 in component v21 (second_order_rate_constant).
C CONSTS(42) is k21_ in component v21 (first_order_rate_constant).
C CONSTS(43) is k22 in component v22 (second_order_rate_constant).
C CONSTS(44) is k22_ in component v22 (first_order_rate_constant).
C CONSTS(45) is k23 in component v23 (first_order_rate_constant).
C CONSTS(46) is k23_ in component v23 (second_order_rate_constant).
C CONSTS(47) is k24 in component v24 (second_order_rate_constant).
C CONSTS(48) is k24_ in component v24 (first_order_rate_constant).
C CONSTS(49) is k25 in component v25 (first_order_rate_constant).
C CONSTS(50) is k25_ in component v25 (first_order_rate_constant).
C ALGBRC(2) is totEGFRphos in component combined_concs (nanomolar).
C ALGBRC(3) is totPLCgammaphos in component combined_concs (nanomolar).
C ALGBRC(4) is totGrb_EGFR in component combined_concs (nanomolar).
C ALGBRC(5) is totGrb_Shc in component combined_concs (nanomolar).
C ALGBRC(6) is totShcphos in component combined_concs (nanomolar).
C ALGBRC(7) is totShc_EGFR in component combined_concs (nanomolar).
C ALGBRC(8) is totSOS_EGFR in component combined_concs (nanomolar).
C RATES(1) is d/dt EGF in component EGF (nanomolar).
C RATES(2) is d/dt R in component R (nanomolar).
C RATES(3) is d/dt Ra in component Ra (nanomolar).
C RATES(4) is d/dt R2 in component R2 (nanomolar).
C RATES(5) is d/dt RP in component RP (nanomolar).
C RATES(6) is d/dt R_PL in component R_PL (nanomolar).
C RATES(7) is d/dt R_PLP in component R_PLP (nanomolar).
C RATES(8) is d/dt R_G in component R_G (nanomolar).
C RATES(9) is d/dt R_G_S in component R_G_S (nanomolar).
C RATES(10) is d/dt R_Sh in component R_Sh (nanomolar).
C RATES(11) is d/dt R_ShP in component R_ShP (nanomolar).
C RATES(12) is d/dt R_Sh_G in component R_Sh_G (nanomolar).
C RATES(13) is d/dt R_Sh_G_S in component R_Sh_G_S (nanomolar).
C RATES(14) is d/dt G_S in component G_S (nanomolar).
C RATES(15) is d/dt ShP in component ShP (nanomolar).
C RATES(16) is d/dt Sh_G in component Sh_G (nanomolar).
C RATES(17) is d/dt Sh_G_S in component Sh_G_S (nanomolar).
C RATES(18) is d/dt PLC_gamma in component PLC_gamma (nanomolar).
C RATES(19) is d/dt PLC_gamma_P in component PLC_gamma_P (nanomolar).
C RATES(20) is d/dt PLC_gamma_P_I in component PLC_gamma_P_I (nanomolar).
C RATES(21) is d/dt Grb in component Grb (nanomolar).
C RATES(22) is d/dt Shc in component Shc (nanomolar).
C RATES(23) is d/dt SOS in component SOS (nanomolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 680
      STATES(2) = 100
      STATES(3) = 0
      STATES(4) = 0
      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) = 0
      STATES(15) = 0
      STATES(16) = 0
      STATES(17) = 0
      STATES(18) = 105
      STATES(19) = 0
      STATES(20) = 0
      STATES(21) = 85
      STATES(22) = 150
      STATES(23) = 34
      CONSTS(1) = 0.003
      CONSTS(2) = 0.06
      CONSTS(3) = 0.01
      CONSTS(4) = 0.1
      CONSTS(5) = 1
      CONSTS(6) = 0.01
      CONSTS(7) = 50
      CONSTS(8) = 450
      CONSTS(9) = 0.06
      CONSTS(10) = 0.2
      CONSTS(11) = 1
      CONSTS(12) = 0.05
      CONSTS(13) = 0.3
      CONSTS(14) = 0.006
      CONSTS(15) = 100
      CONSTS(16) = 1
      CONSTS(17) = 0.003
      CONSTS(18) = 0.05
      CONSTS(19) = 0.01
      CONSTS(20) = 0.06
      CONSTS(21) = 0.03
      CONSTS(22) = 4.5e-3
      CONSTS(23) = 1.5e-3
      CONSTS(24) = 1e-4
      CONSTS(25) = 0.09
      CONSTS(26) = 0.6
      CONSTS(27) = 6
      CONSTS(28) = 0.06
      CONSTS(29) = 0.3
      CONSTS(30) = 9e-4
      CONSTS(31) = 340
      CONSTS(32) = 1.7
      CONSTS(33) = 0.003
      CONSTS(34) = 0.1
      CONSTS(35) = 0.3
      CONSTS(36) = 9e-4
      CONSTS(37) = 0.01
      CONSTS(38) = 2.14e-2
      CONSTS(39) = 0.12
      CONSTS(40) = 2.4e-4
      CONSTS(41) = 0.003
      CONSTS(42) = 0.1
      CONSTS(43) = 0.03
      CONSTS(44) = 0.064
      CONSTS(45) = 0.1
      CONSTS(46) = 0.021
      CONSTS(47) = 0.009
      CONSTS(48) = 4.29e-2
      CONSTS(49) = 1
      CONSTS(50) = 0.03
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(2)*STATES(1) -  CONSTS(2)*STATES(3)
      RATES(1) = - ALGBRC(1)
      RATES(2) = - ALGBRC(1)
      ALGBRC(9) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      RATES(3) = ALGBRC(1) -  2.00000*ALGBRC(9)
      ALGBRC(10) =  CONSTS(5)*STATES(4) -  CONSTS(6)*STATES(5)
      ALGBRC(11) = ( CONSTS(8)*STATES(5))/(CONSTS(7)+STATES(5))
      RATES(4) = (ALGBRC(9)+ALGBRC(11)) - ALGBRC(10)
      ALGBRC(12) =  CONSTS(9)*STATES(5)*STATES(18) -  CONSTS(10)*STATES(6)
      ALGBRC(13) =  CONSTS(11)*STATES(6) -  CONSTS(12)*STATES(7)
      RATES(6) = ALGBRC(12) - ALGBRC(13)
      ALGBRC(14) =  CONSTS(13)*STATES(7) -  CONSTS(14)*STATES(5)*STATES(19)
      RATES(7) = ALGBRC(13) - ALGBRC(14)
      ALGBRC(16) = ( CONSTS(16)*STATES(19))/(CONSTS(15)+STATES(19))
      RATES(18) = ALGBRC(16) - ALGBRC(12)
      ALGBRC(15) =  CONSTS(17)*STATES(5)*STATES(21) -  CONSTS(18)*STATES(8)
      ALGBRC(17) =  CONSTS(19)*STATES(8)*STATES(23) -  CONSTS(20)*STATES(9)
      RATES(8) = ALGBRC(15) - ALGBRC(17)
      ALGBRC(18) =  CONSTS(49)*STATES(19) -  CONSTS(50)*STATES(20)
      RATES(19) = ALGBRC(14) - (ALGBRC(16)+ALGBRC(18))
      RATES(20) = ALGBRC(18)
      ALGBRC(19) =  CONSTS(21)*STATES(9) -  CONSTS(22)*STATES(5)*STATES(14)
      RATES(9) = ALGBRC(17) - ALGBRC(19)
      ALGBRC(21) =  CONSTS(25)*STATES(5)*STATES(22) -  CONSTS(26)*STATES(10)
      ALGBRC(22) =  CONSTS(27)*STATES(10) -  CONSTS(28)*STATES(11)
      RATES(10) = ALGBRC(21) - ALGBRC(22)
      ALGBRC(24) = ( CONSTS(32)*STATES(15))/(CONSTS(31)+STATES(15))
      RATES(22) = ALGBRC(24) - ALGBRC(21)
      ALGBRC(26) =  CONSTS(35)*STATES(12) -  CONSTS(36)*STATES(5)*STATES(16)
      ALGBRC(25) =  CONSTS(33)*STATES(11)*STATES(21) -  CONSTS(34)*STATES(12)
      ALGBRC(27) =  CONSTS(37)*STATES(12)*STATES(23) -  CONSTS(38)*STATES(13)
      RATES(12) = ALGBRC(25) - (ALGBRC(26)+ALGBRC(27))
      ALGBRC(23) =  CONSTS(29)*STATES(11) -  CONSTS(30)*STATES(15)*STATES(5)
      ALGBRC(28) =  CONSTS(39)*STATES(13) -  CONSTS(40)*STATES(17)*STATES(5)
      RATES(5) = (ALGBRC(10)+ALGBRC(14)+ALGBRC(19)+ALGBRC(23)+ALGBRC(26)+ALGBRC(28)) - (ALGBRC(11)+ALGBRC(12)+ALGBRC(15)+ALGBRC(21))
      ALGBRC(20) =  CONSTS(23)*STATES(14) -  CONSTS(24)*STATES(21)*STATES(23)
      ALGBRC(29) =  CONSTS(41)*STATES(15)*STATES(21) -  CONSTS(42)*STATES(16)
      RATES(21) = ALGBRC(20) - (ALGBRC(15)+ALGBRC(25)+ALGBRC(29))
      ALGBRC(30) =  CONSTS(43)*STATES(16)*STATES(23) -  CONSTS(44)*STATES(17)
      RATES(16) = (ALGBRC(26)+ALGBRC(29)) - ALGBRC(30)
      RATES(23) = ALGBRC(20) - (ALGBRC(17)+ALGBRC(27)+ALGBRC(30))
      ALGBRC(31) =  CONSTS(45)*STATES(17) -  CONSTS(46)*STATES(15)*STATES(14)
      RATES(15) = (ALGBRC(23)+ALGBRC(31)) - (ALGBRC(29)+ALGBRC(24))
      RATES(17) = (ALGBRC(28)+ALGBRC(30)) - ALGBRC(31)
      ALGBRC(32) =  CONSTS(47)*STATES(11)*STATES(14) -  CONSTS(48)*STATES(13)
      RATES(11) = ALGBRC(22) - (ALGBRC(32)+ALGBRC(23)+ALGBRC(25))
      RATES(13) = (ALGBRC(27)+ALGBRC(32)) - ALGBRC(28)
      RATES(14) = (ALGBRC(19)+ALGBRC(31)) - (ALGBRC(20)+ALGBRC(32))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(2)*STATES(1) -  CONSTS(2)*STATES(3)
      ALGBRC(9) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      ALGBRC(10) =  CONSTS(5)*STATES(4) -  CONSTS(6)*STATES(5)
      ALGBRC(11) = ( CONSTS(8)*STATES(5))/(CONSTS(7)+STATES(5))
      ALGBRC(12) =  CONSTS(9)*STATES(5)*STATES(18) -  CONSTS(10)*STATES(6)
      ALGBRC(13) =  CONSTS(11)*STATES(6) -  CONSTS(12)*STATES(7)
      ALGBRC(14) =  CONSTS(13)*STATES(7) -  CONSTS(14)*STATES(5)*STATES(19)
      ALGBRC(16) = ( CONSTS(16)*STATES(19))/(CONSTS(15)+STATES(19))
      ALGBRC(15) =  CONSTS(17)*STATES(5)*STATES(21) -  CONSTS(18)*STATES(8)
      ALGBRC(17) =  CONSTS(19)*STATES(8)*STATES(23) -  CONSTS(20)*STATES(9)
      ALGBRC(18) =  CONSTS(49)*STATES(19) -  CONSTS(50)*STATES(20)
      ALGBRC(19) =  CONSTS(21)*STATES(9) -  CONSTS(22)*STATES(5)*STATES(14)
      ALGBRC(21) =  CONSTS(25)*STATES(5)*STATES(22) -  CONSTS(26)*STATES(10)
      ALGBRC(22) =  CONSTS(27)*STATES(10) -  CONSTS(28)*STATES(11)
      ALGBRC(24) = ( CONSTS(32)*STATES(15))/(CONSTS(31)+STATES(15))
      ALGBRC(26) =  CONSTS(35)*STATES(12) -  CONSTS(36)*STATES(5)*STATES(16)
      ALGBRC(25) =  CONSTS(33)*STATES(11)*STATES(21) -  CONSTS(34)*STATES(12)
      ALGBRC(27) =  CONSTS(37)*STATES(12)*STATES(23) -  CONSTS(38)*STATES(13)
      ALGBRC(23) =  CONSTS(29)*STATES(11) -  CONSTS(30)*STATES(15)*STATES(5)
      ALGBRC(28) =  CONSTS(39)*STATES(13) -  CONSTS(40)*STATES(17)*STATES(5)
      ALGBRC(20) =  CONSTS(23)*STATES(14) -  CONSTS(24)*STATES(21)*STATES(23)
      ALGBRC(29) =  CONSTS(41)*STATES(15)*STATES(21) -  CONSTS(42)*STATES(16)
      ALGBRC(30) =  CONSTS(43)*STATES(16)*STATES(23) -  CONSTS(44)*STATES(17)
      ALGBRC(31) =  CONSTS(45)*STATES(17) -  CONSTS(46)*STATES(15)*STATES(14)
      ALGBRC(32) =  CONSTS(47)*STATES(11)*STATES(14) -  CONSTS(48)*STATES(13)
      ALGBRC(2) =  2.00000*(STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13))
      ALGBRC(3) = STATES(7)+STATES(19)
      ALGBRC(4) = STATES(8)+STATES(9)+STATES(12)+STATES(13)
      ALGBRC(5) = STATES(12)+STATES(16)+STATES(13)+STATES(17)
      ALGBRC(6) = STATES(11)+STATES(12)+STATES(13)+STATES(15)+STATES(16)+STATES(17)
      ALGBRC(7) = STATES(11)+STATES(12)+STATES(13)
      ALGBRC(8) = STATES(9)+STATES(13)
      RETURN
      END