Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 66 entries in each of the rate and state variable arrays.
C There are a total of 255 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is AP1 in component AP1 (micromolar).
C CONSTS(1) is v_sap1 in component model_parameters (flux).
C CONSTS(2) is GF in component model_parameters (micromolar).
C CONSTS(3) is K_agf in component model_parameters (micromolar).
C CONSTS(4) is k_dap1 in component model_parameters (first_order_rate_constant).
C CONSTS(5) is eps in component model_parameters (dimensionless).
C STATES(2) is pRB in component pRB (micromolar).
C CONSTS(6) is v_sprb in component model_parameters (flux).
C CONSTS(7) is k_pc1 in component model_parameters (second_order_rate_constant).
C STATES(3) is E2F in component E2F (micromolar).
C CONSTS(8) is k_pc2 in component model_parameters (first_order_rate_constant).
C STATES(4) is pRBc1 in component pRBc1 (micromolar).
C CONSTS(9) is V_1 in component model_parameters (first_order_rate_constant).
C CONSTS(10) is K_1 in component model_parameters (micromolar).
C STATES(5) is Md in component Md (micromolar).
C STATES(6) is Mdp27 in component Mdp27 (micromolar).
C CONSTS(11) is V_2 in component model_parameters (flux).
C STATES(7) is pRBp in component pRBp (micromolar).
C CONSTS(12) is K_2 in component model_parameters (micromolar).
C CONSTS(13) is k_dprb in component model_parameters (first_order_rate_constant).
C CONSTS(14) is V_3 in component model_parameters (first_order_rate_constant).
C CONSTS(15) is K_3 in component model_parameters (micromolar).
C STATES(8) is Me in component Me (micromolar).
C CONSTS(16) is V_4 in component model_parameters (flux).
C STATES(9) is pRBpp in component pRBpp (micromolar).
C CONSTS(17) is K_4 in component model_parameters (micromolar).
C CONSTS(18) is k_pc3 in component model_parameters (second_order_rate_constant).
C CONSTS(19) is k_pc4 in component model_parameters (first_order_rate_constant).
C STATES(10) is pRBc2 in component pRBc2 (micromolar).
C CONSTS(20) is k_dpRBp in component model_parameters (first_order_rate_constant).
C CONSTS(21) is k_dpRBpp in component model_parameters (first_order_rate_constant).
C CONSTS(22) is v_se2f in component model_parameters (flux).
C CONSTS(23) is V_1e2f in component model_parameters (first_order_rate_constant).
C STATES(11) is Ma in component Ma (micromolar).
C CONSTS(24) is K_1e2f in component model_parameters (micromolar).
C CONSTS(25) is V_2e2f in component model_parameters (flux).
C STATES(12) is E2Fp in component E2Fp (micromolar).
C CONSTS(26) is K_2e2f in component model_parameters (micromolar).
C CONSTS(27) is k_de2f in component model_parameters (first_order_rate_constant).
C CONSTS(28) is k_de2fp in component model_parameters (first_order_rate_constant).
C STATES(13) is Cd in component Cd (micromolar).
C CONSTS(29) is k_cd1 in component model_parameters (first_order_rate_constant).
C CONSTS(30) is k_cd2 in component model_parameters (first_order_rate_constant).
C CONSTS(31) is K_i7 in component model_parameters (micromolar).
C CONSTS(32) is K_i8 in component model_parameters (micromolar).
C CONSTS(33) is k_com1 in component model_parameters (second_order_rate_constant).
C CONSTS(34) is Cdk4_tot in component model_parameters (micromolar).
C STATES(14) is Mdi in component Mdi (micromolar).
C CONSTS(35) is k_decom1 in component model_parameters (first_order_rate_constant).
C CONSTS(36) is V_dd in component model_parameters (flux).
C CONSTS(37) is K_dd in component model_parameters (micromolar).
C CONSTS(38) is k_ddd in component model_parameters (first_order_rate_constant).
C CONSTS(39) is V_m2d in component model_parameters (flux).
C CONSTS(40) is K_2d in component model_parameters (micromolar).
C CONSTS(41) is V_m1d in component model_parameters (flux).
C CONSTS(42) is K_1d in component model_parameters (micromolar).
C CONSTS(43) is k_c1 in component model_parameters (second_order_rate_constant).
C STATES(15) is p27 in component p27 (micromolar).
C CONSTS(44) is k_c2 in component model_parameters (first_order_rate_constant).
C STATES(16) is Ce in component Ce (micromolar).
C CONSTS(45) is k_ce in component model_parameters (first_order_rate_constant).
C CONSTS(46) is K_i9 in component model_parameters (micromolar).
C CONSTS(47) is K_i10 in component model_parameters (micromolar).
C CONSTS(48) is k_com2 in component model_parameters (second_order_rate_constant).
C CONSTS(49) is Cdk2_tot in component model_parameters (micromolar).
C STATES(17) is Mei in component Mei (micromolar).
C STATES(18) is Mep27 in component Mep27 (micromolar).
C STATES(19) is Mai in component Mai (micromolar).
C STATES(20) is Map27 in component Map27 (micromolar).
C CONSTS(50) is k_decom2 in component model_parameters (first_order_rate_constant).
C CONSTS(51) is V_de in component model_parameters (flux).
C STATES(21) is Skp2 in component Skp2 (micromolar).
C CONSTS(52) is K_dceskp2 in component model_parameters (micromolar).
C CONSTS(53) is K_de in component model_parameters (micromolar).
C CONSTS(54) is k_dde in component model_parameters (first_order_rate_constant).
C CONSTS(55) is V_m2e in component model_parameters (first_order_rate_constant).
C STATES(22) is Wee1 in component Wee1 (micromolar).
C CONSTS(56) is i_b1 in component model_parameters (micromolar).
C CONSTS(57) is K_2e in component model_parameters (micromolar).
C CONSTS(58) is V_m1e in component model_parameters (first_order_rate_constant).
C STATES(23) is Pe in component Pe (micromolar).
C CONSTS(59) is K_1e in component model_parameters (micromolar).
C CONSTS(60) is k_c3 in component model_parameters (second_order_rate_constant).
C CONSTS(61) is k_c4 in component model_parameters (first_order_rate_constant).
C CONSTS(62) is v_sskp2 in component model_parameters (flux).
C CONSTS(63) is V_dskp2 in component model_parameters (flux).
C CONSTS(64) is K_dskp2 in component model_parameters (micromolar).
C STATES(24) is Cdh1a in component Cdh1a (micromolar).
C CONSTS(65) is K_cdh1 in component model_parameters (micromolar).
C CONSTS(66) is k_ddskp2 in component model_parameters (first_order_rate_constant).
C STATES(25) is Pei in component Pei (micromolar).
C CONSTS(67) is v_spei in component model_parameters (flux).
C CONSTS(68) is V_6e in component model_parameters (flux).
C CONSTS(69) is x_e1 in component model_parameters (dimensionless).
C CONSTS(70) is x_e2 in component model_parameters (per_micromolar).
C STATES(26) is Chk1 in component Chk1 (micromolar).
C CONSTS(71) is K_6e in component model_parameters (micromolar).
C CONSTS(72) is V_m5e in component model_parameters (first_order_rate_constant).
C CONSTS(73) is a_e in component model_parameters (micromolar).
C CONSTS(74) is K_5e in component model_parameters (micromolar).
C CONSTS(75) is k_dpei in component model_parameters (first_order_rate_constant).
C CONSTS(76) is k_dpe in component model_parameters (first_order_rate_constant).
C STATES(27) is Ca in component Ca (micromolar).
C CONSTS(77) is k_ca in component model_parameters (first_order_rate_constant).
C CONSTS(78) is K_i11 in component model_parameters (micromolar).
C CONSTS(79) is K_i12 in component model_parameters (micromolar).
C CONSTS(80) is k_com3 in component model_parameters (second_order_rate_constant).
C CONSTS(81) is k_decom3 in component model_parameters (first_order_rate_constant).
C CONSTS(82) is V_da in component model_parameters (flux).
C CONSTS(83) is K_da in component model_parameters (micromolar).
C STATES(28) is Cdc20a in component Cdc20a (micromolar).
C CONSTS(84) is K_acdc20 in component model_parameters (micromolar).
C CONSTS(85) is k_dda in component model_parameters (first_order_rate_constant).
C CONSTS(86) is V_m2a in component model_parameters (first_order_rate_constant).
C CONSTS(87) is i_b2 in component model_parameters (micromolar).
C CONSTS(88) is K_2a in component model_parameters (micromolar).
C CONSTS(89) is V_m1a in component model_parameters (first_order_rate_constant).
C STATES(29) is Pa in component Pa (micromolar).
C CONSTS(90) is K_1a in component model_parameters (micromolar).
C CONSTS(91) is k_c5 in component model_parameters (second_order_rate_constant).
C CONSTS(92) is k_c6 in component model_parameters (first_order_rate_constant).
C CONSTS(93) is v_s1p27 in component model_parameters (flux).
C CONSTS(94) is v_s2p27 in component model_parameters (first_order_rate_constant).
C CONSTS(95) is K_i13 in component model_parameters (micromolar).
C CONSTS(96) is K_i14 in component model_parameters (micromolar).
C CONSTS(97) is k_c7 in component model_parameters (second_order_rate_constant).
C STATES(30) is Mb in component Mb (micromolar).
C CONSTS(98) is k_c8 in component model_parameters (first_order_rate_constant).
C STATES(31) is Mbp27 in component Mbp27 (micromolar).
C CONSTS(99) is V_1p27 in component model_parameters (first_order_rate_constant).
C CONSTS(100) is K_1p27 in component model_parameters (micromolar).
C CONSTS(101) is V_2p27 in component model_parameters (flux).
C CONSTS(102) is K_2p27 in component model_parameters (micromolar).
C STATES(32) is p27p in component p27p (micromolar).
C CONSTS(103) is k_ddp27 in component model_parameters (first_order_rate_constant).
C CONSTS(104) is V_dp27p in component model_parameters (flux).
C CONSTS(105) is K_dp27skp2 in component model_parameters (micromolar).
C CONSTS(106) is K_dp27p in component model_parameters (micromolar).
C CONSTS(107) is k_ddp27p in component model_parameters (first_order_rate_constant).
C STATES(33) is Cdh1i in component Cdh1i (micromolar).
C CONSTS(108) is V_2cdh1 in component model_parameters (first_order_rate_constant).
C CONSTS(109) is K_2cdh1 in component model_parameters (micromolar).
C CONSTS(110) is V_1cdh1 in component model_parameters (flux).
C CONSTS(111) is K_1cdh1 in component model_parameters (micromolar).
C CONSTS(112) is k_dcdh1i in component model_parameters (first_order_rate_constant).
C CONSTS(113) is v_scdh1a in component model_parameters (flux).
C CONSTS(114) is k_dcdh1a in component model_parameters (first_order_rate_constant).
C STATES(34) is Pai in component Pai (micromolar).
C CONSTS(115) is v_spai in component model_parameters (flux).
C CONSTS(116) is V_6a in component model_parameters (flux).
C CONSTS(117) is x_a1 in component model_parameters (dimensionless).
C CONSTS(118) is x_a2 in component model_parameters (per_micromolar).
C CONSTS(119) is K_6a in component model_parameters (micromolar).
C CONSTS(120) is V_m5a in component model_parameters (first_order_rate_constant).
C CONSTS(121) is a_a in component model_parameters (micromolar).
C CONSTS(122) is K_5a in component model_parameters (micromolar).
C CONSTS(123) is k_dpai in component model_parameters (first_order_rate_constant).
C CONSTS(124) is k_dpa in component model_parameters (first_order_rate_constant).
C STATES(35) is Cb in component Cb (micromolar).
C CONSTS(125) is v_cb in component model_parameters (flux).
C CONSTS(126) is k_com4 in component model_parameters (second_order_rate_constant).
C CONSTS(127) is Cdk1_tot in component model_parameters (micromolar).
C STATES(36) is Mbi in component Mbi (micromolar).
C CONSTS(128) is k_decom4 in component model_parameters (first_order_rate_constant).
C CONSTS(129) is V_db in component model_parameters (flux).
C CONSTS(130) is K_db in component model_parameters (micromolar).
C CONSTS(131) is K_dbcdc20 in component model_parameters (micromolar).
C CONSTS(132) is K_dbcdh1 in component model_parameters (micromolar).
C CONSTS(133) is k_ddb in component model_parameters (first_order_rate_constant).
C CONSTS(134) is V_m2b in component model_parameters (first_order_rate_constant).
C CONSTS(135) is i_b3 in component model_parameters (micromolar).
C CONSTS(136) is K_2b in component model_parameters (micromolar).
C CONSTS(137) is V_m1b in component model_parameters (first_order_rate_constant).
C STATES(37) is Pb in component Pb (micromolar).
C CONSTS(138) is K_1b in component model_parameters (micromolar).
C STATES(38) is Cdc20i in component Cdc20i (micromolar).
C CONSTS(139) is v_scdc20i in component model_parameters (flux).
C CONSTS(140) is V_m3b in component model_parameters (first_order_rate_constant).
C CONSTS(141) is K_3b in component model_parameters (micromolar).
C CONSTS(142) is V_m4b in component model_parameters (flux).
C CONSTS(143) is K_4b in component model_parameters (micromolar).
C CONSTS(144) is k_dcdc20i in component model_parameters (first_order_rate_constant).
C CONSTS(145) is k_dcdc20a in component model_parameters (first_order_rate_constant).
C STATES(39) is Pbi in component Pbi (micromolar).
C CONSTS(146) is v_spbi in component model_parameters (flux).
C CONSTS(147) is V_6b in component model_parameters (flux).
C CONSTS(148) is x_b1 in component model_parameters (dimensionless).
C CONSTS(149) is x_b2 in component model_parameters (per_micromolar).
C CONSTS(150) is K_6b in component model_parameters (micromolar).
C CONSTS(151) is V_m5b in component model_parameters (first_order_rate_constant).
C CONSTS(152) is a_b in component model_parameters (micromolar).
C CONSTS(153) is K_5b in component model_parameters (micromolar).
C CONSTS(154) is k_dpbi in component model_parameters (first_order_rate_constant).
C CONSTS(155) is k_dpb in component model_parameters (first_order_rate_constant).
C CONSTS(156) is v_swee1 in component model_parameters (flux).
C CONSTS(157) is k_sw in component model_parameters (first_order_rate_constant).
C STATES(40) is Mw in component Mw (micromolar).
C CONSTS(158) is V_m7b in component model_parameters (first_order_rate_constant).
C CONSTS(159) is i_b in component model_parameters (micromolar).
C CONSTS(160) is K_7b in component model_parameters (micromolar).
C CONSTS(161) is V_m8b in component model_parameters (flux).
C STATES(41) is Wee1p in component Wee1p (micromolar).
C CONSTS(162) is K_8b in component model_parameters (micromolar).
C CONSTS(163) is k_dwee1 in component model_parameters (first_order_rate_constant).
C CONSTS(164) is k_dwee1p in component model_parameters (first_order_rate_constant).
C STATES(42) is Cdc45 in component Cdc45 (micromolar).
C CONSTS(165) is V_1cdc45 in component model_parameters (first_order_rate_constant).
C CONSTS(166) is Cdc45_tot in component model_parameters (micromolar).
C CONSTS(167) is K_1cdc45 in component model_parameters (micromolar).
C CONSTS(168) is V_2cdc45 in component model_parameters (flux).
C CONSTS(169) is K_2cdc45 in component model_parameters (micromolar).
C CONSTS(170) is k_spol in component model_parameters (second_order_rate_constant).
C CONSTS(171) is Pol_tot in component model_parameters (micromolar).
C STATES(43) is Pol in component Pol (micromolar).
C CONSTS(172) is k_dpol in component model_parameters (first_order_rate_constant).
C STATES(44) is Primer in component Primer (micromolar).
C CONSTS(173) is k_sprim in component model_parameters (first_order_rate_constant).
C CONSTS(174) is k_dprim in component model_parameters (first_order_rate_constant).
C CONSTS(175) is k_aatr in component model_parameters (second_order_rate_constant).
C CONSTS(176) is ATR_tot in component model_parameters (micromolar).
C STATES(45) is ATR in component ATR (micromolar).
C CONSTS(177) is k_datr in component model_parameters (first_order_rate_constant).
C CONSTS(178) is V_1chk in component model_parameters (first_order_rate_constant).
C CONSTS(179) is Chk1_tot in component model_parameters (micromolar).
C CONSTS(180) is K_1chk in component model_parameters (micromolar).
C CONSTS(181) is V_2chk in component model_parameters (flux).
C CONSTS(182) is K_2chk in component model_parameters (micromolar).
C CONSTS(183) is v_sw in component model_parameters (flux).
C STATES(46) is BN in component BN (nanomolar).
C CONSTS(184) is n_gerard in component model_parameters (dimensionless).
C CONSTS(185) is K_iw in component model_parameters (nanomolar).
C CONSTS(186) is v_dw in component model_parameters (flux).
C CONSTS(187) is K_dw in component model_parameters (micromolar).
C STATES(47) is X in component X (micromolar).
C CONSTS(188) is V_1x in component model_parameters (first_order_rate_constant).
C CONSTS(189) is X_tot in component model_parameters (micromolar).
C CONSTS(190) is K_1x in component model_parameters (micromolar).
C CONSTS(191) is V_2x in component model_parameters (flux).
C CONSTS(192) is K_2x in component model_parameters (micromolar).
C STATES(48) is CbA in component CbA (micromolar).
C STATES(49) is MP in component MP (nanomolar).
C CONSTS(193) is vsP in component MP (nano_flux).
C CONSTS(194) is vmP in component MP (nano_flux).
C CONSTS(195) is kdmp in component MP (first_order_rate_constant_nano).
C CONSTS(196) is KAP in component MP (nanomolar).
C CONSTS(197) is KmP in component MP (nanomolar).
C CONSTS(198) is n in component model_parameters (dimensionless).
C STATES(50) is MC in component MC (nanomolar).
C CONSTS(199) is vsC in component MC (nano_flux).
C CONSTS(200) is vmC in component MC (nano_flux).
C CONSTS(201) is kdmc in component MC (first_order_rate_constant_nano).
C CONSTS(202) is KAC in component MC (nanomolar).
C CONSTS(203) is KmC in component MC (nanomolar).
C STATES(51) is MB in component MB (nanomolar).
C CONSTS(204) is vsB in component MB (nano_flux).
C CONSTS(205) is vmB in component MB (nano_flux).
C CONSTS(206) is kdmb in component MB (first_order_rate_constant_nano).
C CONSTS(207) is KIB in component MB (nanomolar).
C CONSTS(208) is KmB in component MB (nanomolar).
C CONSTS(209) is m in component model_parameters (dimensionless).
C STATES(52) is RN in component RN (nanomolar).
C STATES(53) is MR in component MR (nanomolar).
C CONSTS(210) is vsR in component MR (nano_flux).
C CONSTS(211) is vmR in component MR (nano_flux).
C CONSTS(212) is kdmr in component MR (first_order_rate_constant_nano).
C CONSTS(213) is KAR in component MR (nanomolar).
C CONSTS(214) is KmR in component MR (nanomolar).
C CONSTS(215) is h in component model_parameters (dimensionless).
C STATES(54) is PC in component PC (nanomolar).
C CONSTS(216) is ksP in component model_parameters (first_order_rate_constant_nano).
C CONSTS(217) is Kp in component model_parameters (nanomolar).
C CONSTS(218) is Kdp in component model_parameters (nanomolar).
C CONSTS(219) is k3 in component model_parameters (second_order_rate_constant_nano).
C CONSTS(220) is k4 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(221) is kdn in component model_parameters (first_order_rate_constant_nano).
C CONSTS(222) is V1P in component model_parameters (nano_flux).
C CONSTS(223) is V2P in component model_parameters (nano_flux).
C STATES(55) is PCP in component PCP (nanomolar).
C STATES(56) is PCC in component PCC (nanomolar).
C STATES(57) is CC in component CC (nanomolar).
C CONSTS(224) is ksC in component model_parameters (first_order_rate_constant_nano).
C CONSTS(225) is kdnc in component model_parameters (first_order_rate_constant_nano).
C CONSTS(226) is V1C in component model_parameters (nano_flux).
C CONSTS(227) is V2C in component model_parameters (nano_flux).
C STATES(58) is CCP in component CCP (nanomolar).
C STATES(59) is RC in component RC (nanomolar).
C CONSTS(228) is ksR in component model_parameters (first_order_rate_constant_nano).
C CONSTS(229) is Kd in component model_parameters (nanomolar).
C CONSTS(230) is k9 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(231) is k10 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(232) is vdRC in component model_parameters (nano_flux).
C CONSTS(233) is vdPC in component model_parameters (nano_flux).
C CONSTS(234) is vdCC in component model_parameters (nano_flux).
C CONSTS(235) is k1 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(236) is k2 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(237) is V1PC in component model_parameters (nano_flux).
C CONSTS(238) is V2PC in component model_parameters (nano_flux).
C STATES(60) is PCCP in component PCCP (nanomolar).
C STATES(61) is PCN in component PCN (nanomolar).
C CONSTS(239) is k7 in component model_parameters (second_order_rate_constant_nano).
C CONSTS(240) is k8 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(241) is V3PC in component model_parameters (nano_flux).
C CONSTS(242) is V4PC in component model_parameters (nano_flux).
C STATES(62) is PCNP in component PCNP (nanomolar).
C STATES(63) is IN in component IN (nanomolar).
C CONSTS(243) is vdRN in component model_parameters (nano_flux).
C CONSTS(244) is vdPCC in component model_parameters (nano_flux).
C CONSTS(245) is vdPCN in component model_parameters (nano_flux).
C STATES(64) is BC in component BC (nanomolar).
C CONSTS(246) is ksB in component model_parameters (first_order_rate_constant_nano).
C CONSTS(247) is k5 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(248) is k6 in component model_parameters (first_order_rate_constant_nano).
C CONSTS(249) is V1B in component model_parameters (nano_flux).
C CONSTS(250) is V2B in component model_parameters (nano_flux).
C STATES(65) is BCP in component BCP (nanomolar).
C CONSTS(251) is vdBC in component model_parameters (nano_flux).
C CONSTS(252) is V3B in component model_parameters (nano_flux).
C CONSTS(253) is V4B in component model_parameters (nano_flux).
C STATES(66) is BNP in component BNP (nanomolar).
C CONSTS(254) is vdBN in component model_parameters (nano_flux).
C CONSTS(255) is vdIN in component model_parameters (nano_flux).
C RATES(1) is d/dt AP1 in component AP1 (micromolar).
C RATES(2) is d/dt pRB in component pRB (micromolar).
C RATES(4) is d/dt pRBc1 in component pRBc1 (micromolar).
C RATES(7) is d/dt pRBp in component pRBp (micromolar).
C RATES(10) is d/dt pRBc2 in component pRBc2 (micromolar).
C RATES(9) is d/dt pRBpp in component pRBpp (micromolar).
C RATES(3) is d/dt E2F in component E2F (micromolar).
C RATES(12) is d/dt E2Fp in component E2Fp (micromolar).
C RATES(13) is d/dt Cd in component Cd (micromolar).
C RATES(14) is d/dt Mdi in component Mdi (micromolar).
C RATES(5) is d/dt Md in component Md (micromolar).
C RATES(6) is d/dt Mdp27 in component Mdp27 (micromolar).
C RATES(16) is d/dt Ce in component Ce (micromolar).
C RATES(17) is d/dt Mei in component Mei (micromolar).
C RATES(8) is d/dt Me in component Me (micromolar).
C RATES(21) is d/dt Skp2 in component Skp2 (micromolar).
C RATES(18) is d/dt Mep27 in component Mep27 (micromolar).
C RATES(25) is d/dt Pei in component Pei (micromolar).
C RATES(23) is d/dt Pe in component Pe (micromolar).
C RATES(27) is d/dt Ca in component Ca (micromolar).
C RATES(19) is d/dt Mai in component Mai (micromolar).
C RATES(11) is d/dt Ma in component Ma (micromolar).
C RATES(20) is d/dt Map27 in component Map27 (micromolar).
C RATES(15) is d/dt p27 in component p27 (micromolar).
C RATES(32) is d/dt p27p in component p27p (micromolar).
C RATES(33) is d/dt Cdh1i in component Cdh1i (micromolar).
C RATES(24) is d/dt Cdh1a in component Cdh1a (micromolar).
C RATES(34) is d/dt Pai in component Pai (micromolar).
C RATES(29) is d/dt Pa in component Pa (micromolar).
C RATES(35) is d/dt Cb in component Cb (micromolar).
C RATES(36) is d/dt Mbi in component Mbi (micromolar).
C RATES(30) is d/dt Mb in component Mb (micromolar).
C RATES(31) is d/dt Mbp27 in component Mbp27 (micromolar).
C RATES(38) is d/dt Cdc20i in component Cdc20i (micromolar).
C RATES(28) is d/dt Cdc20a in component Cdc20a (micromolar).
C RATES(39) is d/dt Pbi in component Pbi (micromolar).
C RATES(37) is d/dt Pb in component Pb (micromolar).
C RATES(22) is d/dt Wee1 in component Wee1 (micromolar).
C RATES(41) is d/dt Wee1p in component Wee1p (micromolar).
C RATES(42) is d/dt Cdc45 in component Cdc45 (micromolar).
C RATES(43) is d/dt Pol in component Pol (micromolar).
C RATES(44) is d/dt Primer in component Primer (micromolar).
C RATES(45) is d/dt ATR in component ATR (micromolar).
C RATES(26) is d/dt Chk1 in component Chk1 (micromolar).
C RATES(40) is d/dt Mw in component Mw (micromolar).
C RATES(47) is d/dt X in component X (micromolar).
C RATES(48) is d/dt CbA in component CbA (micromolar).
C RATES(49) is d/dt MP in component MP (nanomolar).
C RATES(50) is d/dt MC in component MC (nanomolar).
C RATES(51) is d/dt MB in component MB (nanomolar).
C RATES(53) is d/dt MR in component MR (nanomolar).
C RATES(54) is d/dt PC in component PC (nanomolar).
C RATES(57) is d/dt CC in component CC (nanomolar).
C RATES(59) is d/dt RC in component RC (nanomolar).
C RATES(55) is d/dt PCP in component PCP (nanomolar).
C RATES(58) is d/dt CCP in component CCP (nanomolar).
C RATES(56) is d/dt PCC in component PCC (nanomolar).
C RATES(61) is d/dt PCN in component PCN (nanomolar).
C RATES(52) is d/dt RN in component RN (nanomolar).
C RATES(60) is d/dt PCCP in component PCCP (nanomolar).
C RATES(62) is d/dt PCNP in component PCNP (nanomolar).
C RATES(64) is d/dt BC in component BC (nanomolar).
C RATES(65) is d/dt BCP in component BCP (nanomolar).
C RATES(46) is d/dt BN in component BN (nanomolar).
C RATES(66) is d/dt BNP in component BNP (nanomolar).
C RATES(63) is d/dt IN in component IN (nanomolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.01
CONSTS(1) = 1
CONSTS(2) = 1
CONSTS(3) = 0.1
CONSTS(4) = 0.15
CONSTS(5) = 17
STATES(2) = 1
CONSTS(6) = 0.8
CONSTS(7) = 0.05
STATES(3) = 0.01
CONSTS(8) = 0.5
STATES(4) = 0.1
CONSTS(9) = 2.2
CONSTS(10) = 0.1
STATES(5) = 0.01
STATES(6) = 0.01
CONSTS(11) = 2
STATES(7) = 0.25
CONSTS(12) = 0.1
CONSTS(13) = 0.01
CONSTS(14) = 1
CONSTS(15) = 0.1
STATES(8) = 0.01
CONSTS(16) = 2
STATES(9) = 0.1
CONSTS(17) = 0.1
CONSTS(18) = 0.025
CONSTS(19) = 0.5
STATES(10) = 0.05
CONSTS(20) = 0.06
CONSTS(21) = 0.04
CONSTS(22) = 0.17
CONSTS(23) = 4
STATES(11) = 0.01
CONSTS(24) = 5
CONSTS(25) = 0.75
STATES(12) = 0.01
CONSTS(26) = 5
CONSTS(27) = 0.002
CONSTS(28) = 1.1
STATES(13) = 0.01
CONSTS(29) = 0.4
CONSTS(30) = 0.005
CONSTS(31) = 0.1
CONSTS(32) = 2
CONSTS(33) = 0.175
CONSTS(34) = 1.5
STATES(14) = 0.01
CONSTS(35) = 0.1
CONSTS(36) = 5
CONSTS(37) = 0.1
CONSTS(38) = 0.005
CONSTS(39) = 0.2
CONSTS(40) = 0.1
CONSTS(41) = 1
CONSTS(42) = 0.1
CONSTS(43) = 0.15
STATES(15) = 0.01
CONSTS(44) = 0.05
STATES(16) = 0.01
CONSTS(45) = 0.29
CONSTS(46) = 0.1
CONSTS(47) = 2
CONSTS(48) = 0.2
CONSTS(49) = 2
STATES(17) = 0.01
STATES(18) = 0.01
STATES(19) = 0.01
STATES(20) = 0.01
CONSTS(50) = 0.1
CONSTS(51) = 3
STATES(21) = 0.01
CONSTS(52) = 2
CONSTS(53) = 0.1
CONSTS(54) = 0.005
CONSTS(55) = 1.4
STATES(22) = 0.1
CONSTS(56) = 0.5
CONSTS(57) = 0.1
CONSTS(58) = 2
STATES(23) = 0.01
CONSTS(59) = 0.1
CONSTS(60) = 0.2
CONSTS(61) = 0.1
CONSTS(62) = 0.15
CONSTS(63) = 1.1
CONSTS(64) = 0.5
STATES(24) = 0.01
CONSTS(65) = 0.4
CONSTS(66) = 0.005
STATES(25) = 0.01
CONSTS(67) = 0.13
CONSTS(68) = 0.8
CONSTS(69) = 1
CONSTS(70) = 1
STATES(26) = 0.01
CONSTS(71) = 0.1
CONSTS(72) = 5
CONSTS(73) = 0.25
CONSTS(74) = 0.1
CONSTS(75) = 0.15
CONSTS(76) = 0.075
STATES(27) = 0.01
CONSTS(77) = 0.0375
CONSTS(78) = 0.1
CONSTS(79) = 2
CONSTS(80) = 0.2
CONSTS(81) = 0.1
CONSTS(82) = 2.5
CONSTS(83) = 1.1
STATES(28) = 0.01
CONSTS(84) = 2
CONSTS(85) = 0.005
CONSTS(86) = 1.85
CONSTS(87) = 0.5
CONSTS(88) = 0.1
CONSTS(89) = 2
STATES(29) = 0.01
CONSTS(90) = 0.1
CONSTS(91) = 0.15
CONSTS(92) = 0.125
CONSTS(93) = 0.8
CONSTS(94) = 0.1
CONSTS(95) = 0.1
CONSTS(96) = 2
CONSTS(97) = 0.12
STATES(30) = 0.01
CONSTS(98) = 0.2
STATES(31) = 0.01
CONSTS(99) = 100
CONSTS(100) = 0.5
CONSTS(101) = 0.1
CONSTS(102) = 0.5
STATES(32) = 0.01
CONSTS(103) = 0.06
CONSTS(104) = 5
CONSTS(105) = 0.1
CONSTS(106) = 0.1
CONSTS(107) = 0.01
STATES(33) = 0.01
CONSTS(108) = 8
CONSTS(109) = 0.01
CONSTS(110) = 1.25
CONSTS(111) = 0.01
CONSTS(112) = 0.2
CONSTS(113) = 0.11
CONSTS(114) = 0.1
STATES(34) = 0.01
CONSTS(115) = 0.105
CONSTS(116) = 1
CONSTS(117) = 1
CONSTS(118) = 1
CONSTS(119) = 0.1
CONSTS(120) = 4
CONSTS(121) = 0.2
CONSTS(122) = 0.1
CONSTS(123) = 0.15
CONSTS(124) = 0.075
STATES(35) = 0.01
CONSTS(125) = 0.05
CONSTS(126) = 0.25
CONSTS(127) = 0.5
STATES(36) = 0.01
CONSTS(128) = 0.1
CONSTS(129) = 0.06
CONSTS(130) = 0.005
CONSTS(131) = 0.2
CONSTS(132) = 0.1
CONSTS(133) = 0.005
CONSTS(134) = 2.1
CONSTS(135) = 0.5
CONSTS(136) = 0.1
CONSTS(137) = 3.9
STATES(37) = 0.01
CONSTS(138) = 0.1
STATES(38) = 0.01
CONSTS(139) = 0.1
CONSTS(140) = 8
CONSTS(141) = 0.1
CONSTS(142) = 0.7
CONSTS(143) = 0.1
CONSTS(144) = 0.14
CONSTS(145) = 0.05
STATES(39) = 0.01
CONSTS(146) = 0.12
CONSTS(147) = 1
CONSTS(148) = 1
CONSTS(149) = 1
CONSTS(150) = 0.1
CONSTS(151) = 5
CONSTS(152) = 0.11
CONSTS(153) = 0.1
CONSTS(154) = 0.2
CONSTS(155) = 0.1
CONSTS(156) = 0.06
CONSTS(157) = 5
STATES(40) = 0
CONSTS(158) = 1.2
CONSTS(159) = 0.75
CONSTS(160) = 0.1
CONSTS(161) = 1
STATES(41) = 0.01
CONSTS(162) = 0.1
CONSTS(163) = 0.1
CONSTS(164) = 0.2
STATES(42) = 0.01
CONSTS(165) = 0.8
CONSTS(166) = 0.5
CONSTS(167) = 0.02
CONSTS(168) = 0.12
CONSTS(169) = 0.02
CONSTS(170) = 0.8
CONSTS(171) = 0.5
STATES(43) = 0.01
CONSTS(172) = 0.2
STATES(44) = 0.01
CONSTS(173) = 0.05
CONSTS(174) = 0.15
CONSTS(175) = 0.022
CONSTS(176) = 0.5
STATES(45) = 0.01
CONSTS(177) = 0.15
CONSTS(178) = 4
CONSTS(179) = 0.5
CONSTS(180) = 0.5
CONSTS(181) = 0.1
CONSTS(182) = 0.5
CONSTS(183) = 0
STATES(46) = 0.1
CONSTS(184) = 4
CONSTS(185) = 0.5
CONSTS(186) = 0.12
CONSTS(187) = 0.5
STATES(47) = 0.01
CONSTS(188) = 10
CONSTS(189) = 1
CONSTS(190) = 0.1
CONSTS(191) = 2
CONSTS(192) = 0.1
STATES(48) = 0.01
STATES(49) = 0.1
CONSTS(193) = 2.4
CONSTS(194) = 2.2
CONSTS(195) = 0.02
CONSTS(196) = 0.6
CONSTS(197) = 0.3
CONSTS(198) = 2
STATES(50) = 1.2
CONSTS(199) = 2.2
CONSTS(200) = 2
CONSTS(201) = 0.02
CONSTS(202) = 0.6
CONSTS(203) = 0.4
STATES(51) = 9
CONSTS(204) = 1.8
CONSTS(205) = 1.3
CONSTS(206) = 0.02
CONSTS(207) = 2.2
CONSTS(208) = 0.4
CONSTS(209) = 2
STATES(52) = 0.1
STATES(53) = 1.5
CONSTS(210) = 1.6
CONSTS(211) = 1.6
CONSTS(212) = 0.02
CONSTS(213) = 0.6
CONSTS(214) = 0.4
CONSTS(215) = 2
STATES(54) = 0.1
CONSTS(216) = 1.2
CONSTS(217) = 1.006
CONSTS(218) = 0.1
CONSTS(219) = 0.8
CONSTS(220) = 0.4
CONSTS(221) = 0.02
CONSTS(222) = 9.6
CONSTS(223) = 0.6
STATES(55) = 0.1
STATES(56) = 0.1
STATES(57) = 0.1
CONSTS(224) = 3.2
CONSTS(225) = 0.02
CONSTS(226) = 1.2
CONSTS(227) = 0.2
STATES(58) = 0.1
STATES(59) = 0.1
CONSTS(228) = 1.7
CONSTS(229) = 0.3
CONSTS(230) = 0.8
CONSTS(231) = 0.4
CONSTS(232) = 4.4
CONSTS(233) = 3.4
CONSTS(234) = 1.4
CONSTS(235) = 0.8
CONSTS(236) = 0.4
CONSTS(237) = 2.4
CONSTS(238) = 0.2
STATES(60) = 0.1
STATES(61) = 0.1
CONSTS(239) = 1
CONSTS(240) = 0.2
CONSTS(241) = 2.4
CONSTS(242) = 0.2
STATES(62) = 0.1
STATES(63) = 0.1
CONSTS(243) = 0.8
CONSTS(244) = 1.4
CONSTS(245) = 1.4
STATES(64) = 0.1
CONSTS(246) = 0.32
CONSTS(247) = 0.8
CONSTS(248) = 0.4
CONSTS(249) = 1.4
CONSTS(250) = 0.2
STATES(65) = 0.1
CONSTS(251) = 3
CONSTS(252) = 1.4
CONSTS(253) = 0.4
STATES(66) = 0.1
CONSTS(254) = 3
CONSTS(255) = 1.6
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = (( CONSTS(1)*CONSTS(2))/(CONSTS(3)+CONSTS(2)) - CONSTS(4)*STATES(1))*CONSTS(5)
RATES(2) = (((((CONSTS(6) - CONSTS(7)*STATES(2)*STATES(3))+ CONSTS(8)*STATES(4)) - (( CONSTS(9)*STATES(2))/(CONSTS(10)+STATES(2)))*(STATES(5)+STATES(6)))+( CONSTS(11)*STATES(7))/(CONSTS(12)+STATES(7))) - CONSTS(13)*STATES(2))*CONSTS(5)
RATES(4) = ( CONSTS(7)*STATES(2)*STATES(3) - CONSTS(8)*STATES(4))*CONSTS(5)
RATES(7) = (((((( (( CONSTS(9)*STATES(2))/(CONSTS(10)+STATES(2)))*(STATES(5)+STATES(6)) - ( CONSTS(11)*STATES(7))/(CONSTS(12)+STATES(7))) - (( CONSTS(14)*STATES(7))/(CONSTS(15)+STATES(7)))*STATES(8))+( CONSTS(16)*STATES(9))/(CONSTS(17)+STATES(9))) - CONSTS(18)*STATES(7)*STATES(3))+ CONSTS(19)*STATES(10)) - CONSTS(20)*STATES(7))*CONSTS(5)
RATES(10) = ( CONSTS(18)*STATES(7)*STATES(3) - CONSTS(19)*STATES(10))*CONSTS(5)
RATES(9) = (( (( CONSTS(14)*STATES(7))/(CONSTS(15)+STATES(7)))*STATES(8) - ( CONSTS(16)*STATES(9))/(CONSTS(17)+STATES(9))) - CONSTS(21)*STATES(9))*CONSTS(5)
RATES(3) = (((((((CONSTS(22) - CONSTS(7)*STATES(2)*STATES(3))+ CONSTS(8)*STATES(4)) - CONSTS(18)*STATES(7)*STATES(3))+ CONSTS(19)*STATES(10)) - ( CONSTS(23)*STATES(11)*STATES(3))/(CONSTS(24)+STATES(3)))+( CONSTS(25)*STATES(12))/(CONSTS(26)+STATES(12))) - CONSTS(27)*STATES(3))*CONSTS(5)
RATES(12) = ((( CONSTS(23)*STATES(11)*STATES(3))/(CONSTS(24)+STATES(3)) - ( CONSTS(25)*STATES(12))/(CONSTS(26)+STATES(12))) - CONSTS(28)*STATES(12))*CONSTS(5)
RATES(13) = ((((( CONSTS(29)*STATES(1)+( (( CONSTS(30)*STATES(3)*CONSTS(31))/(CONSTS(31)+STATES(2)))*CONSTS(32))/(CONSTS(32)+STATES(7))) - CONSTS(33)*STATES(13)*(CONSTS(34) - (STATES(14)+STATES(5)+STATES(6))))+ CONSTS(35)*STATES(14)) - ( CONSTS(36)*STATES(13))/(CONSTS(37)+STATES(13))) - CONSTS(38)*STATES(13))*CONSTS(5)
RATES(14) = ((( CONSTS(33)*STATES(13)*(CONSTS(34) - (STATES(14)+STATES(5)+STATES(6))) - CONSTS(35)*STATES(14))+( CONSTS(39)*STATES(5))/(CONSTS(40)+STATES(5))) - ( CONSTS(41)*STATES(14))/(CONSTS(42)+STATES(14)))*CONSTS(5)
RATES(5) = (((( CONSTS(41)*STATES(14))/(CONSTS(42)+STATES(14)) - ( CONSTS(39)*STATES(5))/(CONSTS(40)+STATES(5))) - CONSTS(43)*STATES(5)*STATES(15))+ CONSTS(44)*STATES(6))*CONSTS(5)
RATES(6) = ( CONSTS(43)*STATES(5)*STATES(15) - CONSTS(44)*STATES(6))*CONSTS(5)
RATES(16) = ((((( (( CONSTS(45)*STATES(3)*CONSTS(46))/(CONSTS(46)+STATES(2)))*CONSTS(47))/(CONSTS(47)+STATES(7)) - CONSTS(48)*STATES(16)*(CONSTS(49) - (STATES(17)+STATES(8)+STATES(18)+STATES(19)+STATES(11)+STATES(20))))+ CONSTS(50)*STATES(17)) - ( (( CONSTS(51)*STATES(21))/(CONSTS(52)+STATES(21)))*STATES(16))/(CONSTS(53)+STATES(16))) - CONSTS(54)*STATES(16))*CONSTS(5)
RATES(17) = ((( CONSTS(48)*STATES(16)*(CONSTS(49) - (STATES(17)+STATES(8)+STATES(18)+STATES(19)+STATES(11)+STATES(20))) - CONSTS(50)*STATES(17))+( CONSTS(55)*(STATES(22)+CONSTS(56))*STATES(8))/(CONSTS(57)+STATES(8))) - ( CONSTS(58)*STATES(23)*STATES(17))/(CONSTS(59)+STATES(17)))*CONSTS(5)
RATES(8) = (((( CONSTS(58)*STATES(23)*STATES(17))/(CONSTS(59)+STATES(17)) - ( CONSTS(55)*(STATES(22)+CONSTS(56))*STATES(8))/(CONSTS(57)+STATES(8))) - CONSTS(60)*STATES(8)*STATES(15))+ CONSTS(61)*STATES(18))*CONSTS(5)
RATES(21) = ((CONSTS(62) - ( (( CONSTS(63)*STATES(21))/(CONSTS(64)+STATES(21)))*STATES(24))/(CONSTS(65)+STATES(24))) - CONSTS(66)*STATES(21))*CONSTS(5)
RATES(18) = ( CONSTS(60)*STATES(8)*STATES(15) - CONSTS(61)*STATES(18))*CONSTS(5)
RATES(25) = (((CONSTS(67)+( CONSTS(68)*(CONSTS(69)+ CONSTS(70)*STATES(26))*STATES(23))/(CONSTS(71)+STATES(23))) - ( CONSTS(72)*(STATES(8)+CONSTS(73))*STATES(25))/(CONSTS(74)+STATES(25))) - CONSTS(75)*STATES(25))*CONSTS(5)
RATES(23) = (( CONSTS(72)*(STATES(8)+CONSTS(73))*(STATES(25)/(CONSTS(74)+STATES(25))) - CONSTS(68)*(CONSTS(69)+ CONSTS(70)*STATES(26))*(STATES(23)/(CONSTS(71)+STATES(23)))) - CONSTS(76)*STATES(23))*CONSTS(5)
RATES(27) = ((((( (( CONSTS(77)*STATES(3)*CONSTS(78))/(CONSTS(78)+STATES(2)))*CONSTS(79))/(CONSTS(79)+STATES(7)) - CONSTS(80)*STATES(27)*(CONSTS(49) - (STATES(17)+STATES(8)+STATES(18)+STATES(19)+STATES(11)+STATES(20))))+ CONSTS(81)*STATES(19)) - ( (( CONSTS(82)*STATES(27))/(CONSTS(83)+STATES(27)))*STATES(28))/(CONSTS(84)+STATES(28))) - CONSTS(85)*STATES(27))*CONSTS(5)
RATES(19) = ((( CONSTS(80)*STATES(27)*(CONSTS(49) - (STATES(17)+STATES(8)+STATES(18)+STATES(19)+STATES(11)+STATES(20))) - CONSTS(81)*STATES(19))+( CONSTS(86)*(STATES(22)+CONSTS(87))*STATES(11))/(CONSTS(88)+STATES(11))) - ( CONSTS(89)*STATES(29)*STATES(19))/(CONSTS(90)+STATES(19)))*CONSTS(5)
RATES(11) = (((( CONSTS(89)*STATES(29)*STATES(19))/(CONSTS(90)+STATES(19)) - ( CONSTS(86)*(STATES(22)+CONSTS(87))*STATES(11))/(CONSTS(88)+STATES(11))) - CONSTS(91)*STATES(11)*STATES(15))+ CONSTS(92)*STATES(20))*CONSTS(5)
RATES(20) = ( CONSTS(91)*STATES(11)*STATES(15) - CONSTS(92)*STATES(20))*CONSTS(5)
RATES(15) = ((((((((((((CONSTS(93)+( (( CONSTS(94)*STATES(3)*CONSTS(95))/(CONSTS(95)+STATES(2)))*CONSTS(96))/(CONSTS(96)+STATES(7))) - CONSTS(43)*STATES(5)*STATES(15))+ CONSTS(44)*STATES(6)) - CONSTS(60)*STATES(8)*STATES(15))+ CONSTS(61)*STATES(18)) - CONSTS(91)*STATES(11)*STATES(15))+ CONSTS(92)*STATES(20)) - CONSTS(97)*STATES(30)*STATES(15))+ CONSTS(98)*STATES(31)) - ( CONSTS(99)*STATES(8)*STATES(15))/(CONSTS(100)+STATES(15)))+( CONSTS(101)*STATES(32))/(CONSTS(102)+STATES(32))) - CONSTS(103)*STATES(15))*CONSTS(5)
RATES(32) = (((( CONSTS(99)*STATES(8)*STATES(15))/(CONSTS(100)+STATES(15)) - ( CONSTS(101)*STATES(32))/(CONSTS(102)+STATES(32))) - ( (( CONSTS(104)*STATES(21))/(CONSTS(105)+STATES(21)))*STATES(32))/(CONSTS(106)+STATES(32))) - CONSTS(107)*STATES(32))*CONSTS(5)
RATES(33) = (( (( CONSTS(108)*STATES(24))/(CONSTS(109)+STATES(24)))*(STATES(11)+STATES(30)) - ( CONSTS(110)*STATES(33))/(CONSTS(111)+STATES(33))) - CONSTS(112)*STATES(33))*CONSTS(5)
RATES(24) = (((CONSTS(113)+( CONSTS(110)*STATES(33))/(CONSTS(111)+STATES(33))) - (( CONSTS(108)*STATES(24))/(CONSTS(109)+STATES(24)))*(STATES(11)+STATES(30))) - CONSTS(114)*STATES(24))*CONSTS(5)
RATES(34) = (((CONSTS(115)+( CONSTS(116)*(CONSTS(117)+ CONSTS(118)*STATES(26))*STATES(29))/(CONSTS(119)+STATES(29))) - ( CONSTS(120)*(STATES(11)+CONSTS(121))*STATES(34))/(CONSTS(122)+STATES(34))) - CONSTS(123)*STATES(34))*CONSTS(5)
RATES(29) = ((( CONSTS(120)*(STATES(11)+CONSTS(121))*STATES(34))/(CONSTS(122)+STATES(34)) - ( CONSTS(116)*(CONSTS(117)+ CONSTS(118)*STATES(26))*STATES(29))/(CONSTS(119)+STATES(29))) - CONSTS(124)*STATES(29))*CONSTS(5)
RATES(35) = ((((CONSTS(125) - CONSTS(126)*STATES(35)*(CONSTS(127) - (STATES(36)+STATES(30)+STATES(31))))+ CONSTS(128)*STATES(36)) - (( CONSTS(129)*STATES(35))/(CONSTS(130)+STATES(35)))*(STATES(28)/(CONSTS(131)+STATES(28))+STATES(24)/(CONSTS(132)+STATES(24)))) - CONSTS(133)*STATES(35))*CONSTS(5)
RATES(36) = ((( CONSTS(126)*STATES(35)*(CONSTS(127) - (STATES(36)+STATES(30)+STATES(31))) - CONSTS(128)*STATES(36))+( CONSTS(134)*(STATES(22)+CONSTS(135))*STATES(30))/(CONSTS(136)+STATES(30))) - ( CONSTS(137)*STATES(37)*STATES(36))/(CONSTS(138)+STATES(36)))*CONSTS(5)
RATES(30) = (((( CONSTS(137)*STATES(37)*STATES(36))/(CONSTS(138)+STATES(36)) - ( CONSTS(134)*(STATES(22)+CONSTS(135))*STATES(30))/(CONSTS(136)+STATES(30))) - CONSTS(97)*STATES(30)*STATES(15))+ CONSTS(98)*STATES(31))*CONSTS(5)
RATES(31) = ( CONSTS(97)*STATES(30)*STATES(15) - CONSTS(98)*STATES(31))*CONSTS(5)
RATES(38) = (((CONSTS(139) - ( CONSTS(140)*STATES(30)*STATES(38))/(CONSTS(141)+STATES(38)))+( CONSTS(142)*STATES(28))/(CONSTS(143)+STATES(28))) - CONSTS(144)*STATES(38))*CONSTS(5)
RATES(28) = ((( CONSTS(140)*STATES(30)*STATES(38))/(CONSTS(141)+STATES(38)) - ( CONSTS(142)*STATES(28))/(CONSTS(143)+STATES(28))) - CONSTS(145)*STATES(28))*CONSTS(5)
RATES(39) = (((CONSTS(146)+( CONSTS(147)*(CONSTS(148)+ CONSTS(149)*STATES(26))*STATES(37))/(CONSTS(150)+STATES(37))) - ( CONSTS(151)*(STATES(30)+CONSTS(152))*STATES(39))/(CONSTS(153)+STATES(39))) - CONSTS(154)*STATES(39))*CONSTS(5)
RATES(37) = ((( CONSTS(151)*(STATES(30)+CONSTS(152))*STATES(39))/(CONSTS(153)+STATES(39)) - ( CONSTS(147)*(CONSTS(148)+ CONSTS(149)*STATES(26))*STATES(37))/(CONSTS(150)+STATES(37))) - CONSTS(155)*STATES(37))*CONSTS(5)
RATES(22) = ((((CONSTS(156)+ CONSTS(157)*STATES(40)) - ( CONSTS(158)*(STATES(30)+CONSTS(159))*STATES(22))/(CONSTS(160)+STATES(22)))+( CONSTS(161)*STATES(41))/(CONSTS(162)+STATES(41))) - CONSTS(163)*STATES(22))*CONSTS(5)
RATES(41) = ((( CONSTS(158)*(STATES(30)+CONSTS(159))*STATES(22))/(CONSTS(160)+STATES(22)) - ( CONSTS(161)*STATES(41))/(CONSTS(162)+STATES(41))) - CONSTS(164)*STATES(41))*CONSTS(5)
RATES(42) = (((( CONSTS(165)*STATES(8)*(CONSTS(166) - STATES(42)))/(CONSTS(167)+(CONSTS(166) - STATES(42))) - ( CONSTS(168)*STATES(42))/(CONSTS(169)+STATES(42))) - CONSTS(170)*(CONSTS(171) - STATES(43))*STATES(42))+ CONSTS(172)*STATES(43))*CONSTS(5)
RATES(43) = ( CONSTS(170)*(CONSTS(171) - STATES(43))*STATES(42) - CONSTS(172)*STATES(43))*CONSTS(5)
RATES(44) = ((( CONSTS(173)*STATES(43) - CONSTS(174)*STATES(44)) - CONSTS(175)*(CONSTS(176) - STATES(45))*STATES(44))+ CONSTS(177)*STATES(45))*CONSTS(5)
RATES(45) = ( CONSTS(175)*(CONSTS(176) - STATES(45))*STATES(44) - CONSTS(177)*STATES(45))*CONSTS(5)
RATES(26) = (( CONSTS(178)*STATES(45)*(CONSTS(179) - STATES(26)))/(CONSTS(180)+(CONSTS(179) - STATES(26))) - ( CONSTS(181)*STATES(26))/(CONSTS(182)+STATES(26)))*CONSTS(5)
RATES(40) = ( CONSTS(183)*STATES(46) ** CONSTS(184))/(CONSTS(185) ** CONSTS(184)+STATES(46) ** CONSTS(184)) - ( CONSTS(186)*STATES(40))/(CONSTS(187)+STATES(40))
RATES(47) = ( CONSTS(188)*STATES(11)*((CONSTS(189) - STATES(47))/(CONSTS(190)+(CONSTS(189) - STATES(47)))) - CONSTS(191)*(STATES(47)/(CONSTS(192)+STATES(47))))*CONSTS(5)
RATES(48) = (((( CONSTS(125)*STATES(47)*1.00000 - CONSTS(126)*STATES(48)*(CONSTS(127) - (STATES(36)+STATES(30)+STATES(31))))+ CONSTS(128)*STATES(36)) - CONSTS(129)*(STATES(48)/(CONSTS(130)+STATES(48)))*(STATES(28)/(CONSTS(131)+STATES(28))+STATES(24)/(CONSTS(132)+STATES(24)))) - CONSTS(133)*STATES(48))*CONSTS(5)
RATES(49) = ( CONSTS(193)*STATES(46) ** CONSTS(198))/(CONSTS(196) ** CONSTS(198)+STATES(46) ** CONSTS(198)) - (( CONSTS(194)*STATES(49))/(CONSTS(197)+STATES(49))+ CONSTS(195)*STATES(49))
RATES(50) = ( CONSTS(199)*STATES(46) ** CONSTS(198))/(CONSTS(202) ** CONSTS(198)+STATES(46) ** CONSTS(198)) - (( CONSTS(200)*STATES(50))/(CONSTS(203)+STATES(50))+ CONSTS(201)*STATES(50))
RATES(51) = ( CONSTS(204)*CONSTS(207) ** CONSTS(209))/(CONSTS(207) ** CONSTS(209)+STATES(52) ** CONSTS(209)) - (( CONSTS(205)*STATES(51))/(CONSTS(208)+STATES(51))+ CONSTS(206)*STATES(51))
RATES(53) = ( CONSTS(210)*STATES(46) ** CONSTS(215))/(CONSTS(213) ** CONSTS(215)+STATES(46) ** CONSTS(215)) - (( CONSTS(211)*STATES(53))/(CONSTS(214)+STATES(53))+ CONSTS(212)*STATES(53))
RATES(54) = ( CONSTS(216)*STATES(49)+( CONSTS(223)*STATES(55))/(CONSTS(218)+STATES(55))+ CONSTS(220)*STATES(56)) - (( CONSTS(222)*STATES(54))/(CONSTS(217)+STATES(54))+ CONSTS(219)*STATES(54)*STATES(57)+ CONSTS(221)*STATES(54))
RATES(57) = ( CONSTS(224)*STATES(50)+( CONSTS(227)*STATES(58))/(CONSTS(218)+STATES(58))+ CONSTS(220)*STATES(56)) - (( CONSTS(226)*STATES(57))/(CONSTS(217)+STATES(57))+ CONSTS(219)*STATES(54)*STATES(57)+ CONSTS(225)*STATES(57))
RATES(59) = ( CONSTS(228)*STATES(53)+ CONSTS(231)*STATES(52)) - ( CONSTS(230)*STATES(59)+( CONSTS(232)*STATES(59))/(CONSTS(229)+STATES(59))+ CONSTS(221)*STATES(59))
RATES(55) = ( CONSTS(222)*STATES(54))/(CONSTS(217)+STATES(54)) - (( CONSTS(223)*STATES(55))/(CONSTS(218)+STATES(55))+( CONSTS(233)*STATES(55))/(CONSTS(229)+STATES(55))+ CONSTS(221)*STATES(55))
RATES(58) = ( CONSTS(226)*STATES(57))/(CONSTS(217)+STATES(57)) - (( CONSTS(227)*STATES(58))/(CONSTS(218)+STATES(58))+( CONSTS(234)*STATES(58))/(CONSTS(229)+STATES(58))+ CONSTS(221)*STATES(58))
RATES(56) = (( CONSTS(238)*STATES(60))/(CONSTS(218)+STATES(60))+ CONSTS(219)*STATES(54)*STATES(57)+ CONSTS(236)*STATES(61)) - (( CONSTS(237)*STATES(56))/(CONSTS(217)+STATES(56))+ CONSTS(220)*STATES(56)+ CONSTS(235)*STATES(56)+ CONSTS(221)*STATES(56))
RATES(61) = (( CONSTS(242)*STATES(62))/(CONSTS(218)+STATES(62))+ CONSTS(235)*STATES(56)+ CONSTS(240)*STATES(63)) - (( CONSTS(241)*STATES(61))/(CONSTS(217)+STATES(61))+ CONSTS(236)*STATES(61)+ CONSTS(239)*STATES(46)*STATES(61)+ CONSTS(221)*STATES(61))
RATES(52) = CONSTS(230)*STATES(59) - ( CONSTS(231)*STATES(52)+( CONSTS(243)*STATES(52))/(CONSTS(229)+STATES(52))+ CONSTS(221)*STATES(52))
RATES(60) = ( CONSTS(237)*STATES(56))/(CONSTS(217)+STATES(56)) - (( CONSTS(238)*STATES(60))/(CONSTS(218)+STATES(60))+( CONSTS(244)*STATES(60))/(CONSTS(229)+STATES(60))+ CONSTS(221)*STATES(60))
RATES(62) = ( CONSTS(241)*STATES(61))/(CONSTS(217)+STATES(61)) - (( CONSTS(242)*STATES(62))/(CONSTS(218)+STATES(62))+( CONSTS(245)*STATES(62))/(CONSTS(229)+STATES(62))+ CONSTS(221)*STATES(62))
RATES(64) = (( CONSTS(250)*STATES(65))/(CONSTS(218)+STATES(65))+ CONSTS(248)*STATES(46)+ CONSTS(246)*STATES(51)) - (( CONSTS(249)*STATES(64))/(CONSTS(217)+STATES(64))+ CONSTS(247)*STATES(64)+ CONSTS(221)*STATES(64))
RATES(65) = ( CONSTS(249)*STATES(64))/(CONSTS(217)+STATES(64)) - (( CONSTS(250)*STATES(65))/(CONSTS(218)+STATES(65))+( CONSTS(251)*STATES(65))/(CONSTS(229)+STATES(65))+ CONSTS(221)*STATES(65))
RATES(46) = (( CONSTS(253)*STATES(66))/(CONSTS(218)+STATES(66))+ CONSTS(247)*STATES(64)+ CONSTS(240)*STATES(63)) - (( CONSTS(252)*STATES(46))/(CONSTS(217)+STATES(46))+ CONSTS(248)*STATES(46)+ CONSTS(239)*STATES(46)*STATES(61)+ CONSTS(221)*STATES(46))
RATES(66) = ( CONSTS(252)*STATES(46))/(CONSTS(217)+STATES(46)) - (( CONSTS(253)*STATES(66))/(CONSTS(218)+STATES(66))+( CONSTS(254)*STATES(66))/(CONSTS(229)+STATES(66))+ CONSTS(221)*STATES(66))
RATES(63) = CONSTS(239)*STATES(46)*STATES(61) - ( CONSTS(240)*STATES(63)+( CONSTS(255)*STATES(63))/(CONSTS(229)+STATES(63))+ CONSTS(221)*STATES(63))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END