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 77 entries in each of the rate and state variable arrays.
C There are a total of 200 entries in the constant variable array.
C
C
C VOI is t in component environment (s).
C CONSTS(1) is Ca in component Calcium (um).
C CONSTS(2) is GTP in component undefinedvariables (um).
C CONSTS(3) is PKA in component undefinedvariables (um).
C STATES(1) is PKC in component KEL_DAGIP3_Model_Modules (um).
C STATES(2) is GAP in component B_Model_Modules (um).
C STATES(3) is Raf in component B_Model_Modules (um).
C STATES(4) is GEF_star in component DGBM_Model_Modules (um).
C STATES(5) is CaM_GEF in component DGBM_Model_Modules (um).
C STATES(6) is Gbg_GEF in component DGBM_Model_Modules (um).
C STATES(7) is CaGqPLC in component DGBM_Model_Modules (um).
C STATES(8) is CaPLC in component DGBM_Model_Modules (um).
C STATES(9) is NgCaM in component DGBM_Model_Modules (um).
C STATES(10) is GEF in component DGBM_Model_Modules (um).
C STATES(11) is Ng in component DGBM_Model_Modules (um).
C STATES(12) is Gbg in component DGBM_Model_Modules (um).
C STATES(13) is Ga_GDP in component DGBM_Model_Modules (um).
C STATES(14) is mGluR in component DGBM_Model_Modules (um).
C STATES(15) is PLC in component DGBM_Model_Modules (um).
C STATES(16) is Ca2CaM in component DGBM_Model_Modules (um).
C STATES(17) is Ca3CaM in component DGBM_Model_Modules (um).
C STATES(18) is Ng_star in component DGBM_Model_Modules (um).
C STATES(19) is Gabg in component DGBM_Model_Modules (um).
C STATES(20) is R in component DGBM_Model_Modules (um).
C STATES(21) is Gabg_R in component DGBM_Model_Modules (um).
C STATES(22) is Glu_synapse in component DGBM_Model_Modules (um).
C STATES(23) is Gabg_GluR in component DGBM_Model_Modules (um).
C STATES(24) is Glu in component DGBM_Model_Modules (um).
C STATES(25) is Ga_GTP in component DGBM_Model_Modules (um).
C STATES(26) is GqPLC in component DGBM_Model_Modules (um).
C STATES(27) is GEF_inact in component DGBM_Model_Modules (um).
C STATES(28) is Ca4CAM in component DGBM_Model_Modules (um).
C STATES(29) is CaM in component DGBM_Model_Modules (um).
C CONSTS(4) is one in component DGBM_Model_Modules (dimensionless).
C CONSTS(5) is kfD1 in component DGBM_Model_Modules (per_s).
C CONSTS(6) is kbD1 in component DGBM_Model_Modules (per_s).
C CONSTS(7) is kfD2 in component DGBM_Model_Modules (per_s).
C CONSTS(8) is kbD2 in component DGBM_Model_Modules (per_um_s).
C CONSTS(9) is kfD3 in component DGBM_Model_Modules (per_s).
C CONSTS(10) is kbD3 in component DGBM_Model_Modules (per_um_s).
C CONSTS(11) is kfD4 in component DGBM_Model_Modules (per_um_s).
C CONSTS(12) is kbD4 in component DGBM_Model_Modules (per_s).
C CONSTS(13) is kfD5 in component DGBM_Model_Modules (per_s).
C CONSTS(14) is kbD5 in component DGBM_Model_Modules (per_um_s).
C CONSTS(15) is kfD6 in component DGBM_Model_Modules (per_s).
C CONSTS(16) is kfD7 in component DGBM_Model_Modules (per_s).
C CONSTS(17) is kfD8 in component DGBM_Model_Modules (per_um_s).
C CONSTS(18) is kfD9 in component DGBM_Model_Modules (per_um_s).
C CONSTS(19) is kfG1 in component DGBM_Model_Modules (per_um_s).
C CONSTS(20) is kbG1 in component DGBM_Model_Modules (per_s).
C CONSTS(21) is kfG2 in component DGBM_Model_Modules (per_um_s).
C CONSTS(22) is kbG2 in component DGBM_Model_Modules (per_s).
C CONSTS(23) is kfG3 in component DGBM_Model_Modules (per_um_s).
C CONSTS(24) is kbG3 in component DGBM_Model_Modules (per_s).
C CONSTS(25) is kfG4 in component DGBM_Model_Modules (per_um_s).
C CONSTS(26) is kbG4 in component DGBM_Model_Modules (per_s).
C CONSTS(27) is kfG5 in component DGBM_Model_Modules (per_s).
C CONSTS(28) is kfB2 in component DGBM_Model_Modules (per_s).
C CONSTS(29) is kbB2 in component DGBM_Model_Modules (per_s).
C CONSTS(30) is kmB3 in component DGBM_Model_Modules (um).
C CONSTS(31) is kfB5 in component DGBM_Model_Modules (per_um_s).
C CONSTS(32) is kbB5 in component DGBM_Model_Modules (per_s).
C CONSTS(33) is kfB7 in component DGBM_Model_Modules (per_s).
C CONSTS(34) is kbB7 in component DGBM_Model_Modules (per_s).
C CONSTS(35) is kfB8 in component DGBM_Model_Modules (per_um_s).
C CONSTS(36) is kbB8 in component DGBM_Model_Modules (per_s).
C CONSTS(37) is kmB1 in component DGBM_Model_Modules (um).
C CONSTS(38) is VmaxB1 in component DGBM_Model_Modules (per_s).
C CONSTS(39) is kmB6 in component DGBM_Model_Modules (um).
C CONSTS(40) is VmaxB6 in component DGBM_Model_Modules (per_s).
C CONSTS(41) is kfM1 in component DGBM_Model_Modules (per_um2_s).
C CONSTS(42) is kbM1 in component DGBM_Model_Modules (per_s).
C CONSTS(43) is kfM2 in component DGBM_Model_Modules (per_um_s).
C CONSTS(44) is kbM2 in component DGBM_Model_Modules (per_s).
C CONSTS(45) is kfM3 in component DGBM_Model_Modules (per_um_s).
C CONSTS(46) is kbM3 in component DGBM_Model_Modules (per_s).
C CONSTS(47) is kfM4 in component DGBM_Model_Modules (per_um_s).
C CONSTS(48) is kbM4 in component DGBM_Model_Modules (per_s).
C CONSTS(49) is kfM6 in component DGBM_Model_Modules (per_s).
C CONSTS(50) is kmM5 in component DGBM_Model_Modules (um).
C CONSTS(51) is VmaxM5 in component DGBM_Model_Modules (um_per_s).
C CONSTS(52) is kmM7 in component DGBM_Model_Modules (um).
C CONSTS(53) is VmaxM7 in component DGBM_Model_Modules (per_s).
C CONSTS(54) is kmM8 in component DGBM_Model_Modules (um).
C CONSTS(55) is VmaxM8 in component DGBM_Model_Modules (per_s).
C CONSTS(56) is kmH1 in component DGBM_Model_Modules (um).
C STATES(30) is SHCstar_SOS_GRB2 in component A_Model_Module_SHC (um).
C STATES(31) is MAPK_star in component H_Model_Module_MAPK (um).
C STATES(32) is SOS in component A_Model_Module_SHC (um).
C STATES(33) is MAPKK_star_star in component H_Model_Module_MAPKK (um).
C STATES(34) is MAPKK_star in component H_Model_Module_MAPKK (um).
C STATES(35) is Raf_star in component B_Model_Modules (um).
C STATES(36) is GTPRasRaf_star in component B_Model_Modules (um).
C STATES(37) is Raf_star_star in component B_Model_Modules (um).
C STATES(38) is GTP_RAS in component B_Model_Modules (um).
C STATES(39) is GDP_RAS in component B_Model_Modules (um).
C STATES(40) is GAPstar in component B_Model_Modules (um).
C CONSTS(57) is PP2A in component B_Model_Modules (um).
C CONSTS(58) is one in component B_Model_Modules (dimensionless).
C CONSTS(59) is kfB12 in component B_Model_Modules (per_s).
C CONSTS(60) is kbB12 in component B_Model_Modules (per_s).
C CONSTS(61) is kmB9 in component B_Model_Modules (um).
C CONSTS(62) is VmaxB9 in component B_Model_Modules (per_s).
C CONSTS(63) is kmB10 in component B_Model_Modules (um).
C CONSTS(64) is VmaxB10 in component B_Model_Modules (per_s).
C CONSTS(65) is kmB11 in component B_Model_Modules (um).
C CONSTS(66) is VmaxB11 in component B_Model_Modules (per_s).
C CONSTS(67) is kmB13 in component B_Model_Modules (um).
C CONSTS(68) is VmaxB13 in component B_Model_Modules (per_s).
C CONSTS(69) is kfB4 in component B_Model_Modules (per_s).
C CONSTS(70) is kbB4 in component B_Model_Modules (per_s).
C CONSTS(71) is kmB3 in component B_Model_Modules (um).
C CONSTS(72) is VmaxB3 in component B_Model_Modules (per_s).
C CONSTS(73) is kmA7 in component B_Model_Modules (um).
C CONSTS(74) is VmaxA7 in component B_Model_Modules (per_s).
C CONSTS(75) is kfH5 in component B_Model_Modules (per_um_s).
C CONSTS(76) is kbH5 in component B_Model_Modules (per_s).
C CONSTS(77) is kmH1 in component B_Model_Modules (um).
C CONSTS(78) is kmH2 in component B_Model_Modules (um).
C CONSTS(79) is kmH3 in component B_Model_Modules (um).
C CONSTS(80) is kmH4 in component B_Model_Modules (um).
C CONSTS(81) is VmaxH1 in component B_Model_Modules (per_s).
C CONSTS(82) is VmaxH2 in component B_Model_Modules (per_s).
C CONSTS(83) is VmaxH3 in component B_Model_Modules (per_s).
C CONSTS(84) is VmaxH4 in component B_Model_Modules (per_s).
C CONSTS(85) is kmM8 in component B_Model_Modules (um).
C CONSTS(86) is kmB6 in component B_Model_Modules (um).
C CONSTS(87) is kmM7 in component B_Model_Modules (um).
C CONSTS(88) is kmA9 in component B_Model_Modules (um).
C CONSTS(89) is kmH8 in component B_Model_Modules (um).
C CONSTS(90) is kmH9 in component B_Model_Modules (um).
C STATES(41) is MAPKK in component H_Model_Module_MAPKK (um).
C CONSTS(91) is PP2A in component H_Model_Module_MAPKK (um).
C CONSTS(92) is one in component H_Model_Module_MAPKK (dimensionless).
C CONSTS(93) is kmH6 in component H_Model_Module_MAPKK (um).
C CONSTS(94) is kmH7 in component H_Model_Module_MAPKK (um).
C CONSTS(95) is kmH8 in component H_Model_Module_MAPKK (um).
C CONSTS(96) is kmH9 in component H_Model_Module_MAPKK (um).
C CONSTS(97) is VmaxH6 in component H_Model_Module_MAPKK (per_s).
C CONSTS(98) is VmaxH7 in component H_Model_Module_MAPKK (per_s).
C CONSTS(99) is VmaxH8 in component H_Model_Module_MAPKK (per_s).
C CONSTS(100) is VmaxH9 in component H_Model_Module_MAPKK (per_s).
C CONSTS(101) is kmH3 in component H_Model_Module_MAPKK (um).
C CONSTS(102) is kmH4 in component H_Model_Module_MAPKK (um).
C STATES(42) is MAPK in component H_Model_Module_MAPK (um).
C STATES(43) is MAPK_tyr in component H_Model_Module_MAPK (um).
C CONSTS(103) is MKP1 in component H_Model_Module_MAPK (um).
C CONSTS(104) is one in component H_Model_Module_MAPK (dimensionless).
C CONSTS(105) is kmH10 in component H_Model_Module_MAPK (um).
C CONSTS(106) is kmH11 in component H_Model_Module_MAPK (um).
C CONSTS(107) is kmH12 in component H_Model_Module_MAPK (um).
C CONSTS(108) is kmH13 in component H_Model_Module_MAPK (um).
C CONSTS(109) is VmaxH10 in component H_Model_Module_MAPK (per_s).
C CONSTS(110) is VmaxH11 in component H_Model_Module_MAPK (per_s).
C CONSTS(111) is VmaxH12 in component H_Model_Module_MAPK (per_s).
C CONSTS(112) is VmaxH13 in component H_Model_Module_MAPK (per_s).
C STATES(44) is EGF_EGFR in component A_Model_Module_EGFR (um).
C STATES(45) is CaPLCg in component F_Model_Module_PLC (um).
C STATES(46) is SHC in component A_Model_Module_SHC (um).
C STATES(47) is SHCstar in component A_Model_Module_SHC (um).
C STATES(48) is SOSstar_GRB2 in component A_Model_Module_SHC (um).
C STATES(49) is SOS_GRB2 in component A_Model_Module_SHC (um).
C STATES(50) is GRB2 in component A_Model_Module_SHC (um).
C STATES(51) is SOSstar in component A_Model_Module_SHC (um).
C CONSTS(113) is one in component A_Model_Module_SHC (dimensionless).
C CONSTS(114) is kmA3 in component A_Model_Module_SHC (um).
C CONSTS(115) is kfA4 in component A_Model_Module_SHC (per_s).
C CONSTS(116) is kbA4 in component A_Model_Module_SHC (per_s).
C CONSTS(117) is kfA5 in component A_Model_Module_SHC (per_um_s).
C CONSTS(118) is kbA5 in component A_Model_Module_SHC (per_s).
C CONSTS(119) is kfA6 in component A_Model_Module_SHC (per_um_s).
C CONSTS(120) is kbA6 in component A_Model_Module_SHC (per_s).
C CONSTS(121) is kfA8 in component A_Model_Module_SHC (per_s).
C CONSTS(122) is kbA8 in component A_Model_Module_SHC (per_s).
C CONSTS(123) is kfA10 in component A_Model_Module_SHC (per_um_s).
C CONSTS(124) is kbA10 in component A_Model_Module_SHC (per_s).
C CONSTS(125) is kmA7 in component A_Model_Module_SHC (um).
C CONSTS(126) is kmA9 in component A_Model_Module_SHC (um).
C CONSTS(127) is VmaxA7 in component A_Model_Module_SHC (per_s).
C CONSTS(128) is VmaxA3 in component A_Model_Module_SHC (per_s).
C CONSTS(129) is VmaxA9 in component A_Model_Module_SHC (per_s).
C CONSTS(130) is kmH2 in component A_Model_Module_SHC (um).
C CONSTS(131) is kmF4 in component A_Model_Module_SHC (um).
C STATES(52) is EGF_EGFR_INTERNAL in component A_Model_Module_EGFR (um).
C STATES(53) is EGF in component A_Model_Module_EGFR (um).
C STATES(54) is EGFR in component A_Model_Module_EGFR (um).
C CONSTS(132) is kfA1 in component A_Model_Module_EGFR (per_um_s).
C CONSTS(133) is kbA1 in component A_Model_Module_EGFR (per_s).
C CONSTS(134) is kfA2 in component A_Model_Module_EGFR (per_s).
C CONSTS(135) is kbA2 in component A_Model_Module_EGFR (per_s).
C STATES(55) is CaPLCg_star in component F_Model_Module_PLC (um).
C STATES(56) is PLCg in component F_Model_Module_PLC (um).
C STATES(57) is PLCg_star in component F_Model_Module_PLC (um).
C CONSTS(136) is kfF1 in component F_Model_Module_PLC (per_um_s).
C CONSTS(137) is kbF1 in component F_Model_Module_PLC (per_s).
C CONSTS(138) is kfF3 in component F_Model_Module_PLC (per_s).
C CONSTS(139) is kfF5 in component F_Model_Module_PLC (per_um_s).
C CONSTS(140) is kbF5 in component F_Model_Module_PLC (per_s).
C CONSTS(141) is kmF4 in component F_Model_Module_PLC (um).
C CONSTS(142) is VmaxF4 in component F_Model_Module_PLC (per_s).
C CONSTS(143) is kmA3 in component F_Model_Module_PLC (um).
C CONSTS(144) is one in component F_Model_Module_PLC (dimensionless).
C STATES(58) is PKC_i in component KEL_DAGIP3_Model_Modules (um).
C STATES(59) is DAGPKC in component KEL_DAGIP3_Model_Modules (um).
C STATES(60) is AADAGPKC in component KEL_DAGIP3_Model_Modules (um).
C STATES(61) is DAGCaPKC in component KEL_DAGIP3_Model_Modules (um).
C STATES(62) is CaPKC in component KEL_DAGIP3_Model_Modules (um).
C STATES(63) is AA in component KEL_DAGIP3_Model_Modules (um).
C STATES(64) is APC in component KEL_DAGIP3_Model_Modules (um).
C STATES(65) is PIP2_star in component KEL_DAGIP3_Model_Modules (um).
C STATES(66) is PLA2_cyt in component KEL_DAGIP3_Model_Modules (um).
C STATES(67) is PLA2_star in component KEL_DAGIP3_Model_Modules (um).
C STATES(68) is CaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C STATES(69) is PIP2CaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C STATES(70) is PIP2PLA2 in component KEL_DAGIP3_Model_Modules (um).
C STATES(71) is DAGCaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C STATES(72) is CaPLA2_star in component KEL_DAGIP3_Model_Modules (um).
C STATES(73) is PIP2 in component KEL_DAGIP3_Model_Modules (um).
C STATES(74) is DAG in component KEL_DAGIP3_Model_Modules (um).
C STATES(75) is IP3 in component KEL_DAGIP3_Model_Modules (um).
C STATES(76) is Inositol in component KEL_DAGIP3_Model_Modules (um).
C STATES(77) is PC in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(145) is kfK1 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(146) is kbK1 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(147) is kfK2 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(148) is kbK2 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(149) is kfK3 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(150) is kbK3 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(151) is kfK4 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(152) is kbK4 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(153) is kfK5 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(154) is kbK5 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(155) is kfK6 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(156) is kbK6 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(157) is kfK7 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(158) is kbK7 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(159) is kfK8 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(160) is kbK8 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(161) is kfK9 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(162) is kbK9 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(163) is kfK10 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(164) is kbK10 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(165) is kfE1 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(166) is kbE1 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(167) is kfE3 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(168) is kbE3 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(169) is kfE5 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(170) is kbE5 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(171) is kfE7 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(172) is kbE7 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(173) is kfE10 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(174) is kbE10 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(175) is kfE11 in component KEL_DAGIP3_Model_Modules (per_um_s).
C CONSTS(176) is kbE11 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(177) is kmE9 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(178) is VmaxE9 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(179) is kfE13 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(180) is kbE13 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(181) is kmE2 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(182) is kmE4 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(183) is kmE6 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(184) is kmE8 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(185) is kmE12 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(186) is VmaxE2 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(187) is VmaxE4 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(188) is VmaxE6 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(189) is VmaxE8 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(190) is VmaxE12 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(191) is kmF2 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(192) is kmF6 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(193) is VmaxF2 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(194) is VmaxF6 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(195) is kfG8 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(196) is kfG9 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(197) is kmG6 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(198) is VmaxG6 in component KEL_DAGIP3_Model_Modules (per_s).
C CONSTS(199) is kmG7 in component KEL_DAGIP3_Model_Modules (um).
C CONSTS(200) is VmaxG7 in component KEL_DAGIP3_Model_Modules (per_s).
C RATES(13) is d/dt Ga_GDP in component DGBM_Model_Modules (um).
C RATES(14) is d/dt mGluR in component DGBM_Model_Modules (um).
C RATES(19) is d/dt Gabg in component DGBM_Model_Modules (um).
C RATES(20) is d/dt R in component DGBM_Model_Modules (um).
C RATES(12) is d/dt Gbg in component DGBM_Model_Modules (um).
C RATES(21) is d/dt Gabg_R in component DGBM_Model_Modules (um).
C RATES(22) is d/dt Glu_synapse in component DGBM_Model_Modules (um).
C RATES(23) is d/dt Gabg_GluR in component DGBM_Model_Modules (um).
C RATES(24) is d/dt Glu in component DGBM_Model_Modules (um).
C RATES(25) is d/dt Ga_GTP in component DGBM_Model_Modules (um).
C RATES(15) is d/dt PLC in component DGBM_Model_Modules (um).
C RATES(26) is d/dt GqPLC in component DGBM_Model_Modules (um).
C RATES(8) is d/dt CaPLC in component DGBM_Model_Modules (um).
C RATES(7) is d/dt CaGqPLC in component DGBM_Model_Modules (um).
C RATES(27) is d/dt GEF_inact in component DGBM_Model_Modules (um).
C RATES(10) is d/dt GEF in component DGBM_Model_Modules (um).
C RATES(5) is d/dt CaM_GEF in component DGBM_Model_Modules (um).
C RATES(28) is d/dt Ca4CAM in component DGBM_Model_Modules (um).
C RATES(6) is d/dt Gbg_GEF in component DGBM_Model_Modules (um).
C RATES(4) is d/dt GEF_star in component DGBM_Model_Modules (um).
C RATES(29) is d/dt CaM in component DGBM_Model_Modules (um).
C RATES(9) is d/dt NgCaM in component DGBM_Model_Modules (um).
C RATES(16) is d/dt Ca2CaM in component DGBM_Model_Modules (um).
C RATES(17) is d/dt Ca3CaM in component DGBM_Model_Modules (um).
C RATES(18) is d/dt Ng_star in component DGBM_Model_Modules (um).
C RATES(11) is d/dt Ng in component DGBM_Model_Modules (um).
C RATES(38) is d/dt GTP_RAS in component B_Model_Modules (um).
C RATES(39) is d/dt GDP_RAS in component B_Model_Modules (um).
C RATES(2) is d/dt GAP in component B_Model_Modules (um).
C RATES(40) is d/dt GAPstar in component B_Model_Modules (um).
C RATES(36) is d/dt GTPRasRaf_star in component B_Model_Modules (um).
C RATES(3) is d/dt Raf in component B_Model_Modules (um).
C RATES(35) is d/dt Raf_star in component B_Model_Modules (um).
C RATES(37) is d/dt Raf_star_star in component B_Model_Modules (um).
C RATES(41) is d/dt MAPKK in component H_Model_Module_MAPKK (um).
C RATES(34) is d/dt MAPKK_star in component H_Model_Module_MAPKK (um).
C RATES(33) is d/dt MAPKK_star_star in component H_Model_Module_MAPKK (um).
C RATES(42) is d/dt MAPK in component H_Model_Module_MAPK (um).
C RATES(43) is d/dt MAPK_tyr in component H_Model_Module_MAPK (um).
C RATES(31) is d/dt MAPK_star in component H_Model_Module_MAPK (um).
C RATES(47) is d/dt SHCstar in component A_Model_Module_SHC (um).
C RATES(46) is d/dt SHC in component A_Model_Module_SHC (um).
C RATES(30) is d/dt SHCstar_SOS_GRB2 in component A_Model_Module_SHC (um).
C RATES(49) is d/dt SOS_GRB2 in component A_Model_Module_SHC (um).
C RATES(32) is d/dt SOS in component A_Model_Module_SHC (um).
C RATES(51) is d/dt SOSstar in component A_Model_Module_SHC (um).
C RATES(48) is d/dt SOSstar_GRB2 in component A_Model_Module_SHC (um).
C RATES(50) is d/dt GRB2 in component A_Model_Module_SHC (um).
C RATES(53) is d/dt EGF in component A_Model_Module_EGFR (um).
C RATES(54) is d/dt EGFR in component A_Model_Module_EGFR (um).
C RATES(44) is d/dt EGF_EGFR in component A_Model_Module_EGFR (um).
C RATES(52) is d/dt EGF_EGFR_INTERNAL in component A_Model_Module_EGFR (um).
C RATES(56) is d/dt PLCg in component F_Model_Module_PLC (um).
C RATES(45) is d/dt CaPLCg in component F_Model_Module_PLC (um).
C RATES(55) is d/dt CaPLCg_star in component F_Model_Module_PLC (um).
C RATES(57) is d/dt PLCg_star in component F_Model_Module_PLC (um).
C RATES(58) is d/dt PKC_i in component KEL_DAGIP3_Model_Modules (um).
C RATES(59) is d/dt DAGPKC in component KEL_DAGIP3_Model_Modules (um).
C RATES(60) is d/dt AADAGPKC in component KEL_DAGIP3_Model_Modules (um).
C RATES(62) is d/dt CaPKC in component KEL_DAGIP3_Model_Modules (um).
C RATES(61) is d/dt DAGCaPKC in component KEL_DAGIP3_Model_Modules (um).
C RATES(1) is d/dt PKC in component KEL_DAGIP3_Model_Modules (um).
C RATES(63) is d/dt AA in component KEL_DAGIP3_Model_Modules (um).
C RATES(66) is d/dt PLA2_cyt in component KEL_DAGIP3_Model_Modules (um).
C RATES(70) is d/dt PIP2PLA2 in component KEL_DAGIP3_Model_Modules (um).
C RATES(65) is d/dt PIP2_star in component KEL_DAGIP3_Model_Modules (um).
C RATES(67) is d/dt PLA2_star in component KEL_DAGIP3_Model_Modules (um).
C RATES(72) is d/dt CaPLA2_star in component KEL_DAGIP3_Model_Modules (um).
C RATES(68) is d/dt CaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C RATES(69) is d/dt PIP2CaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C RATES(71) is d/dt DAGCaPLA2 in component KEL_DAGIP3_Model_Modules (um).
C RATES(64) is d/dt APC in component KEL_DAGIP3_Model_Modules (um).
C RATES(73) is d/dt PIP2 in component KEL_DAGIP3_Model_Modules (um).
C RATES(74) is d/dt DAG in component KEL_DAGIP3_Model_Modules (um).
C RATES(75) is d/dt IP3 in component KEL_DAGIP3_Model_Modules (um).
C RATES(77) is d/dt PC in component KEL_DAGIP3_Model_Modules (um).
C RATES(76) is d/dt Inositol in component KEL_DAGIP3_Model_Modules (um).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 1
CONSTS(2) = 0.5
CONSTS(3) = 0.1
STATES(1) = 0
STATES(2) = 0.002
STATES(3) = 0.2
STATES(4) = 0
STATES(5) = 0
STATES(6) = 0
STATES(7) = 0
STATES(8) = 0
STATES(9) = 0
STATES(10) = 0
STATES(11) = 10
STATES(12) = 0.1
STATES(13) = 1
STATES(14) = 0.3
STATES(15) = 0.8
STATES(16) = 0
STATES(17) = 0
STATES(18) = 0
STATES(19) = 0
STATES(20) = 0
STATES(21) = 0
STATES(22) = 0
STATES(23) = 0
STATES(24) = 0
STATES(25) = 0
STATES(26) = 0
STATES(27) = 0.1
STATES(28) = 0
STATES(29) = 20
CONSTS(4) = 1
CONSTS(5) = 500
CONSTS(6) = 1000
CONSTS(7) = 0.000028
CONSTS(8) = 0.1
CONSTS(9) = 0.000028
CONSTS(10) = 10
CONSTS(11) = 0.000001
CONSTS(12) = 1
CONSTS(13) = 0.00000001
CONSTS(14) = 0.0001
CONSTS(15) = 0.01
CONSTS(16) = 0.0133
CONSTS(17) = 0.00001
CONSTS(18) = 0.0001
CONSTS(19) = 0.000005
CONSTS(20) = 1
CONSTS(21) = 0.0000042
CONSTS(22) = 1
CONSTS(23) = 0.00005
CONSTS(24) = 1
CONSTS(25) = 0.0000042
CONSTS(26) = 1
CONSTS(27) = 0.0133
CONSTS(28) = 1
CONSTS(29) = 0
CONSTS(30) = 3.333333333333
CONSTS(31) = 0.0001
CONSTS(32) = 1
CONSTS(33) = 1
CONSTS(34) = 0
CONSTS(35) = 0.00001
CONSTS(36) = 1
CONSTS(37) = 7.5
CONSTS(38) = 9
CONSTS(39) = 3.33333333333
CONSTS(40) = 4
CONSTS(41) = 0.0000000002
CONSTS(42) = 72
CONSTS(43) = 0.000006
CONSTS(44) = 10
CONSTS(45) = 0.000000775
CONSTS(46) = 10
CONSTS(47) = 0.0000005
CONSTS(48) = 1
CONSTS(49) = 0.005
CONSTS(50) = 10.012
CONSTS(51) = 0.67
CONSTS(52) = 28.626667
CONSTS(53) = 0.58
CONSTS(54) = 28.595
CONSTS(55) = 0.35
CONSTS(56) = 66.666666667
STATES(30) = 0
STATES(31) = 0
STATES(32) = 0.1
STATES(33) = 0
STATES(34) = 0
STATES(35) = 0
STATES(36) = 0
STATES(37) = 0
STATES(38) = 0
STATES(39) = 0.2
STATES(40) = 0
CONSTS(57) = 0.224
CONSTS(58) = 1
CONSTS(59) = 0.0001
CONSTS(60) = 0
CONSTS(61) = 0.50505
CONSTS(62) = 0.02
CONSTS(63) = 0.50505
CONSTS(64) = 0.02
CONSTS(65) = 0.50505
CONSTS(66) = 0.02
CONSTS(67) = 1.0104
CONSTS(68) = 10
CONSTS(69) = 0.1
CONSTS(70) = 0
CONSTS(71) = 3.33333333333
CONSTS(72) = 4
CONSTS(73) = 0.50505
CONSTS(74) = 0.02
CONSTS(75) = 0.00004
CONSTS(76) = 0.5
CONSTS(77) = 66.666666667
CONSTS(78) = 25.64166667
CONSTS(79) = 15.6565
CONSTS(80) = 15.6565
CONSTS(81) = 4
CONSTS(82) = 10
CONSTS(83) = 6
CONSTS(84) = 6
CONSTS(85) = 28.595
CONSTS(86) = 3.33333333333
CONSTS(87) = 28.626667
CONSTS(88) = 2.564166667
CONSTS(89) = 15.6565
CONSTS(90) = 15.6565
STATES(41) = 0.18
CONSTS(91) = 0.224
CONSTS(92) = 1
CONSTS(93) = 0.159091667
CONSTS(94) = 0.159091667
CONSTS(95) = 15.6565
CONSTS(96) = 15.6565
CONSTS(97) = 0.105
CONSTS(98) = 0.105
CONSTS(99) = 6
CONSTS(100) = 6
CONSTS(101) = 15.6565
CONSTS(102) = 15.6565
STATES(42) = 0.36
STATES(43) = 0
CONSTS(103) = 0.032
CONSTS(104) = 1
CONSTS(105) = 0.046296667
CONSTS(106) = 0.046296667
CONSTS(107) = 0.066666667
CONSTS(108) = 0.066666667
CONSTS(109) = 0.15
CONSTS(110) = 0.15
CONSTS(111) = 1
CONSTS(112) = 1
STATES(44) = 0
STATES(45) = 0
STATES(46) = 0.5
STATES(47) = 0
STATES(48) = 0
STATES(49) = 0
STATES(50) = 1
STATES(51) = 0
CONSTS(113) = 1
CONSTS(114) = 0.83333333333
CONSTS(115) = 0.0016667
CONSTS(116) = 0
CONSTS(117) = 0.000000041667
CONSTS(118) = 0.0168
CONSTS(119) = 0.0000008333
CONSTS(120) = 0.1
CONSTS(121) = 0.001
CONSTS(122) = 0
CONSTS(123) = 0.000000041667
CONSTS(124) = 0.0168
CONSTS(125) = 0.40404
CONSTS(126) = 2.564166667
CONSTS(127) = 0.02
CONSTS(128) = 0.2
CONSTS(129) = 10
CONSTS(130) = 25.64166667
CONSTS(131) = 0.33333333333
STATES(52) = 0
STATES(53) = 166.67
STATES(54) = 0.16667
CONSTS(132) = 0.000007
CONSTS(133) = 0.25
CONSTS(134) = 0.002
CONSTS(135) = 0.00033
STATES(55) = 0
STATES(56) = 0.82
STATES(57) = 0
CONSTS(136) = 0.0003
CONSTS(137) = 10
CONSTS(138) = 0.05
CONSTS(139) = 0.00002
CONSTS(140) = 10
CONSTS(141) = 0.33333333333
CONSTS(142) = 0.2
CONSTS(143) = 0.83333333333
CONSTS(144) = 1
STATES(58) = 1
STATES(59) = 0
STATES(60) = 0
STATES(61) = 0
STATES(62) = 0
STATES(63) = 0
STATES(64) = 30
STATES(65) = 2.5
STATES(66) = 0.4
STATES(67) = 0
STATES(68) = 0
STATES(69) = 0
STATES(70) = 0
STATES(71) = 0
STATES(72) = 0
STATES(73) = 10
STATES(74) = 0
STATES(75) = 0
STATES(76) = 0
STATES(77) = 0
CONSTS(145) = 1
CONSTS(146) = 50
CONSTS(147) = 0.0000000002
CONSTS(148) = 0.1
CONSTS(149) = 1.2705
CONSTS(150) = 3.5026
CONSTS(151) = 0.000000002
CONSTS(152) = 0.1
CONSTS(153) = 1
CONSTS(154) = 0.1
CONSTS(155) = 2
CONSTS(156) = 0.2
CONSTS(157) = 0.000001
CONSTS(158) = 0.5
CONSTS(159) = 0.000000013333
CONSTS(160) = 8.6348
CONSTS(161) = 0.000000001
CONSTS(162) = 0.1
CONSTS(163) = 0.00000003
CONSTS(164) = 2
CONSTS(165) = 0.000000002
CONSTS(166) = 0.5
CONSTS(167) = 0.00000016667
CONSTS(168) = 0.1
CONSTS(169) = 0.00000002
CONSTS(170) = 0.1
CONSTS(171) = 0.000000005
CONSTS(172) = 4
CONSTS(173) = 0.17
CONSTS(174) = 0
CONSTS(175) = 0.00001
CONSTS(176) = 0.1
CONSTS(177) = 25.64166667
CONSTS(178) = 20
CONSTS(179) = 0.4
CONSTS(180) = 0
CONSTS(181) = 20
CONSTS(182) = 20
CONSTS(183) = 20
CONSTS(184) = 20
CONSTS(185) = 20
CONSTS(186) = 11.04
CONSTS(187) = 5.4
CONSTS(188) = 36
CONSTS(189) = 60
CONSTS(190) = 120
CONSTS(191) = 97
CONSTS(192) = 19.79166667
CONSTS(193) = 14
CONSTS(194) = 57
CONSTS(195) = 0.15
CONSTS(196) = 2.5
CONSTS(197) = 5
CONSTS(198) = 48
CONSTS(199) = 19.84166667
CONSTS(200) = 10
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(13) = ( STATES(25)*CONSTS(16) - STATES(13)*CONSTS(17)*STATES(12))+ STATES(7)*CONSTS(27)
RATES(14) = (( CONSTS(13)*STATES(23) - CONSTS(14)*STATES(19)*STATES(14)) - ( CONSTS(9)*STATES(14) - CONSTS(10)*STATES(24)*STATES(20)))+ CONSTS(15)*STATES(23)
RATES(19) = ((( STATES(13)*STATES(12)*CONSTS(17) - ( STATES(19)*STATES(20)*CONSTS(11) - STATES(21)*CONSTS(12)))+ CONSTS(13)*STATES(23)) - CONSTS(14)*STATES(19)*STATES(14)) - CONSTS(18)*STATES(19)*CONSTS(2)
RATES(20) = ( CONSTS(9)*STATES(14) - CONSTS(10)*STATES(24)*STATES(20)) - ( STATES(19)*STATES(20)*CONSTS(11) - STATES(21)*CONSTS(12))
RATES(12) = ( - STATES(13)*STATES(12)*CONSTS(17)+ CONSTS(15)*STATES(23)+ CONSTS(18)*STATES(19)*CONSTS(2)) - ( CONSTS(35)*STATES(10)*STATES(12) - CONSTS(36)*STATES(6))
RATES(21) = (( STATES(19)*STATES(20)*CONSTS(11) - STATES(21)*CONSTS(12))+ CONSTS(7)*STATES(23)) - STATES(21)*STATES(24)*CONSTS(8)
RATES(22) = - ( CONSTS(5)*STATES(22) - CONSTS(6)*STATES(24))
RATES(23) = (- ( CONSTS(7)*STATES(23) - STATES(21)*STATES(24)*CONSTS(8)) - ( CONSTS(13)*STATES(23) - CONSTS(14)*STATES(19)*STATES(14))) - CONSTS(15)*STATES(23)
RATES(24) = (((( CONSTS(5)*STATES(22) - CONSTS(6)*STATES(24))+ CONSTS(7)*STATES(23)) - STATES(21)*STATES(24)*CONSTS(8))+ CONSTS(9)*STATES(14)) - CONSTS(10)*STATES(24)*STATES(20)
RATES(25) = (( CONSTS(18)*STATES(19)*CONSTS(2) - STATES(25)*CONSTS(16)) - ( STATES(15)*STATES(25)*CONSTS(21) - STATES(26)*CONSTS(22))) - ( STATES(7)*CONSTS(26)+ STATES(8)*STATES(25)*CONSTS(25))
RATES(15) = - ( STATES(15)*CONSTS(1)*CONSTS(19) - STATES(8)*CONSTS(20)) - ( STATES(15)*STATES(25)*CONSTS(21) - STATES(26)*CONSTS(22))
RATES(26) = ( STATES(15)*STATES(25)*CONSTS(21) - STATES(26)*CONSTS(22)) - ( CONSTS(1)*STATES(26)*CONSTS(23) - STATES(7)*CONSTS(24))
RATES(8) = (( STATES(15)*CONSTS(1)*CONSTS(19) - STATES(8)*CONSTS(20))+ STATES(7)*CONSTS(27)) - ( - STATES(7)*CONSTS(26)+ STATES(8)*STATES(25)*CONSTS(25))
RATES(7) = ((( STATES(26)*CONSTS(1)*CONSTS(23) - STATES(7)*CONSTS(24))+ STATES(8)*STATES(25)*CONSTS(25)) - STATES(7)*CONSTS(26)) - STATES(7)*CONSTS(27)
RATES(27) = ( CONSTS(3)*CONSTS(38)*STATES(10))/(CONSTS(37)+STATES(10)) - CONSTS(28)*STATES(27)
RATES(10) = (((- (((( CONSTS(3)*CONSTS(38)*STATES(10))/(CONSTS(37)+STATES(10))+( STATES(1)*CONSTS(40)*STATES(10))/( CONSTS(39)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(11)/CONSTS(52)+STATES(3)/CONSTS(56)+STATES(9)/CONSTS(54))+STATES(10))) - CONSTS(28)*STATES(27)) - CONSTS(33)*STATES(4)) - CONSTS(31)*STATES(28)*STATES(10))+ CONSTS(32)*STATES(5)) - CONSTS(35)*STATES(10)*STATES(12))+ CONSTS(36)*STATES(6)
RATES(5) = - CONSTS(32)*STATES(5)+ CONSTS(31)*STATES(10)*STATES(28)
RATES(28) = (- ( - CONSTS(32)*STATES(5)+ CONSTS(31)*STATES(10)*STATES(28))+ CONSTS(1)*STATES(17)*CONSTS(45)) - CONSTS(46)*STATES(28)
RATES(6) = - ( - CONSTS(35)*STATES(10)*STATES(12)+ CONSTS(36)*STATES(6))
RATES(4) = ( STATES(1)*CONSTS(40)*STATES(10))/( CONSTS(39)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(11)/CONSTS(52)+STATES(3)/CONSTS(56)+STATES(9)/CONSTS(54))+STATES(10)) - CONSTS(33)*STATES(4)
RATES(29) = (- ( CONSTS(41)*CONSTS(1)*CONSTS(1)*STATES(29) - CONSTS(42)*STATES(16)) - ( CONSTS(47)*STATES(11)*STATES(29) - CONSTS(48)*STATES(9)))+( STATES(9)*STATES(1)*CONSTS(55))/(STATES(9)+CONSTS(54))
RATES(9) = ( CONSTS(47)*STATES(11)*STATES(29) - CONSTS(48)*STATES(9)) - ( STATES(9)*STATES(1)*CONSTS(55))/(STATES(9)+ CONSTS(54)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(10)/CONSTS(39)+STATES(3)/CONSTS(56)+STATES(11)/CONSTS(52)))
RATES(16) = ( CONSTS(41)*CONSTS(1)*CONSTS(1)*STATES(29) - CONSTS(42)*STATES(16)) - ( CONSTS(1)*STATES(16)*CONSTS(43) - CONSTS(44)*STATES(17))
RATES(17) = ( CONSTS(1)*STATES(16)*CONSTS(43) - CONSTS(44)*STATES(17)) - ( CONSTS(1)*STATES(17)*CONSTS(45) - CONSTS(46)*STATES(28))
RATES(18) = ((( STATES(9)*STATES(1)*CONSTS(55))/(STATES(9)+ CONSTS(54)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(10)/CONSTS(39)+STATES(3)/CONSTS(56)+STATES(11)/CONSTS(52))) - CONSTS(49)*STATES(18)) - ( CONSTS(51)*STATES(18))/(STATES(18)+CONSTS(50)))+( STATES(11)*STATES(1)*CONSTS(53))/(STATES(11)+ CONSTS(52)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(10)/CONSTS(39)+STATES(3)/CONSTS(56)+STATES(9)/CONSTS(54)))
RATES(11) = ( CONSTS(49)*STATES(18)+( CONSTS(51)*STATES(18))/(STATES(18)+CONSTS(50))) - ( STATES(11)*STATES(1)*CONSTS(53))/(STATES(11)+ CONSTS(52)*(CONSTS(4)+STATES(2)/CONSTS(30)+STATES(10)/CONSTS(39)+STATES(3)/CONSTS(56)+STATES(9)/CONSTS(54)))
RATES(38) = (((( - STATES(2)*CONSTS(68)*STATES(38))/(CONSTS(67)+STATES(38)) - CONSTS(59)*STATES(38))+( STATES(4)*CONSTS(64)*STATES(39))/(CONSTS(63)+STATES(39))+( STATES(6)*CONSTS(66)*STATES(39))/(CONSTS(65)+STATES(39))+( STATES(5)*CONSTS(62)*STATES(39))/(CONSTS(61)+STATES(39))+( STATES(30)*CONSTS(74)*STATES(39))/(CONSTS(73)+STATES(39))) - STATES(38)*STATES(35)*CONSTS(75))+ STATES(36)*CONSTS(76)
RATES(39) = - ((( - STATES(2)*CONSTS(68)*STATES(38))/(CONSTS(67)+STATES(38)) - CONSTS(59)*STATES(38))+( STATES(4)*CONSTS(64)*STATES(39))/(CONSTS(63)+STATES(39))+( STATES(6)*CONSTS(66)*STATES(39))/(CONSTS(65)+STATES(39))+( STATES(5)*CONSTS(62)*STATES(39))/(CONSTS(61)+STATES(39))) - ( STATES(30)*CONSTS(74)*STATES(39))/(CONSTS(73)+STATES(39))
RATES(2) = ( - STATES(1)*CONSTS(72)*STATES(2))/( CONSTS(71)*(CONSTS(58)+STATES(3)/CONSTS(77)+STATES(10)/CONSTS(86)+STATES(11)/CONSTS(87)+STATES(9)/CONSTS(85))+STATES(2))+ CONSTS(69)*STATES(40)
RATES(40) = ( STATES(1)*CONSTS(72)*STATES(2))/( CONSTS(71)*(CONSTS(58)+STATES(3)/CONSTS(77)+STATES(10)/CONSTS(86)+STATES(11)/CONSTS(87)+STATES(9)/CONSTS(85))+STATES(2)) - CONSTS(69)*STATES(40)
RATES(36) = STATES(38)*STATES(35)*CONSTS(75) - STATES(36)*CONSTS(76)
RATES(3) = ( - STATES(1)*STATES(3)*CONSTS(81))/( CONSTS(77)*(CONSTS(58)+STATES(2)/CONSTS(71)+STATES(10)/CONSTS(86)+STATES(11)/CONSTS(87)+STATES(9)/CONSTS(85))+STATES(3))+( STATES(35)*CONSTS(57)*CONSTS(83))/( CONSTS(79)*(CONSTS(58)+STATES(37)/CONSTS(80)+STATES(33)/CONSTS(90)+STATES(34)/CONSTS(89))+STATES(35))
RATES(35) = ((( STATES(1)*STATES(3)*CONSTS(81))/( CONSTS(77)*(CONSTS(58)+STATES(2)/CONSTS(71)+STATES(10)/CONSTS(86)+STATES(11)/CONSTS(87)+STATES(9)/CONSTS(85))+STATES(3)) - ( STATES(35)*CONSTS(57)*CONSTS(83))/( CONSTS(79)*(CONSTS(58)+STATES(37)/CONSTS(80)+STATES(33)/CONSTS(90)+STATES(34)/CONSTS(89))+STATES(35))) - ( STATES(35)*STATES(31)*CONSTS(82))/( CONSTS(78)*(CONSTS(58)+STATES(32)/CONSTS(88))+STATES(35)))+( STATES(37)*CONSTS(57)*CONSTS(84))/( CONSTS(80)*(CONSTS(58)+STATES(35)/CONSTS(79)+STATES(33)/CONSTS(90)+STATES(34)/CONSTS(89))+STATES(37))
RATES(37) = ( STATES(35)*STATES(31)*CONSTS(82))/( CONSTS(78)*(CONSTS(58)+STATES(32)/CONSTS(88))+STATES(35)) - ( STATES(37)*CONSTS(57)*CONSTS(84))/( CONSTS(80)*(CONSTS(58)+STATES(35)/CONSTS(79)+STATES(33)/CONSTS(90)+STATES(34)/CONSTS(89))+STATES(37))
RATES(41) = ( - STATES(41)*CONSTS(97)*STATES(36))/( CONSTS(93)*(CONSTS(92)+STATES(34)/CONSTS(94))+STATES(41))+( STATES(34)*CONSTS(99)*CONSTS(91))/( CONSTS(95)*(CONSTS(92)+STATES(35)/CONSTS(101)+STATES(37)/CONSTS(102)+STATES(33)/CONSTS(96))+STATES(34))
RATES(34) = ((( STATES(41)*CONSTS(97)*STATES(36))/( CONSTS(93)*(CONSTS(92)+STATES(34)/CONSTS(94))+STATES(41)) - ( STATES(34)*CONSTS(99)*CONSTS(91))/( CONSTS(95)*(CONSTS(92)+STATES(35)/CONSTS(101)+STATES(37)/CONSTS(102)+STATES(33)/CONSTS(96))+STATES(34))) - ( STATES(34)*CONSTS(98)*STATES(36))/( CONSTS(94)*(CONSTS(92)+STATES(41)/CONSTS(93))+STATES(34)))+( STATES(33)*CONSTS(100)*CONSTS(91))/( CONSTS(96)*(CONSTS(92)+STATES(35)/CONSTS(101)+STATES(37)/CONSTS(102)+STATES(34)/CONSTS(95))+STATES(33))
RATES(33) = ( STATES(34)*CONSTS(98)*STATES(36))/( CONSTS(94)*(CONSTS(92)+STATES(41)/CONSTS(93))+STATES(34)) - ( STATES(33)*CONSTS(100)*CONSTS(91))/( CONSTS(96)*(CONSTS(92)+STATES(35)/CONSTS(101)+STATES(37)/CONSTS(102)+STATES(34)/CONSTS(95))+STATES(33))
RATES(42) = ( - STATES(42)*CONSTS(109)*STATES(33))/( CONSTS(105)*(CONSTS(104)+STATES(43)/CONSTS(106))+STATES(42))+( STATES(43)*CONSTS(111)*CONSTS(103))/( CONSTS(107)*(CONSTS(104)+STATES(31)/CONSTS(108))+STATES(43))
RATES(43) = ((( STATES(42)*CONSTS(109)*STATES(33))/(( CONSTS(105)*STATES(43))/CONSTS(106)+STATES(42)) - ( STATES(43)*CONSTS(111)*CONSTS(103))/( CONSTS(107)*(CONSTS(104)+STATES(31)/CONSTS(108))+STATES(43))) - ( STATES(43)*CONSTS(110)*STATES(33))/( CONSTS(106)*(CONSTS(104)+STATES(43)/CONSTS(107))+STATES(43)))+( STATES(31)*CONSTS(112)*CONSTS(103))/( CONSTS(108)*(CONSTS(104)+STATES(43)/CONSTS(107))+STATES(31))
RATES(31) = ( STATES(43)*CONSTS(110)*STATES(33))/( CONSTS(106)*(CONSTS(104)+STATES(42)/CONSTS(105))+STATES(43)) - ( STATES(31)*CONSTS(112)*CONSTS(103))/( CONSTS(108)*(CONSTS(104)+STATES(43)/CONSTS(107))+STATES(31))
RATES(47) = ((( STATES(44)*CONSTS(128)*STATES(46))/( CONSTS(114)*(CONSTS(113)+STATES(45)/CONSTS(131))+STATES(46)) - CONSTS(119)*STATES(47)*STATES(49))+ CONSTS(120)*STATES(30)) - CONSTS(115)*STATES(47)
RATES(46) = ( CONSTS(115)*STATES(47) - CONSTS(116)*STATES(46)) - ( STATES(44)*CONSTS(128)*STATES(46))/( CONSTS(114)*(CONSTS(113)+STATES(45)/CONSTS(131))+STATES(46))
RATES(30) = CONSTS(119)*STATES(47)*STATES(49) - CONSTS(120)*STATES(30)
RATES(49) = (( - CONSTS(119)*STATES(47)*STATES(49)+ CONSTS(120)*STATES(30)) - CONSTS(118)*STATES(49))+ CONSTS(117)*STATES(32)*STATES(50)
RATES(32) = ((( CONSTS(121)*STATES(51) - CONSTS(122)*STATES(32)) - ( STATES(31)*CONSTS(129)*STATES(32))/( CONSTS(126)*(CONSTS(113)+STATES(35)/CONSTS(130))+STATES(32)))+ CONSTS(118)*STATES(49)) - CONSTS(117)*STATES(32)*STATES(50)
RATES(51) = (( - CONSTS(121)*STATES(51)+ CONSTS(122)*STATES(32)+( STATES(31)*CONSTS(129)*STATES(32))/( CONSTS(126)*(CONSTS(113)+STATES(35)/CONSTS(130))+STATES(32))) - CONSTS(123)*STATES(51)*STATES(50))+ CONSTS(124)*STATES(48)
RATES(48) = CONSTS(123)*STATES(51)*STATES(50) - CONSTS(124)*STATES(48)
RATES(50) = (( CONSTS(118)*STATES(49) - CONSTS(117)*STATES(32)*STATES(50)) - CONSTS(123)*STATES(51)*STATES(50))+ CONSTS(124)*STATES(48)
RATES(53) = - CONSTS(132)*STATES(53)*STATES(54)+ CONSTS(133)*STATES(44)
RATES(54) = - CONSTS(132)*STATES(53)*STATES(54)+ CONSTS(133)*STATES(44)
RATES(44) = (( CONSTS(132)*STATES(53)*STATES(54) - CONSTS(133)*STATES(44)) - CONSTS(134)*STATES(44))+ CONSTS(135)*STATES(52)
RATES(52) = CONSTS(134)*STATES(44) - CONSTS(135)*STATES(52)
RATES(56) = - CONSTS(1)*STATES(56)*CONSTS(136)+ STATES(45)*CONSTS(137)
RATES(45) = (- ( - CONSTS(1)*STATES(56)*CONSTS(136)+ STATES(45)*CONSTS(137))+ CONSTS(138)*STATES(55)) - ( CONSTS(142)*STATES(44)*STATES(45))/( CONSTS(141)*(CONSTS(144)+STATES(46)/CONSTS(143))+STATES(45))
RATES(55) = (- ( CONSTS(138)*STATES(55) - ( CONSTS(142)*STATES(44)*STATES(45))/(CONSTS(141)+STATES(45)))+ CONSTS(1)*STATES(57)*CONSTS(139)) - CONSTS(140)*STATES(55)
RATES(57) = - CONSTS(1)*STATES(57)*CONSTS(139)+ CONSTS(140)*STATES(55)
RATES(58) = (((((( - STATES(58)*CONSTS(145)+ STATES(1)*CONSTS(146)) - STATES(58)*STATES(63)*CONSTS(147))+ STATES(1)*CONSTS(148)) - STATES(58)*STATES(74)*CONSTS(161))+ STATES(59)*CONSTS(162)) - STATES(58)*CONSTS(1)*CONSTS(157))+ STATES(62)*CONSTS(158)
RATES(59) = (( STATES(58)*STATES(74)*CONSTS(161) - STATES(59)*CONSTS(162)) - STATES(59)*STATES(63)*CONSTS(163))+ STATES(60)*CONSTS(164)
RATES(60) = (( STATES(59)*STATES(63)*CONSTS(163) - STATES(60)*CONSTS(164)) - STATES(60)*CONSTS(155))+ STATES(1)*CONSTS(156)
RATES(62) = (((((( STATES(58)*CONSTS(1)*CONSTS(157) - STATES(62)*CONSTS(158)) - STATES(62)*CONSTS(149))+ STATES(1)*CONSTS(150)) - STATES(62)*STATES(63)*CONSTS(151))+ STATES(1)*CONSTS(152)) - STATES(74)*STATES(62)*CONSTS(159))+ STATES(61)*CONSTS(160)
RATES(61) = (( STATES(74)*STATES(62)*CONSTS(159) - STATES(61)*CONSTS(160)) - STATES(61)*CONSTS(153))+ STATES(1)*CONSTS(154)
RATES(1) = (((((((((( STATES(58)*CONSTS(145) - STATES(1)*CONSTS(146))+ STATES(58)*STATES(63)*CONSTS(147)) - STATES(1)*CONSTS(148))+ STATES(62)*CONSTS(149)) - STATES(1)*CONSTS(150))+ STATES(62)*STATES(63)*CONSTS(151)) - STATES(1)*CONSTS(152))+ STATES(61)*CONSTS(153)) - STATES(1)*CONSTS(154))+ STATES(60)*CONSTS(155)) - STATES(1)*CONSTS(156)
RATES(63) = (((((( - STATES(58)*STATES(63)*CONSTS(147)+ STATES(1)*CONSTS(148)) - STATES(62)*STATES(63)*CONSTS(151))+ STATES(1)*CONSTS(152)) - STATES(59)*STATES(63)*CONSTS(163))+ STATES(60)*CONSTS(164)) - STATES(63)*CONSTS(179))+( STATES(72)*CONSTS(190)*STATES(64))/(CONSTS(185)+STATES(64))+( STATES(68)*CONSTS(187)*STATES(64))/(CONSTS(182)+STATES(64))+( STATES(71)*CONSTS(189)*STATES(64))/(CONSTS(184)+STATES(64))+( STATES(69)*CONSTS(188)*STATES(64))/(CONSTS(183)+STATES(64))+( STATES(70)*CONSTS(186)*STATES(64))/(CONSTS(181)+STATES(64))
RATES(66) = (((( STATES(67)*CONSTS(173) - STATES(65)*STATES(66)*CONSTS(165))+ STATES(70)*CONSTS(166)) - CONSTS(1)*STATES(66)*CONSTS(167))+ STATES(68)*CONSTS(168)) - ( STATES(42)*CONSTS(178)*STATES(66))/(CONSTS(177)+STATES(66))
RATES(70) = STATES(65)*STATES(66)*CONSTS(165) - STATES(70)*CONSTS(166)
RATES(65) = (( - STATES(65)*STATES(66)*CONSTS(165)+ STATES(70)*CONSTS(166)) - STATES(65)*STATES(68)*CONSTS(169))+ STATES(69)*CONSTS(170)
RATES(67) = (( - STATES(67)*CONSTS(173)+( STATES(42)*CONSTS(178)*STATES(66))/(CONSTS(177)+STATES(66))) - CONSTS(1)*STATES(67)*CONSTS(175))+ STATES(72)*CONSTS(176)
RATES(72) = CONSTS(1)*STATES(67)*CONSTS(175) - STATES(72)*CONSTS(176)
RATES(68) = (((( CONSTS(1)*STATES(66)*CONSTS(167) - STATES(68)*CONSTS(168)) - STATES(74)*STATES(68)*CONSTS(171))+ STATES(71)*CONSTS(172)) - STATES(65)*STATES(68)*CONSTS(169))+ STATES(69)*CONSTS(170)
RATES(69) = STATES(65)*STATES(68)*CONSTS(169) - STATES(69)*CONSTS(170)
RATES(71) = STATES(74)*STATES(68)*CONSTS(171) - STATES(71)*CONSTS(172)
RATES(64) = STATES(63)*CONSTS(179) - (( STATES(72)*CONSTS(190)*STATES(64))/(CONSTS(185)+STATES(64))+( STATES(68)*CONSTS(187)*STATES(64))/(CONSTS(182)+STATES(64))+( STATES(71)*CONSTS(189)*STATES(64))/(CONSTS(184)+STATES(64))+( STATES(69)*CONSTS(188)*STATES(64))/(CONSTS(183)+STATES(64))+( STATES(70)*CONSTS(186)*STATES(64))/(CONSTS(181)+STATES(64)))
RATES(73) = ((( - STATES(45)*CONSTS(193)*STATES(73))/(CONSTS(191)+STATES(73)) - ( STATES(55)*CONSTS(194)*STATES(73))/(CONSTS(192)+STATES(73))) - ( STATES(8)*CONSTS(200)*STATES(73))/(CONSTS(199)+STATES(73))) - ( STATES(7)*CONSTS(198)*STATES(73))/(CONSTS(197)+STATES(73))
RATES(74) = ((((((((( STATES(45)*CONSTS(193)*STATES(73))/(CONSTS(191)+STATES(73))+( STATES(55)*CONSTS(194)*STATES(73))/(CONSTS(192)+STATES(73))+( STATES(8)*CONSTS(200)*STATES(73))/(CONSTS(199)+STATES(73))) - CONSTS(195)*STATES(74))+( STATES(7)*CONSTS(198)*STATES(73))/(CONSTS(197)+STATES(73))) - STATES(58)*STATES(74)*CONSTS(161))+ STATES(59)*CONSTS(162)) - STATES(74)*STATES(62)*CONSTS(159))+ STATES(61)*CONSTS(160)) - STATES(74)*STATES(68)*CONSTS(171))+ STATES(71)*CONSTS(172)
RATES(75) = ((( STATES(45)*CONSTS(193)*STATES(73))/(CONSTS(191)+STATES(73))+( STATES(55)*CONSTS(194)*STATES(73))/(CONSTS(192)+STATES(73))+( STATES(7)*CONSTS(198)*STATES(73))/(CONSTS(197)+STATES(73))) - CONSTS(196)*STATES(75))+( STATES(8)*CONSTS(200)*STATES(73))/(CONSTS(199)+STATES(73))
RATES(77) = CONSTS(195)*STATES(74)
RATES(76) = CONSTS(196)*STATES(75)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END