C C There are a total of 207 entries in the algebraic variable array. C There are a total of 60 entries in each of the rate and state variable arrays. C There are a total of 388 entries in the constant variable array. C C C VOI is time in component Environment (second). C CONSTS(1) is ExpType in component Environment (dimensionless). C CONSTS(2) is StateType in component Environment (dimensionless). C CONSTS(3) is F in component parameters (kilojoule_per_mole_per_millivolt). C CONSTS(324) is V_cyto in component parameters (cyto_per_cell). C CONSTS(308) is V_mito in component parameters (mito_per_cell). C CONSTS(4) is W_x in component parameters (l_water_per_l_mito). C CONSTS(5) is W_i in component parameters (l_water_per_l_mito). C CONSTS(287) is W_c in component parameters (l_water_per_l_cyto). C CONSTS(6) is gamma in component parameters (per_micrometer). C CONSTS(7) is rho_m in component parameters (mg_per_mitochondria). C CONSTS(8) is n_A in component parameters (dimensionless). C CONSTS(9) is p_PI in component parameters (micrometer_per_second). C CONSTS(10) is p_A in component parameters (micrometer_per_second). C CONSTS(11) is p_TI in component parameters (micrometer_per_second). C CONSTS(12) is k_mADP in component parameters (molar). C CONSTS(13) is theta in component parameters (dimensionless). C CONSTS(14) is k_O2 in component parameters (molar). C CONSTS(15) is NAD_tot in component parameters (molar). C CONSTS(16) is Q_tot in component parameters (molar). C CONSTS(17) is cytC_tot in component parameters (molar). C CONSTS(18) is A_tot in component parameters (molar). C CONSTS(19) is CR_tot in component parameters (molar). C CONSTS(20) is CO2_tot_x in component parameters (molar). C CONSTS(21) is FAD_tot in component parameters (molar). C CONSTS(22) is C_IM in component parameters (mole_per_l_mito_volume_per_mv). C CONSTS(23) is X_pdh in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(24) is X_cits in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(25) is X_acon in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(26) is X_isod in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(27) is X_akgd in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(28) is K_ir_akgd in component parameters (molar). C CONSTS(29) is X_scoas in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(30) is X_sdh in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(31) is X_fum in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(32) is X_mdh in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(33) is X_ndk in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(34) is X_got in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(35) is X_PYRH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(36) is X_GLUH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(37) is X_CITMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(38) is X_AKGMAL in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(39) is X_SUCMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(40) is X_MALPI in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(41) is X_ASPGLU in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(42) is X_C1 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(43) is X_C3 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_half_molar). C CONSTS(44) is k_PI_1 in component parameters (molar). C CONSTS(45) is k_PI_2 in component parameters (molar). C CONSTS(46) is X_C4 in component parameters (mole_per_second_per_l_mito_volume_per_molar). C CONSTS(47) is X_F1 in component parameters (mole_per_second_per_l_mito_volume_per_molar). C CONSTS(48) is X_ANT in component parameters (mole_per_second_per_l_mito_volume). C CONSTS(49) is X_PiHt in component parameters (mole_per_second_per_l_mito_volume_per_molar). C CONSTS(50) is k_PiHt in component parameters (molar). C CONSTS(51) is X_KH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(52) is X_Hle in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_mv). C CONSTS(274) is X_HK in component parameters (mole_per_second_per_l_cyto_volume). C CONSTS(53) is J_AtC in component parameters (mole_per_second_per_l_cyto_volume). C CONSTS(275) is RT in component thermodynamic_data (kilojoule_per_mole). C CONSTS(54) is dGf1_H2O in component thermodynamic_data (kilojoule_per_mole). C CONSTS(55) is dGf1_O2 in component thermodynamic_data (kilojoule_per_mole). C CONSTS(56) is dGf1_NADH in component thermodynamic_data (kilojoule_per_mole). C CONSTS(57) is dGf1_NAD in component thermodynamic_data (kilojoule_per_mole). C CONSTS(58) is dGf1_QH2 in component thermodynamic_data (kilojoule_per_mole). C CONSTS(59) is dGf1_COQ in component thermodynamic_data (kilojoule_per_mole). C CONSTS(60) is dGf1_ATP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(61) is dGf1_ADP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(62) is dGf1_AMP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(63) is dGf1_GTP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(64) is dGf1_GDP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(65) is dGf1_Cred in component thermodynamic_data (kilojoule_per_mole). C CONSTS(66) is dGf1_Cox in component thermodynamic_data (kilojoule_per_mole). C CONSTS(67) is dGf1_Pi in component thermodynamic_data (kilojoule_per_mole). C CONSTS(68) is dGf1_PCr in component thermodynamic_data (kilojoule_per_mole). C CONSTS(69) is dGf1_Cr in component thermodynamic_data (kilojoule_per_mole). C CONSTS(70) is dGf1_FADH2 in component thermodynamic_data (kilojoule_per_mole). C CONSTS(71) is dGf1_FAD in component thermodynamic_data (kilojoule_per_mole). C CONSTS(72) is dGf1_COASH in component thermodynamic_data (kilojoule_per_mole). C CONSTS(73) is dGf1_ACCOA in component thermodynamic_data (kilojoule_per_mole). C CONSTS(74) is dGf1_OAA in component thermodynamic_data (kilojoule_per_mole). C CONSTS(75) is dGf1_CIT in component thermodynamic_data (kilojoule_per_mole). C CONSTS(76) is dGf1_ICIT in component thermodynamic_data (kilojoule_per_mole). C CONSTS(77) is dGf1_AKG in component thermodynamic_data (kilojoule_per_mole). C CONSTS(78) is dGf1_SCOA in component thermodynamic_data (kilojoule_per_mole). C CONSTS(79) is dGf1_SUC in component thermodynamic_data (kilojoule_per_mole). C CONSTS(80) is dGf1_FUM in component thermodynamic_data (kilojoule_per_mole). C CONSTS(81) is dGf1_MAL in component thermodynamic_data (kilojoule_per_mole). C CONSTS(82) is dGf1_ASP in component thermodynamic_data (kilojoule_per_mole). C CONSTS(83) is dGf1_GLU in component thermodynamic_data (kilojoule_per_mole). C CONSTS(84) is dGf1_CO2_tot in component thermodynamic_data (kilojoule_per_mole). C CONSTS(85) is dGf1_PYR in component thermodynamic_data (kilojoule_per_mole). C CONSTS(86) is dGf1_GLC in component thermodynamic_data (kilojoule_per_mole). C CONSTS(87) is dGf1_G6P in component thermodynamic_data (kilojoule_per_mole). C ALGBRC(1) is P_ATP_x in component P_ATP (dimensionless). C ALGBRC(164) is P_ATP_i in component P_ATP (dimensionless). C ALGBRC(18) is P_ATP_c in component P_ATP (dimensionless). C ALGBRC(19) is uATP_x in component P_ATP (molar). C ALGBRC(21) is mATP_c in component P_ATP (molar). C CONSTS(297) is K_H in component P_ATP (molar). C CONSTS(312) is K_Mg in component P_ATP (molar). C CONSTS(326) is K_K in component P_ATP (molar). C STATES(1) is ATP_x in component dATP_x_dt (molar). C STATES(2) is H_x in component dH_x_dt (molar). C STATES(3) is Mg_x in component dMg_x_dt (molar). C STATES(4) is K_x in component dK_x_dt (molar). C ALGBRC(133) is H_i in component dH_i_dt (molar). C ALGBRC(150) is Mg_i in component dMg_i_dt (molar). C ALGBRC(163) is K_i in component dK_i_dt (molar). C STATES(5) is ATP_c in component dATP_c_dt (molar). C STATES(6) is H_c in component dH_c_dt (molar). C STATES(7) is Mg_c in component dMg_c_dt (molar). C STATES(8) is K_c in component dK_c_dt (molar). C ALGBRC(22) is P_ADP_x in component P_ADP (dimensionless). C ALGBRC(166) is P_ADP_i in component P_ADP (dimensionless). C ALGBRC(23) is P_ADP_c in component P_ADP (dimensionless). C ALGBRC(24) is uADP_x in component P_ADP (molar). C ALGBRC(26) is mADP_c in component P_ADP (molar). C CONSTS(333) is K_H in component P_ADP (molar). C CONSTS(337) is K_Mg in component P_ADP (molar). C CONSTS(341) is K_K in component P_ADP (molar). C STATES(9) is ADP_x in component dADP_x_dt (molar). C STATES(10) is ADP_c in component dADP_c_dt (molar). C ALGBRC(27) is P_AMP_x in component P_AMP (dimensionless). C ALGBRC(168) is P_AMP_i in component P_AMP (dimensionless). C ALGBRC(28) is P_AMP_c in component P_AMP (dimensionless). C ALGBRC(29) is uAMP_x in component P_AMP (molar). C CONSTS(345) is K_H in component P_AMP (molar). C CONSTS(348) is K_Mg in component P_AMP (molar). C CONSTS(350) is K_K in component P_AMP (molar). C STATES(11) is AMP_x in component dAMP_x_dt (molar). C ALGBRC(30) is P_GTP_x in component P_GTP (dimensionless). C ALGBRC(165) is P_GTP_i in component P_GTP (dimensionless). C ALGBRC(20) is P_GTP_c in component P_GTP (dimensionless). C ALGBRC(31) is uGTP_x in component P_GTP (molar). C CONSTS(314) is K_H in component P_GTP (molar). C CONSTS(327) is K_Mg in component P_GTP (molar). C CONSTS(334) is K_K in component P_GTP (molar). C STATES(12) is GTP_x in component dGTP_x_dt (molar). C ALGBRC(32) is P_GDP_x in component P_GDP (dimensionless). C ALGBRC(167) is P_GDP_i in component P_GDP (dimensionless). C ALGBRC(25) is P_GDP_c in component P_GDP (dimensionless). C ALGBRC(33) is uGDP_x in component P_GDP (molar). C CONSTS(338) is K_H in component P_GDP (molar). C CONSTS(342) is K_Mg in component P_GDP (molar). C CONSTS(346) is K_K in component P_GDP (molar). C STATES(13) is GDP_x in component dGDP_x_dt (molar). C ALGBRC(34) is P_PI_x in component P_PI (dimensionless). C ALGBRC(169) is P_PI_i in component P_PI (dimensionless). C ALGBRC(2) is P_PI_c in component P_PI (dimensionless). C ALGBRC(35) is hPI_x in component P_PI (molar). C ALGBRC(170) is hPI_i in component P_PI (molar). C STATES(14) is PI_x in component dPI_x_dt (molar). C STATES(15) is PI_i in component dPI_i_dt (molar). C CONSTS(352) is K_H in component P_PI (molar). C CONSTS(354) is K_Mg in component P_PI (molar). C CONSTS(356) is K_K in component P_PI (molar). C ALGBRC(36) is P_COASH_x in component P_COASH (dimensionless). C ALGBRC(134) is P_COASH_i in component P_COASH (dimensionless). C ALGBRC(3) is P_COASH_c in component P_COASH (dimensionless). C CONSTS(358) is K_H in component P_COASH (molar). C CONSTS(88) is K_Mg in component P_COASH (molar). C CONSTS(89) is K_K in component P_COASH (molar). C ALGBRC(37) is P_OAA_x in component P_OAA (dimensionless). C ALGBRC(151) is P_OAA_i in component P_OAA (dimensionless). C ALGBRC(4) is P_OAA_c in component P_OAA (dimensionless). C CONSTS(360) is K_Mg in component P_OAA (molar). C CONSTS(90) is K_H in component P_OAA (molar). C CONSTS(91) is K_K in component P_OAA (molar). C ALGBRC(38) is P_CIT_x in component P_CIT (dimensionless). C ALGBRC(171) is P_CIT_i in component P_CIT (dimensionless). C ALGBRC(5) is P_CIT_c in component P_CIT (dimensionless). C ALGBRC(39) is uCIT_x in component P_CIT (molar). C ALGBRC(172) is uCIT_i in component P_CIT (molar). C ALGBRC(40) is hCIT_x in component P_CIT (molar). C ALGBRC(173) is hCIT_i in component P_CIT (molar). C STATES(16) is CIT_x in component dCIT_x_dt (molar). C STATES(17) is CIT_i in component dCIT_i_dt (molar). C CONSTS(361) is K_H in component P_CIT (molar). C CONSTS(362) is K_Mg in component P_CIT (molar). C CONSTS(363) is K_K in component P_CIT (molar). C ALGBRC(41) is P_ICIT_x in component P_ICIT (dimensionless). C ALGBRC(152) is P_ICIT_i in component P_ICIT (dimensionless). C ALGBRC(6) is P_ICIT_c in component P_ICIT (dimensionless). C CONSTS(364) is K_H in component P_ICIT (molar). C CONSTS(365) is K_Mg in component P_ICIT (molar). C CONSTS(92) is K_K in component P_ICIT (molar). C ALGBRC(42) is P_SCOA_x in component P_SCOA (dimensionless). C ALGBRC(135) is P_SCOA_i in component P_SCOA (dimensionless). C ALGBRC(7) is P_SCOA_c in component P_SCOA (dimensionless). C CONSTS(366) is K_H in component P_SCOA (molar). C CONSTS(93) is K_Mg in component P_SCOA (molar). C CONSTS(94) is K_K in component P_SCOA (molar). C ALGBRC(43) is P_SUC_x in component P_SUC (dimensionless). C ALGBRC(174) is P_SUC_i in component P_SUC (dimensionless). C ALGBRC(8) is P_SUC_c in component P_SUC (dimensionless). C CONSTS(367) is K_H in component P_SUC (molar). C CONSTS(368) is K_Mg in component P_SUC (molar). C CONSTS(369) is K_K in component P_SUC (molar). C ALGBRC(44) is P_FUM_x in component P_FUM (dimensionless). C ALGBRC(136) is P_FUM_i in component P_FUM (dimensionless). C ALGBRC(9) is P_FUM_c in component P_FUM (dimensionless). C CONSTS(370) is K_H in component P_FUM (molar). C CONSTS(95) is K_Mg in component P_FUM (molar). C CONSTS(96) is K_K in component P_FUM (molar). C ALGBRC(45) is P_MAL_x in component P_MAL (dimensionless). C ALGBRC(175) is P_MAL_i in component P_MAL (dimensionless). C ALGBRC(10) is P_MAL_c in component P_MAL (dimensionless). C CONSTS(371) is K_H in component P_MAL (molar). C CONSTS(372) is K_Mg in component P_MAL (molar). C CONSTS(373) is K_K in component P_MAL (molar). C ALGBRC(46) is P_CO2_tot_x in component P_CO2_tot (dimensionless). C ALGBRC(137) is P_CO2_tot_i in component P_CO2_tot (dimensionless). C ALGBRC(11) is P_CO2_tot_c in component P_CO2_tot (dimensionless). C CONSTS(374) is K_H in component P_CO2_tot (molar). C CONSTS(97) is K_Mg in component P_CO2_tot (molar). C CONSTS(98) is K_K in component P_CO2_tot (molar). C ALGBRC(47) is P_PYR_x in component P_PYR (dimensionless). C ALGBRC(153) is P_PYR_i in component P_PYR (dimensionless). C ALGBRC(12) is P_PYR_c in component P_PYR (dimensionless). C CONSTS(375) is K_Mg in component P_PYR (molar). C CONSTS(99) is K_H in component P_PYR (molar). C CONSTS(100) is K_K in component P_PYR (molar). C ALGBRC(13) is P_G6P_x in component P_G6P (dimensionless). C ALGBRC(138) is P_G6P_i in component P_G6P (dimensionless). C ALGBRC(48) is P_G6P_c in component P_G6P (dimensionless). C CONSTS(376) is K_H in component P_G6P (molar). C CONSTS(101) is K_Mg in component P_G6P (molar). C CONSTS(102) is K_K in component P_G6P (molar). C ALGBRC(49) is P_GLU_x in component P_GLU (dimensionless). C ALGBRC(154) is P_GLU_i in component P_GLU (dimensionless). C ALGBRC(14) is P_GLU_c in component P_GLU (dimensionless). C CONSTS(377) is K_H in component P_GLU (molar). C CONSTS(378) is K_Mg in component P_GLU (molar). C CONSTS(103) is K_K in component P_GLU (molar). C ALGBRC(50) is P_ASP_x in component P_ASP (dimensionless). C ALGBRC(155) is P_ASP_i in component P_ASP (dimensionless). C ALGBRC(15) is P_ASP_c in component P_ASP (dimensionless). C CONSTS(379) is K_H in component P_ASP (molar). C CONSTS(380) is K_Mg in component P_ASP (molar). C CONSTS(104) is K_K in component P_ASP (molar). C CONSTS(276) is P_GLC_c in component P_GLC (dimensionless). C CONSTS(105) is K_H in component P_GLC (molar). C CONSTS(106) is K_Mg in component P_GLC (molar). C CONSTS(107) is K_K in component P_GLC (molar). C CONSTS(277) is P_NAD_x in component P_NAD (dimensionless). C CONSTS(108) is K_H in component P_NAD (molar). C CONSTS(109) is K_Mg in component P_NAD (molar). C CONSTS(110) is K_K in component P_NAD (molar). C CONSTS(278) is P_NADH_x in component P_NADH (dimensionless). C CONSTS(111) is K_H in component P_NADH (molar). C CONSTS(112) is K_Mg in component P_NADH (molar). C CONSTS(113) is K_K in component P_NADH (molar). C CONSTS(305) is P_QH2_x in component P_QH2 (dimensionless). C CONSTS(114) is K_H in component P_QH2 (molar). C CONSTS(115) is K_Mg in component P_QH2 (molar). C CONSTS(116) is K_K in component P_QH2 (molar). C CONSTS(321) is P_COQ_x in component P_COQ (dimensionless). C CONSTS(117) is K_H in component P_COQ (molar). C CONSTS(118) is K_Mg in component P_COQ (molar). C CONSTS(119) is K_K in component P_COQ (molar). C CONSTS(279) is P_ACCOA_x in component P_ACCOA (dimensionless). C CONSTS(120) is K_H in component P_ACCOA (molar). C CONSTS(121) is K_Mg in component P_ACCOA (molar). C CONSTS(122) is K_K in component P_ACCOA (molar). C CONSTS(307) is P_AKG_x in component P_AKG (dimensionless). C CONSTS(123) is K_H in component P_AKG (molar). C CONSTS(124) is K_Mg in component P_AKG (molar). C CONSTS(125) is K_K in component P_AKG (molar). C CONSTS(280) is P_O2_x in component P_O2 (dimensionless). C CONSTS(126) is K_H in component P_O2 (molar). C CONSTS(127) is K_Mg in component P_O2 (molar). C CONSTS(128) is K_K in component P_O2 (molar). C CONSTS(281) is P_FADH2_x in component P_FADH2 (dimensionless). C CONSTS(129) is K_H in component P_FADH2 (molar). C CONSTS(130) is K_Mg in component P_FADH2 (molar). C CONSTS(131) is K_K in component P_FADH2 (molar). C CONSTS(282) is P_FAD_x in component P_FAD (dimensionless). C CONSTS(132) is K_H in component P_FAD (molar). C CONSTS(133) is K_Mg in component P_FAD (molar). C CONSTS(134) is K_K in component P_FAD (molar). C CONSTS(283) is P_H_x in component P_H (dimensionless). C CONSTS(135) is K_H in component P_H (molar). C CONSTS(136) is K_Mg in component P_H (molar). C CONSTS(137) is K_K in component P_H (molar). C CONSTS(284) is P_Mg_x in component P_Mg (dimensionless). C CONSTS(138) is K_H in component P_Mg (molar). C CONSTS(139) is K_Mg in component P_Mg (molar). C CONSTS(140) is K_K in component P_Mg (molar). C CONSTS(285) is P_K_x in component P_K (dimensionless). C CONSTS(141) is K_H in component P_K (molar). C CONSTS(142) is K_Mg in component P_K (molar). C CONSTS(143) is K_K in component P_K (molar). C ALGBRC(115) is J_pdh in component J_pdh (mole_per_second_per_l_mito_volume). C STATES(18) is PYR_x in component dPYR_x_dt (molar). C STATES(19) is COASH_x in component dCOASH_x_dt (molar). C ALGBRC(114) is NAD_x in component dNAD_x_dt (molar). C STATES(20) is CO2_tot_x in component dCO2_tot_x_dt (molar). C STATES(21) is ACCOA_x in component dACCOA_x_dt (molar). C STATES(22) is NADH_x in component dNADH_x_dt (molar). C ALGBRC(51) is K_eq_pdh in component J_pdh (dimensionless). C CONSTS(315) is K0_eq_pdh in component J_pdh (molar). C CONSTS(299) is dGr0_pdh in component J_pdh (kilojoule_per_mole). C CONSTS(144) is K_mA in component J_pdh (molar). C CONSTS(145) is K_mB in component J_pdh (molar). C CONSTS(146) is K_mC in component J_pdh (molar). C ALGBRC(52) is alpha_i1 in component J_pdh (dimensionless). C ALGBRC(53) is alpha_i2 in component J_pdh (dimensionless). C CONSTS(147) is K_iACCOA in component J_pdh (molar). C CONSTS(148) is K_iNADH in component J_pdh (molar). C CONSTS(149) is minCon in component J_pdh (molar). C ALGBRC(57) is J_cits in component J_cits (mole_per_second_per_l_mito_volume). C STATES(23) is OAA_x in component dOAA_x_dt (molar). C STATES(24) is SCOA_x in component dSCOA_x_dt (molar). C ALGBRC(54) is K_eq_cits in component J_cits (dimensionless). C CONSTS(316) is K0_eq_cits in component J_cits (molar_squared). C CONSTS(300) is dGr0_cits in component J_cits (kilojoule_per_mole). C CONSTS(150) is K_mA in component J_cits (molar). C CONSTS(151) is K_mB in component J_cits (molar). C CONSTS(152) is K_ia in component J_cits (molar). C CONSTS(153) is K_iCIT in component J_cits (molar). C CONSTS(154) is K_iATP in component J_cits (molar). C CONSTS(155) is K_iADP in component J_cits (molar). C CONSTS(156) is K_iAMP in component J_cits (molar). C CONSTS(157) is K_iCOASH in component J_cits (molar). C CONSTS(158) is K_iSCOA in component J_cits (molar). C ALGBRC(55) is alpha_i1 in component J_cits (dimensionless). C ALGBRC(56) is alpha_i2 in component J_cits (dimensionless). C ALGBRC(60) is J_acon in component J_acon (mole_per_second_per_l_mito_volume). C STATES(25) is ICIT_x in component dICIT_x_dt (molar). C ALGBRC(59) is V_mr in component J_acon (mole_per_second_per_l_mito_volume). C ALGBRC(58) is K_eq_acon in component J_acon (dimensionless). C CONSTS(317) is K0_eq_acon in component J_acon (dimensionless). C CONSTS(301) is dGr0_acon in component J_acon (kilojoule_per_mole). C CONSTS(159) is K_mA in component J_acon (molar). C CONSTS(160) is K_mP in component J_acon (molar). C ALGBRC(116) is J_isod in component J_isod (mole_per_second_per_l_mito_volume). C STATES(26) is AKG_x in component dAKG_x_dt (molar). C ALGBRC(61) is K_eq_isod in component J_isod (molar). C CONSTS(318) is K0_eq_isod in component J_isod (molar_cubed). C CONSTS(302) is dGr0_isod in component J_isod (kilojoule_per_mole). C CONSTS(161) is K_mA in component J_isod (molar). C CONSTS(162) is K_mB in component J_isod (molar). C CONSTS(163) is K_ib in component J_isod (molar). C CONSTS(164) is K_iq in component J_isod (molar). C ALGBRC(62) is alpha_i in component J_isod (dimensionless). C CONSTS(165) is n_H in component J_isod (dimensionless). C CONSTS(166) is K_iATP in component J_isod (molar). C CONSTS(167) is K_aADP in component J_isod (molar). C CONSTS(168) is minCon in component J_isod (molar). C ALGBRC(117) is J_akgd in component J_akgd (mole_per_second_per_l_mito_volume). C ALGBRC(63) is K_eq_akgd in component J_akgd (dimensionless). C CONSTS(319) is K0_eq_akgd in component J_akgd (molar). C CONSTS(303) is dGr0_akgd in component J_akgd (kilojoule_per_mole). C CONSTS(169) is K_mA in component J_akgd (molar). C CONSTS(170) is K_mB in component J_akgd (molar). C CONSTS(171) is K_mC in component J_akgd (molar). C CONSTS(172) is K_iq in component J_akgd (molar). C CONSTS(173) is K_ir1 in component J_akgd (molar). C CONSTS(174) is K_ir2 in component J_akgd (molar). C CONSTS(175) is K_iATP in component J_akgd (molar). C CONSTS(176) is K_aADP in component J_akgd (molar). C ALGBRC(64) is alpha_i in component J_akgd (dimensionless). C CONSTS(177) is minCon in component J_akgd (molar). C ALGBRC(67) is J_scoas in component J_scoas (mole_per_second_per_l_mito_volume). C STATES(27) is SUC_x in component dSUC_x_dt (molar). C ALGBRC(65) is K_eq_scoas in component J_scoas (dimensionless). C CONSTS(320) is K0_eq_scoas in component J_scoas (molar). C CONSTS(304) is dGr0_scoas in component J_scoas (kilojoule_per_mole). C ALGBRC(66) is V_mr in component J_scoas (mole_per_second_per_l_mito_volume). C CONSTS(178) is K_ia in component J_scoas (molar). C CONSTS(179) is K_ib in component J_scoas (molar). C CONSTS(180) is K_ic in component J_scoas (molar). C CONSTS(181) is K_ip in component J_scoas (molar). C CONSTS(182) is K_iq in component J_scoas (molar). C CONSTS(183) is K_ir in component J_scoas (molar). C CONSTS(184) is K_mA in component J_scoas (molar). C CONSTS(185) is K_mB in component J_scoas (molar). C CONSTS(186) is K_mC in component J_scoas (molar). C CONSTS(187) is K_mP in component J_scoas (molar). C CONSTS(188) is K_mQ in component J_scoas (molar). C CONSTS(189) is K_mR in component J_scoas (molar). C ALGBRC(120) is J_sdh in component J_sdh (mole_per_second_per_l_mito_volume). C ALGBRC(119) is COQ_x in component dCOQ_x_dt (molar). C STATES(28) is QH2_x in component dQH2_x_dt (molar). C STATES(29) is FUM_x in component dFUM_x_dt (molar). C ALGBRC(68) is K_eq_sdh in component J_sdh (dimensionless). C CONSTS(335) is K0_eq_sdh in component J_sdh (dimensionless). C CONSTS(328) is dGr0_sdh in component J_sdh (kilojoule_per_mole). C ALGBRC(69) is V_mr in component J_sdh (mole_per_second_per_l_mito_volume). C CONSTS(190) is K_ia in component J_sdh (molar). C CONSTS(191) is K_iq in component J_sdh (molar). C CONSTS(192) is K_mSUC in component J_sdh (molar). C CONSTS(193) is K_mCOQ in component J_sdh (molar). C CONSTS(194) is K_mQH2 in component J_sdh (molar). C CONSTS(195) is K_mFUM in component J_sdh (molar). C ALGBRC(70) is alpha_i in component J_sdh (dimensionless). C CONSTS(196) is K_iOAA in component J_sdh (molar). C CONSTS(197) is K_aSUC in component J_sdh (molar). C CONSTS(198) is K_aFUM in component J_sdh (molar). C ALGBRC(74) is J_fum in component J_fum (mole_per_second_per_l_mito_volume). C STATES(30) is MAL_x in component dMAL_x_dt (molar). C ALGBRC(71) is K_eq_fum in component J_fum (dimensionless). C CONSTS(343) is K0_eq_fum in component J_fum (dimensionless). C CONSTS(339) is dGr0_fum in component J_fum (kilojoule_per_mole). C ALGBRC(72) is V_mr in component J_fum (mole_per_second_per_l_mito_volume). C CONSTS(199) is K_mFUM in component J_fum (molar). C CONSTS(200) is K_mMAL in component J_fum (molar). C ALGBRC(73) is alpha_i in component J_fum (dimensionless). C CONSTS(201) is K_iCIT in component J_fum (molar). C CONSTS(202) is K_iATP in component J_fum (molar). C CONSTS(203) is K_iADP in component J_fum (molar). C CONSTS(204) is K_iGTP in component J_fum (molar). C CONSTS(205) is K_iGDP in component J_fum (molar). C ALGBRC(118) is J_mdh in component J_mdh (mole_per_second_per_l_mito_volume). C ALGBRC(75) is K_eq_mdh in component J_mdh (dimensionless). C CONSTS(322) is K0_eq_mdh in component J_mdh (molar). C CONSTS(306) is dGr0_mdh in component J_mdh (kilojoule_per_mole). C ALGBRC(76) is V_mr in component J_mdh (mole_per_second_per_l_mito_volume). C CONSTS(206) is K_ia in component J_mdh (molar). C CONSTS(207) is K_ib in component J_mdh (molar). C CONSTS(208) is K_ip in component J_mdh (molar). C CONSTS(209) is K_iq in component J_mdh (molar). C CONSTS(210) is K_mNAD in component J_mdh (molar). C CONSTS(211) is K_mMAL in component J_mdh (molar). C CONSTS(212) is K_mOAA in component J_mdh (molar). C CONSTS(213) is K_mNADH in component J_mdh (molar). C ALGBRC(77) is alpha_i in component J_mdh (dimensionless). C CONSTS(214) is K_iATP in component J_mdh (molar). C CONSTS(215) is K_iADP in component J_mdh (molar). C CONSTS(216) is K_iAMP in component J_mdh (molar). C ALGBRC(81) is J_ndk_f in component J_ndk (mole_per_second_per_l_mito_volume). C ALGBRC(82) is J_ndk in component J_ndk (mole_per_second_per_l_mito_volume). C ALGBRC(78) is K_eq_ndk in component J_ndk (dimensionless). C CONSTS(313) is K0_eq_ndk in component J_ndk (dimensionless). C CONSTS(298) is dGr0_ndk in component J_ndk (kilojoule_per_mole). C ALGBRC(79) is V_mr in component J_ndk (mole_per_second_per_l_mito_volume). C CONSTS(217) is K_ia in component J_ndk (molar). C CONSTS(218) is K_ib in component J_ndk (molar). C CONSTS(219) is K_ip in component J_ndk (molar). C CONSTS(220) is K_iq in component J_ndk (molar). C CONSTS(221) is K_mA in component J_ndk (molar). C CONSTS(222) is K_mB in component J_ndk (molar). C CONSTS(223) is K_mP in component J_ndk (molar). C CONSTS(224) is K_mQ in component J_ndk (molar). C ALGBRC(80) is alpha_i in component J_ndk (dimensionless). C CONSTS(225) is K_iAMP in component J_ndk (molar). C CONSTS(226) is minCon in component J_ndk (molar). C ALGBRC(86) is J_got_f in component J_got (mole_per_second_per_l_mito_volume). C ALGBRC(87) is J_got in component J_got (mole_per_second_per_l_mito_volume). C STATES(31) is ASP_x in component dASP_x_dt (molar). C STATES(32) is GLU_x in component dGLU_x_dt (molar). C ALGBRC(83) is K_eq_got in component J_got (dimensionless). C CONSTS(329) is K0_eq_got in component J_got (dimensionless). C CONSTS(323) is dGr0_got in component J_got (kilojoule_per_mole). C ALGBRC(84) is V_mr in component J_got (mole_per_second_per_l_mito_volume). C CONSTS(227) is K_ia in component J_got (molar). C CONSTS(228) is K_ib in component J_got (molar). C CONSTS(229) is K_ip in component J_got (molar). C CONSTS(230) is K_iq in component J_got (molar). C CONSTS(231) is K_mA in component J_got (molar). C CONSTS(232) is K_mB in component J_got (molar). C CONSTS(233) is K_mP in component J_got (molar). C CONSTS(234) is K_mQ in component J_got (molar). C ALGBRC(85) is alpha_i in component J_got (dimensionless). C CONSTS(235) is K_iAKG in component J_got (molar). C CONSTS(236) is minCon in component J_got (molar). C ALGBRC(156) is J_PYRH in component J_PYRH (mole_per_second_per_l_mito_volume). C STATES(33) is PYR_i in component dPYR_i_dt (molar). C ALGBRC(157) is J_GLUH in component J_GLUH (mole_per_second_per_l_mito_volume). C STATES(34) is GLU_i in component dGLU_i_dt (molar). C ALGBRC(176) is J_CITMAL in component J_CITMAL (mole_per_second_per_l_mito_volume). C STATES(35) is MAL_i in component dMAL_i_dt (molar). C ALGBRC(88) is J_AKGMAL in component J_AKGMAL (mole_per_second_per_l_mito_volume). C CONSTS(237) is K_mAKGi in component J_AKGMAL (molar). C CONSTS(238) is K_mAKGx in component J_AKGMAL (molar). C CONSTS(239) is K_mMALi in component J_AKGMAL (molar). C CONSTS(240) is K_mMALx in component J_AKGMAL (molar). C STATES(36) is AKG_i in component dAKG_i_dt (molar). C ALGBRC(177) is J_MALPI in component J_MALPI (mole_per_second_per_l_mito_volume). C ALGBRC(159) is J_ASPGLU in component J_ASPGLU (mole_per_second_per_l_mito_volume). C ALGBRC(158) is K_eqASPGLU in component J_ASPGLU (dimensionless). C CONSTS(286) is K_hASPGLU in component J_ASPGLU (molar). C CONSTS(241) is K_iASPi in component J_ASPGLU (molar). C CONSTS(242) is K_iASPx in component J_ASPGLU (molar). C CONSTS(243) is K_iGLUi in component J_ASPGLU (molar). C CONSTS(244) is K_iGLUx in component J_ASPGLU (molar). C STATES(37) is ASP_i in component dASP_i_dt (molar). C CONSTS(245) is m in component J_ASPGLU (dimensionless). C CONSTS(246) is pKa_ASPGLU in component J_ASPGLU (dimensionless). C STATES(38) is dPsi in component ddPsi_dt (millivolt). C ALGBRC(178) is J_SUCMAL in component J_SUCMAL (mole_per_second_per_l_mito_volume). C STATES(39) is SUC_i in component dSUC_i_dt (molar). C ALGBRC(89) is J_PYRt in component J_PYRt (mole_per_second_per_l_mito_volume). C STATES(40) is PYR_c in component dPYR_c_dt (molar). C ALGBRC(91) is J_CITt in component J_CITt (mole_per_second_per_l_mito_volume). C STATES(41) is CIT_c in component dCIT_c_dt (molar). C ALGBRC(92) is J_MALt in component J_MALt (mole_per_second_per_l_mito_volume). C STATES(42) is MAL_c in component dMAL_c_dt (molar). C ALGBRC(90) is J_AKGt in component J_AKGt (mole_per_second_per_l_mito_volume). C STATES(43) is AKG_c in component dAKG_c_dt (molar). C ALGBRC(93) is J_SUCt in component J_SUCt (mole_per_second_per_l_mito_volume). C STATES(44) is SUC_c in component dSUC_c_dt (molar). C ALGBRC(94) is J_GLUt in component J_GLUt (mole_per_second_per_l_mito_volume). C STATES(45) is GLU_c in component dGLU_c_dt (molar). C ALGBRC(95) is J_ASPt in component J_ASPt (mole_per_second_per_l_mito_volume). C STATES(46) is ASP_c in component dASP_c_dt (molar). C ALGBRC(16) is J_FUMt in component J_FUMt (mole_per_second_per_l_mito_volume). C CONSTS(247) is zero in component J_FUMt (micrometer_per_second). C STATES(47) is FUM_c in component dFUM_c_dt (molar). C STATES(48) is FUM_i in component dFUM_i_dt (molar). C ALGBRC(17) is J_ICITt in component J_ICITt (mole_per_second_per_l_mito_volume). C STATES(49) is ICIT_c in component dICIT_c_dt (molar). C STATES(50) is ICIT_i in component dICIT_i_dt (molar). C ALGBRC(96) is J_ADPt in component J_ADPt (mole_per_second_per_l_mito_volume). C STATES(51) is ADP_i in component dADP_i_dt (molar). C ALGBRC(97) is J_ATPt in component J_ATPt (mole_per_second_per_l_mito_volume). C STATES(52) is ATP_i in component dATP_i_dt (molar). C ALGBRC(98) is J_AMPt in component J_AMPt (mole_per_second_per_l_mito_volume). C STATES(53) is AMP_c in component dAMP_c_dt (molar). C STATES(54) is AMP_i in component dAMP_i_dt (molar). C ALGBRC(99) is J_PIt in component J_PIt (mole_per_second_per_l_mito_volume). C STATES(55) is PI_c in component dPI_c_dt (molar). C ALGBRC(140) is J_C1 in component J_C1 (mole_per_second_per_l_mito_volume). C CONSTS(288) is dGr0_C1 in component J_C1 (kilojoule_per_mole). C ALGBRC(101) is K_eq_C1 in component J_C1 (per_molar). C ALGBRC(139) is K_app_C1 in component J_C1 (dimensionless). C ALGBRC(123) is J_C3 in component J_C3 (mole_per_second_per_l_mito_volume). C CONSTS(289) is dGr0_C3 in component J_C3 (kilojoule_per_mole). C ALGBRC(103) is K_eq_C3 in component J_C3 (molar_squared). C ALGBRC(105) is K_app_C3 in component J_C3 (dimensionless). C ALGBRC(107) is QH2_x in component J_C3 (molar). C ALGBRC(121) is COQ_x in component J_C3 (molar). C STATES(56) is Cred_i in component dCred_i_dt (molar). C ALGBRC(122) is Cox_i in component dCox_i_dt (molar). C CONSTS(248) is minCon in component J_C3 (molar). C ALGBRC(142) is J_C4 in component J_C4 (mole_per_second_per_l_mito_volume). C CONSTS(290) is dGr0_C4 in component J_C4 (kilojoule_per_mole). C ALGBRC(109) is K_eq_C4 in component J_C4 (dimensionless). C ALGBRC(141) is K_app_C4 in component J_C4 (dimensionless). C ALGBRC(111) is O2_x in component J_C4 (molar). C STATES(57) is O2_x in component dO2_x_dt (molar). C CONSTS(249) is minCon in component J_C4 (molar). C ALGBRC(144) is J_F1 in component J_F1 (mole_per_second_per_l_mito_volume). C CONSTS(291) is dGr0_F1 in component J_F1 (kilojoule_per_mole). C ALGBRC(113) is K_eq_F1 in component J_F1 (per_molar_squared). C ALGBRC(143) is K_app_F1 in component J_F1 (per_molar). C ALGBRC(179) is J_ANT in component J_ANT (mole_per_second_per_l_mito_volume). C CONSTS(250) is minCond in component J_ANT (molar). C ALGBRC(180) is J_PIHt in component J_PIHt (mole_per_second_per_l_mito_volume). C ALGBRC(145) is J_Hle in component J_Hle (mole_per_second_per_l_mito_volume). C CONSTS(251) is minCond in component J_Hle (millivolt). C ALGBRC(182) is J_KH in component J_KH (mole_per_second_per_l_mito_volume). C ALGBRC(104) is J_HK in component J_HK (mole_per_second_per_l_cyto_volume). C CONSTS(309) is dGr0_HK in component J_HK (kilojoule_per_mole). C STATES(58) is GLC_c in component dGLC_c_dt (molar). C STATES(59) is G6P_c in component dG6P_c_dt (molar). C CONSTS(252) is K_ia in component J_HK (molar). C CONSTS(253) is K_ib in component J_HK (molar). C CONSTS(254) is K_ip in component J_HK (molar). C CONSTS(255) is K_iq in component J_HK (molar). C CONSTS(256) is K_mA in component J_HK (molar). C CONSTS(257) is K_mB in component J_HK (molar). C CONSTS(258) is K_mP in component J_HK (molar). C CONSTS(259) is K_mQ in component J_HK (molar). C CONSTS(260) is Ki_G6P in component J_HK (molar). C CONSTS(325) is K_eq_HK in component J_HK (molar). C ALGBRC(100) is Kapp_HK in component J_HK (dimensionless). C ALGBRC(102) is Kapp_HK_m in component J_HK (dimensionless). C ALGBRC(108) is J_CKc in component J_CKc (mole_per_second_per_l_cyto_volume). C CONSTS(261) is X_CK in component J_CKc (mole_per_second_per_l_cyto_volume_per_molar_per_molar). C CONSTS(310) is K_CK in component J_CKc (per_molar). C STATES(60) is PCr_c in component dPCr_c_dt (molar). C ALGBRC(106) is Cr_c in component J_CKc (molar). C ALGBRC(183) is J_AKi in component J_AK (mole_per_second_per_l_mito_volume). C ALGBRC(112) is J_AKc in component J_AK (mole_per_second_per_l_mito_volume). C CONSTS(262) is X_AK in component J_AK (mole_per_second_per_l_mito_volume_per_molar_per_molar). C CONSTS(381) is dGr0_AK in component J_AK (kilojoule_per_mole). C CONSTS(382) is K_eq_AK in component J_AK (dimensionless). C ALGBRC(181) is Kapp_AKi in component J_AK (dimensionless). C ALGBRC(110) is Kapp_AKc in component J_AK (dimensionless). C ALGBRC(147) is dNAD_x_dt in component dNAD_x_dt (molar_per_second). C ALGBRC(146) is dNADH_x_dt in component dNADH_x_dt (molar_per_second). C ALGBRC(149) is dCOQ_x_dt in component dCOQ_x_dt (molar_per_second). C ALGBRC(148) is dQH2_x_dt in component dQH2_x_dt (molar_per_second). C ALGBRC(184) is dATP_x_dt in component dATP_x_dt (molar_per_second). C ALGBRC(185) is dADP_x_dt in component dADP_x_dt (molar_per_second). C CONSTS(263) is dAMP_x_dt in component dAMP_x_dt (molar_per_second). C ALGBRC(124) is dGTP_x_dt in component dGTP_x_dt (molar_per_second). C ALGBRC(125) is dGDP_x_dt in component dGDP_x_dt (molar_per_second). C ALGBRC(186) is dPI_x_dt in component dPI_x_dt (molar_per_second). C ALGBRC(160) is dPYR_x_dt in component dPYR_x_dt (molar_per_second). C ALGBRC(126) is dACCOA_x_dt in component dACCOA_x_dt (molar_per_second). C ALGBRC(187) is dCIT_x_dt in component dCIT_x_dt (molar_per_second). C ALGBRC(127) is dICIT_x_dt in component dICIT_x_dt (molar_per_second). C ALGBRC(128) is dAKG_x_dt in component dAKG_x_dt (molar_per_second). C ALGBRC(129) is dSCOA_x_dt in component dSCOA_x_dt (molar_per_second). C ALGBRC(130) is dCOASH_x_dt in component dCOASH_x_dt (molar_per_second). C ALGBRC(188) is dSUC_x_dt in component dSUC_x_dt (molar_per_second). C ALGBRC(131) is dFUM_x_dt in component dFUM_x_dt (molar_per_second). C ALGBRC(189) is dMAL_x_dt in component dMAL_x_dt (molar_per_second). C ALGBRC(132) is dOAA_x_dt in component dOAA_x_dt (molar_per_second). C ALGBRC(161) is dGLU_x_dt in component dGLU_x_dt (molar_per_second). C ALGBRC(162) is dASP_x_dt in component dASP_x_dt (molar_per_second). C CONSTS(264) is dO2_x_dt in component dO2_x_dt (molar_per_second). C CONSTS(265) is dCO2_tot_x_dt in component dCO2_tot_x_dt (molar_per_second). C CONSTS(332) is Rm_cyto in component dPYR_c_dt (mito_per_cyto). C CONSTS(336) is Rm_cyto in component dCIT_c_dt (mito_per_cyto). C CONSTS(331) is Rm_cyto in component dICIT_c_dt (mito_per_cyto). C CONSTS(340) is Rm_cyto in component dAKG_c_dt (mito_per_cyto). C CONSTS(344) is Rm_cyto in component dSUC_c_dt (mito_per_cyto). C CONSTS(330) is Rm_cyto in component dFUM_c_dt (mito_per_cyto). C CONSTS(347) is Rm_cyto in component dMAL_c_dt (mito_per_cyto). C CONSTS(349) is Rm_cyto in component dGLU_c_dt (mito_per_cyto). C CONSTS(351) is Rm_cyto in component dASP_c_dt (mito_per_cyto). C CONSTS(353) is Rm_cyto in component dPI_c_dt (mito_per_cyto). C CONSTS(266) is conversion_factor in component dPI_c_dt (mito_per_cyto). C CONSTS(355) is Rm_cyto in component dATP_c_dt (mito_per_cyto). C CONSTS(267) is conversion_factor in component dATP_c_dt (mito_per_cyto). C CONSTS(357) is Rm_cyto in component dADP_c_dt (mito_per_cyto). C CONSTS(268) is conversion_factor in component dADP_c_dt (mito_per_cyto). C CONSTS(359) is Rm_cyto in component dAMP_c_dt (mito_per_cyto). C CONSTS(269) is conversion_factor in component dAMP_c_dt (mito_per_cyto). C CONSTS(311) is FAD_x in component dFAD_x_dt (molar). C CONSTS(270) is dFAD_x_dt in component dFAD_x_dt (molar_per_second). C CONSTS(271) is dFADH2_x_dt in component dFAD_x_dt (molar_per_second). C CONSTS(292) is FADH2_x in component dFAD_x_dt (molar). C CONSTS(293) is dH_x_dt in component dH_x_dt (molar_per_second). C ALGBRC(207) is D in component BINDING_IONS (dimensionless). C ALGBRC(201) is phi_H in component BINDING_IONS (molar_per_second). C ALGBRC(202) is phi_Mg in component BINDING_IONS (molar_per_second). C ALGBRC(203) is phi_K in component BINDING_IONS (molar_per_second). C ALGBRC(190) is dHB_dMg in component BINDING_IONS (dimensionless). C ALGBRC(191) is dHB_dK in component BINDING_IONS (dimensionless). C ALGBRC(194) is dMgB_dK in component BINDING_IONS (dimensionless). C ALGBRC(197) is dKB_dMg in component BINDING_IONS (dimensionless). C ALGBRC(206) is alpha_K in component BINDING_IONS (dimensionless). C ALGBRC(205) is alpha_Mg in component BINDING_IONS (dimensionless). C CONSTS(294) is dMg_x_dt in component dMg_x_dt (molar_per_second). C ALGBRC(204) is alpha_H in component BINDING_IONS (dimensionless). C ALGBRC(196) is dKB_dH in component BINDING_IONS (dimensionless). C ALGBRC(193) is dMgB_dH in component BINDING_IONS (dimensionless). C CONSTS(295) is dK_x_dt in component dK_x_dt (molar_per_second). C ALGBRC(192) is dHB_dH in component BINDING_IONS (dimensionless). C ALGBRC(195) is dMgB_dMg in component BINDING_IONS (dimensionless). C ALGBRC(198) is dKB_dK in component BINDING_IONS (dimensionless). C ALGBRC(199) is J_H_t in component BINDING_IONS (molar_per_second). C CONSTS(296) is J_Mg_t in component BINDING_IONS (molar_per_second). C ALGBRC(200) is J_K_t in component BINDING_IONS (molar_per_second). C CONSTS(272) is B_x in component BINDING_IONS (molar). C CONSTS(273) is K_BX in component BINDING_IONS (molar). C RATES(38) is d/dt dPsi in component ddPsi_dt (millivolt). C RATES(1) is d/dt ATP_x in component dATP_x_dt (molar). C RATES(9) is d/dt ADP_x in component dADP_x_dt (molar). C RATES(11) is d/dt AMP_x in component dAMP_x_dt (molar). C RATES(12) is d/dt GTP_x in component dGTP_x_dt (molar). C RATES(13) is d/dt GDP_x in component dGDP_x_dt (molar). C RATES(14) is d/dt PI_x in component dPI_x_dt (molar). C RATES(22) is d/dt NADH_x in component dNADH_x_dt (molar). C RATES(28) is d/dt QH2_x in component dQH2_x_dt (molar). C RATES(18) is d/dt PYR_x in component dPYR_x_dt (molar). C RATES(21) is d/dt ACCOA_x in component dACCOA_x_dt (molar). C RATES(16) is d/dt CIT_x in component dCIT_x_dt (molar). C RATES(25) is d/dt ICIT_x in component dICIT_x_dt (molar). C RATES(26) is d/dt AKG_x in component dAKG_x_dt (molar). C RATES(24) is d/dt SCOA_x in component dSCOA_x_dt (molar). C RATES(19) is d/dt COASH_x in component dCOASH_x_dt (molar). C RATES(27) is d/dt SUC_x in component dSUC_x_dt (molar). C RATES(29) is d/dt FUM_x in component dFUM_x_dt (molar). C RATES(30) is d/dt MAL_x in component dMAL_x_dt (molar). C RATES(23) is d/dt OAA_x in component dOAA_x_dt (molar). C RATES(32) is d/dt GLU_x in component dGLU_x_dt (molar). C RATES(31) is d/dt ASP_x in component dASP_x_dt (molar). C RATES(57) is d/dt O2_x in component dO2_x_dt (molar). C RATES(20) is d/dt CO2_tot_x in component dCO2_tot_x_dt (molar). C RATES(56) is d/dt Cred_i in component dCred_i_dt (molar). C RATES(52) is d/dt ATP_i in component dATP_i_dt (molar). C RATES(51) is d/dt ADP_i in component dADP_i_dt (molar). C RATES(54) is d/dt AMP_i in component dAMP_i_dt (molar). C RATES(15) is d/dt PI_i in component dPI_i_dt (molar). C RATES(33) is d/dt PYR_i in component dPYR_i_dt (molar). C RATES(17) is d/dt CIT_i in component dCIT_i_dt (molar). C RATES(50) is d/dt ICIT_i in component dICIT_i_dt (molar). C RATES(36) is d/dt AKG_i in component dAKG_i_dt (molar). C RATES(39) is d/dt SUC_i in component dSUC_i_dt (molar). C RATES(48) is d/dt FUM_i in component dFUM_i_dt (molar). C RATES(35) is d/dt MAL_i in component dMAL_i_dt (molar). C RATES(34) is d/dt GLU_i in component dGLU_i_dt (molar). C RATES(37) is d/dt ASP_i in component dASP_i_dt (molar). C RATES(40) is d/dt PYR_c in component dPYR_c_dt (molar). C RATES(41) is d/dt CIT_c in component dCIT_c_dt (molar). C RATES(49) is d/dt ICIT_c in component dICIT_c_dt (molar). C RATES(43) is d/dt AKG_c in component dAKG_c_dt (molar). C RATES(44) is d/dt SUC_c in component dSUC_c_dt (molar). C RATES(47) is d/dt FUM_c in component dFUM_c_dt (molar). C RATES(42) is d/dt MAL_c in component dMAL_c_dt (molar). C RATES(45) is d/dt GLU_c in component dGLU_c_dt (molar). C RATES(46) is d/dt ASP_c in component dASP_c_dt (molar). C RATES(6) is d/dt H_c in component dH_c_dt (molar). C RATES(7) is d/dt Mg_c in component dMg_c_dt (molar). C RATES(8) is d/dt K_c in component dK_c_dt (molar). C RATES(55) is d/dt PI_c in component dPI_c_dt (molar). C RATES(5) is d/dt ATP_c in component dATP_c_dt (molar). C RATES(10) is d/dt ADP_c in component dADP_c_dt (molar). C RATES(58) is d/dt GLC_c in component dGLC_c_dt (molar). C RATES(59) is d/dt G6P_c in component dG6P_c_dt (molar). C RATES(60) is d/dt PCr_c in component dPCr_c_dt (molar). C RATES(53) is d/dt AMP_c in component dAMP_c_dt (molar). C RATES(2) is d/dt H_x in component dH_x_dt (molar). C RATES(3) is d/dt Mg_x in component dMg_x_dt (molar). C RATES(4) is d/dt K_x in component dK_x_dt (molar). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 1 CONSTS(2) = 2 CONSTS(3) = 0.096484 CONSTS(4) = 0.651384 CONSTS(5) = 0.072376 CONSTS(6) = 5.99 CONSTS(7) = 2.725e5 CONSTS(8) = 3 CONSTS(9) = 327 CONSTS(10) = 85 CONSTS(11) = 85 CONSTS(12) = 3.5e-6 CONSTS(13) = 0.6 CONSTS(14) = 1.2e-4 CONSTS(15) = 2.97e-3 CONSTS(16) = 1.35e-3 CONSTS(17) = 2.7e-3 CONSTS(18) = 10e-3 CONSTS(19) = 42.7e-3 CONSTS(20) = 21.4e-3 CONSTS(21) = 1e-4 CONSTS(22) = 6.75e-6 CONSTS(23) = 0.205093643103884 CONSTS(24) = 9.82530750761898 CONSTS(25) = 0.027665411180519 CONSTS(26) = 0.492047553306503 CONSTS(27) = 0.087630003425524 CONSTS(28) = 6.04e-7 CONSTS(29) = 0.448520452008943 CONSTS(30) = 0.08578867843174 CONSTS(31) = 0.00707728025683 CONSTS(32) = 0.077060004692409 CONSTS(33) = 0.025543515147654 CONSTS(34) = 4.551391837332316 CONSTS(35) = 2.996779273769102e8 CONSTS(36) = 2.725114215177993e8 CONSTS(37) = 80.432110322866379 CONSTS(38) = 0.282091450119444 CONSTS(39) = 96.995618247389871 CONSTS(40) = 16.062441848715917 CONSTS(41) = 5.865621634754698e-5 CONSTS(42) = 3.236481869934019e4 CONSTS(43) = 0.790813746051711 CONSTS(44) = 3.248913595933763e-5 CONSTS(45) = 0.003474542589709 CONSTS(46) = 1.076088893810943e-4 CONSTS(47) = 8.120393248498723e3 CONSTS(48) = 0.007522124149343 CONSTS(49) = 3.335568500871044e7 CONSTS(50) = 0.00161251709168 CONSTS(51) = 4.758023126099895e6 CONSTS(52) = 3.02092849850504e2 CONSTS(53) = 0 CONSTS(54) = -235.74 CONSTS(55) = 16.4 CONSTS(56) = 39.31 CONSTS(57) = 18.1 CONSTS(58) = -23.3 CONSTS(59) = 65.17 CONSTS(60) = -2771 CONSTS(61) = -1903.96 CONSTS(62) = -1034.66 CONSTS(63) = -2771 CONSTS(64) = -1903.96 CONSTS(65) = -27.41 CONSTS(66) = -6.52 CONSTS(67) = -1098.27 CONSTS(68) = 0 CONSTS(69) = -252.68 CONSTS(70) = -67.6 CONSTS(71) = 19.55 CONSTS(72) = -0.72 CONSTS(73) = -178.19 CONSTS(74) = -794.74 CONSTS(75) = -1165.59 CONSTS(76) = -1158.94 CONSTS(77) = -793.41 CONSTS(78) = -507.55 CONSTS(79) = -690.44 CONSTS(80) = -603.32 CONSTS(81) = -842.66 CONSTS(82) = -692.26 CONSTS(83) = -692.4 CONSTS(84) = -530.71 CONSTS(85) = -470.82 CONSTS(86) = -907.21 CONSTS(87) = -1758.87 STATES(1) = 1.025e-11 STATES(2) = 5.7919e-8 STATES(3) = 0.0046388 STATES(4) = 0.086878 STATES(5) = 1e-9 STATES(6) = 1e-7 STATES(7) = 0.005 STATES(8) = 0.15 STATES(9) = 0.01 STATES(10) = 1e-9 STATES(11) = 1e-6 STATES(12) = 2.7157e-10 STATES(13) = 0.005 STATES(14) = 2.9918e-4 STATES(15) = 1.1574e-4 CONSTS(88) = 1e20 CONSTS(89) = 1e20 CONSTS(90) = 1e20 CONSTS(91) = 1e20 STATES(16) = 6.815100000000001e-5 STATES(17) = 1e-9 CONSTS(92) = 1e20 CONSTS(93) = 1e20 CONSTS(94) = 1e20 CONSTS(95) = 1e20 CONSTS(96) = 1e20 CONSTS(97) = 1e20 CONSTS(98) = 1e20 CONSTS(99) = 1e20 CONSTS(100) = 1e20 CONSTS(101) = 1e20 CONSTS(102) = 1e20 CONSTS(103) = 1e20 CONSTS(104) = 1e20 CONSTS(105) = 1e20 CONSTS(106) = 1e20 CONSTS(107) = 1e20 CONSTS(108) = 1e20 CONSTS(109) = 1e20 CONSTS(110) = 1e20 CONSTS(111) = 1e20 CONSTS(112) = 1e20 CONSTS(113) = 1e20 CONSTS(114) = 1e20 CONSTS(115) = 1e20 CONSTS(116) = 1e20 CONSTS(117) = 1e20 CONSTS(118) = 1e20 CONSTS(119) = 1e20 CONSTS(120) = 1e20 CONSTS(121) = 1e20 CONSTS(122) = 1e20 CONSTS(123) = 1e20 CONSTS(124) = 1e20 CONSTS(125) = 1e20 CONSTS(126) = 1e20 CONSTS(127) = 1e20 CONSTS(128) = 1e20 CONSTS(129) = 1e20 CONSTS(130) = 1e20 CONSTS(131) = 1e20 CONSTS(132) = 1e20 CONSTS(133) = 1e20 CONSTS(134) = 1e20 CONSTS(135) = 1e20 CONSTS(136) = 1e20 CONSTS(137) = 1e20 CONSTS(138) = 1e20 CONSTS(139) = 1e20 CONSTS(140) = 1e20 CONSTS(141) = 1e20 CONSTS(142) = 1e20 CONSTS(143) = 1e20 STATES(18) = 1.2076e-11 STATES(19) = 0.0030107 STATES(20) = 0.0214 STATES(21) = 2.5507e-12 STATES(22) = 8.484200000000001e-21 CONSTS(144) = 38.3e-6 CONSTS(145) = 9.9e-6 CONSTS(146) = 60.7e-6 CONSTS(147) = 40.2e-6 CONSTS(148) = 40e-6 CONSTS(149) = 1e-32 STATES(23) = 8.9236e-4 STATES(24) = 6.9627e-10 CONSTS(150) = 4e-6 CONSTS(151) = 14e-6 CONSTS(152) = 3.33e-6 CONSTS(153) = 1600e-6 CONSTS(154) = 900e-6 CONSTS(155) = 1800e-6 CONSTS(156) = 6000e-6 CONSTS(157) = 67e-6 CONSTS(158) = 140e-6 STATES(25) = 9.3558e-7 CONSTS(159) = 1161e-6 CONSTS(160) = 434e-6 STATES(26) = 4.7979e-11 CONSTS(161) = 74e-6 CONSTS(162) = 183e-6 CONSTS(163) = 23.8e-6 CONSTS(164) = 29e-6 CONSTS(165) = 3 CONSTS(166) = 91e-6 CONSTS(167) = 50e-6 CONSTS(168) = 1e-32 CONSTS(169) = 80e-6 CONSTS(170) = 55e-6 CONSTS(171) = 21e-6 CONSTS(172) = 6.9e-6 CONSTS(173) = 6.010527470838273e-7 CONSTS(174) = 1e3 CONSTS(175) = 50e-6 CONSTS(176) = 100e-6 CONSTS(177) = 1e-32 STATES(27) = 2.1681e-7 CONSTS(178) = 5.5e-6 CONSTS(179) = 100e-6 CONSTS(180) = 2000e-6 CONSTS(181) = 20e-6 CONSTS(182) = 3000e-6 CONSTS(183) = 11.1e-6 CONSTS(184) = 16e-6 CONSTS(185) = 55e-6 CONSTS(186) = 660e-6 CONSTS(187) = 20e-6 CONSTS(188) = 880e-6 CONSTS(189) = 11.1e-6 STATES(28) = 0 STATES(29) = 7.8476e-9 CONSTS(190) = 120e-6 CONSTS(191) = 1275e-6 CONSTS(192) = 467e-6 CONSTS(193) = 480e-6 CONSTS(194) = 2.45e-6 CONSTS(195) = 1200e-6 CONSTS(196) = 1.5e-6 CONSTS(197) = 450e-6 CONSTS(198) = 375e-6 STATES(30) = 3.2191e-8 CONSTS(199) = 44.7e-6 CONSTS(200) = 197.7e-6 CONSTS(201) = 3500e-6 CONSTS(202) = 40e-6 CONSTS(203) = 400e-6 CONSTS(204) = 80e-6 CONSTS(205) = 330e-6 CONSTS(206) = 279e-6 CONSTS(207) = 360e-6 CONSTS(208) = 5.5e-6 CONSTS(209) = 3.18e-6 CONSTS(210) = 90.55e-6 CONSTS(211) = 250e-6 CONSTS(212) = 6.128e-6 CONSTS(213) = 2.58e-6 CONSTS(214) = 183.2e-6 CONSTS(215) = 394.4e-6 CONSTS(216) = 420e-6 CONSTS(217) = 170e-6 CONSTS(218) = 143.6e-6 CONSTS(219) = 146.6e-6 CONSTS(220) = 156.5e-6 CONSTS(221) = 111e-6 CONSTS(222) = 100e-6 CONSTS(223) = 260e-6 CONSTS(224) = 278e-6 CONSTS(225) = 650e-6 CONSTS(226) = 1e-32 STATES(31) = 1.0217e-4 STATES(32) = 9.9311e-12 CONSTS(227) = 3480e-6 CONSTS(228) = 710e-6 CONSTS(229) = 50e-6 CONSTS(230) = 8400e-6 CONSTS(231) = 3900e-6 CONSTS(232) = 430e-6 CONSTS(233) = 88e-6 CONSTS(234) = 8900e-6 CONSTS(235) = 16.6e-3 CONSTS(236) = 1e-32 STATES(33) = 9.1804e-10 STATES(34) = 5.7539e-12 STATES(35) = 1e-9 CONSTS(237) = 0.3e-3 CONSTS(238) = 0.17e-3 CONSTS(239) = 1.4e-3 CONSTS(240) = 0.7e-3 STATES(36) = 9.999500000000001e-10 CONSTS(241) = 28e-6 CONSTS(242) = 2.8e-3 CONSTS(243) = 180e-6 CONSTS(244) = 1.6e-3 STATES(37) = 4.4404e-9 CONSTS(245) = 1.8 CONSTS(246) = 6.5 STATES(38) = -14.087 STATES(39) = 1e-9 STATES(40) = 1e-9 STATES(41) = 1e-9 STATES(42) = 1e-9 STATES(43) = 1e-9 STATES(44) = 1e-9 STATES(45) = 5.754e-12 STATES(46) = 4.4404e-9 CONSTS(247) = 0 STATES(47) = 1e-9 STATES(48) = 1e-9 STATES(49) = 1e-9 STATES(50) = 1e-9 STATES(51) = 1e-9 STATES(52) = 1e-9 STATES(53) = 1e-9 STATES(54) = 1e-9 STATES(55) = 1.1574e-4 STATES(56) = 9.6996e-8 CONSTS(248) = 1e-32 STATES(57) = 6.499999999999999e-5 CONSTS(249) = 1e-32 CONSTS(250) = 1e-9 CONSTS(251) = 1e-9 STATES(58) = 1e-9 STATES(59) = 1e-9 CONSTS(252) = 1e-3 CONSTS(253) = 47e-6 CONSTS(254) = 47e-6 CONSTS(255) = 1e-3 CONSTS(256) = 1e-3 CONSTS(257) = 47e-6 CONSTS(258) = 47e-6 CONSTS(259) = 1e-3 CONSTS(260) = 10e-6 CONSTS(261) = 1e7 STATES(60) = 1e-9 CONSTS(262) = 1e7 CONSTS(263) = 0 CONSTS(264) = 0 CONSTS(265) = 0 CONSTS(266) = 1 CONSTS(267) = 1 CONSTS(268) = 1 CONSTS(269) = 1 CONSTS(270) = 0 CONSTS(271) = 0 CONSTS(272) = 0.02 CONSTS(273) = 1e-7 CONSTS(274) = TERNRY(CONSTS(1).EQ.1.00000.AND.CONSTS(2).EQ.3.00000, 0.107108, 0.00000) CONSTS(275) = TERNRY(CONSTS(1).EQ.1.00000, ( 8.31400*(28.0000+273.150))/1000.00, ( 8.31400*(37.0000+273.150))/1000.00) CONSTS(276) = 1.00000 CONSTS(277) = 1.00000 CONSTS(278) = 1.00000 CONSTS(279) = 1.00000 CONSTS(280) = 1.00000 CONSTS(281) = 1.00000 CONSTS(282) = 1.00000 CONSTS(283) = 1.00000 CONSTS(284) = 1.00000 CONSTS(285) = 1.00000 CONSTS(286) = 1.00000*10.0000 ** - CONSTS(246) CONSTS(287) = TERNRY(CONSTS(1).EQ.1.00000, 80.0000, TERNRY(CONSTS(1).EQ.2.00000, 800.000, TERNRY(CONSTS(1).EQ.3.00000.OR.CONSTS(1).EQ.4.00000, 0.807000*1.04400, 80.0000) CONSTS(288) = ((CONSTS(57)+CONSTS(58)) - CONSTS(56)) - CONSTS(59) CONSTS(289) = ((CONSTS(59)+ 2.00000*CONSTS(65)) - CONSTS(58)) - 2.00000*CONSTS(66) CONSTS(290) = (( 2.00000*CONSTS(66)+CONSTS(54)) - 2.00000*CONSTS(65)) - 0.500000*CONSTS(55) CONSTS(291) = ((CONSTS(60)+CONSTS(54)) - CONSTS(61)) - CONSTS(67) CONSTS(292) = CONSTS(21)/2.00000 CONSTS(293) = 0.00000 CONSTS(294) = 0.00000 CONSTS(295) = 0.00000 CONSTS(296) = 0.00000 CONSTS(382) = 0.00000 CONSTS(383) = 0.00000 CONSTS(384) = 0.00000 CONSTS(385) = 0.00000 CONSTS(386) = 0.00000 CONSTS(387) = 0.00000 CONSTS(297) = 1.00000*10.0000 ** - 6.59000 CONSTS(298) = ((CONSTS(64)+CONSTS(60)) - CONSTS(63)) - CONSTS(61) CONSTS(299) = ((((CONSTS(84)+CONSTS(73)+CONSTS(56)) - CONSTS(85)) - CONSTS(72)) - CONSTS(57)) - CONSTS(54) CONSTS(300) = (((CONSTS(72)+CONSTS(75)) - CONSTS(73)) - CONSTS(74)) - CONSTS(54) CONSTS(301) = CONSTS(76) - CONSTS(75) CONSTS(302) = (((CONSTS(77)+CONSTS(56)+CONSTS(84)) - CONSTS(76)) - CONSTS(57)) - CONSTS(54) CONSTS(303) = ((((CONSTS(84)+CONSTS(78)+CONSTS(56)) - CONSTS(77)) - CONSTS(72)) - CONSTS(57)) - CONSTS(54) CONSTS(304) = (((CONSTS(72)+CONSTS(79)+CONSTS(63)) - CONSTS(64)) - CONSTS(78)) - CONSTS(67) CONSTS(305) = 1.00000 CONSTS(306) = ((CONSTS(74)+CONSTS(56)) - CONSTS(57)) - CONSTS(81) CONSTS(307) = 1.00000 CONSTS(308) = TERNRY(CONSTS(1).EQ.3.00000, 0.288200, TERNRY(CONSTS(1).EQ.4.00000, 0.0560000, 0.00000) CONSTS(309) = ((CONSTS(87)+CONSTS(61)) - CONSTS(86)) - CONSTS(60) CONSTS(310) = EXP(50.7800/CONSTS(275))/1.00000 CONSTS(311) = CONSTS(21) - CONSTS(292) CONSTS(312) = 1.00000*10.0000 ** - 3.82000 CONSTS(313) = EXP(- CONSTS(298)/CONSTS(275)) CONSTS(314) = CONSTS(297) CONSTS(315) = 1.00000*EXP(- CONSTS(299)/CONSTS(275)) CONSTS(316) = 1.00000*EXP(- CONSTS(300)/CONSTS(275)) CONSTS(317) = EXP(- CONSTS(301)/CONSTS(275)) CONSTS(318) = 1.00000*EXP(- CONSTS(302)/CONSTS(275)) CONSTS(319) = 1.00000*EXP(- CONSTS(303)/CONSTS(275)) CONSTS(320) = 1.00000*EXP(- CONSTS(304)/CONSTS(275)) CONSTS(321) = 1.00000 CONSTS(322) = 1.00000*EXP(- CONSTS(306)/CONSTS(275)) CONSTS(323) = ((CONSTS(74)+CONSTS(83)) - CONSTS(82)) - CONSTS(77) CONSTS(324) = TERNRY(CONSTS(1).EQ.3.00000, 0.680100, TERNRY(CONSTS(1).EQ.4.00000, 0.894000, 1.00000) CONSTS(325) = 1.00000*EXP(- CONSTS(309)/CONSTS(275)) CONSTS(326) = 1.00000*10.0000 ** - 1.01300 CONSTS(327) = CONSTS(312) CONSTS(328) = ((CONSTS(58)+CONSTS(80)) - CONSTS(79)) - CONSTS(59) CONSTS(329) = EXP(- CONSTS(323)/CONSTS(275)) CONSTS(330) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(331) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(332) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.3.00000.OR.CONSTS(1).EQ.4.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(333) = 1.00000*10.0000 ** - 6.42000 CONSTS(334) = CONSTS(326) CONSTS(335) = EXP(- CONSTS(328)/CONSTS(275)) CONSTS(336) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(337) = 1.00000*10.0000 ** - 2.79000 CONSTS(338) = CONSTS(333) CONSTS(339) = (CONSTS(81) - CONSTS(80)) - CONSTS(54) CONSTS(340) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(341) = 1.00000*10.0000 ** - 0.882000 CONSTS(342) = CONSTS(337) CONSTS(343) = EXP(- CONSTS(339)/CONSTS(275)) CONSTS(344) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(345) = 1.00000*10.0000 ** - 6.22000 CONSTS(346) = CONSTS(341) CONSTS(347) = TERNRY(CONSTS(2).EQ.0.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(348) = 1.00000*10.0000 ** - 1.86000 CONSTS(349) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(350) = 1.00000*10.0000 ** - 0.621500 CONSTS(351) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, 1.00000, CONSTS(308)/CONSTS(324)) CONSTS(352) = 1.00000*10.0000 ** - 6.71000 CONSTS(353) = CONSTS(308)/CONSTS(324) CONSTS(354) = 1.00000*10.0000 ** - 1.69000 CONSTS(355) = CONSTS(308)/CONSTS(324) CONSTS(356) = 1.00000*10.0000 ** 0.00740000 CONSTS(357) = CONSTS(308)/CONSTS(324) CONSTS(358) = 1.00000*10.0000 ** - 8.13000 CONSTS(359) = CONSTS(308)/CONSTS(324) CONSTS(360) = 1.00000*10.0000 ** - 0.862900 CONSTS(361) = 1.00000*10.0000 ** - 5.63000 CONSTS(362) = 1.00000*10.0000 ** - 3.37000 CONSTS(363) = 1.00000*10.0000 ** - 0.339000 CONSTS(364) = 1.00000*10.0000 ** - 5.64000 CONSTS(365) = 1.00000*10.0000 ** - 2.46000 CONSTS(366) = 1.00000*10.0000 ** - 3.96000 CONSTS(367) = 1.00000*10.0000 ** - 5.13000 CONSTS(368) = 1.00000*10.0000 ** - 1.17000 CONSTS(369) = 1.00000*10.0000 ** - 0.352500 CONSTS(370) = 1.00000*10.0000 ** - 4.10000 CONSTS(371) = 1.00000*10.0000 ** - 4.75000 CONSTS(372) = 1.00000*10.0000 ** - 1.55000 CONSTS(373) = 1.00000*10.0000 ** 0.107000 CONSTS(374) = 1.00000*10.0000 ** - 9.82000 CONSTS(375) = 1.00000*10.0000 ** - 1.02000 CONSTS(376) = 1.00000*10.0000 ** - 5.91000 CONSTS(377) = 1.00000*10.0000 ** - 4.06000 CONSTS(378) = 1.00000*10.0000 ** - 1.82000 CONSTS(379) = 1.00000*10.0000 ** - 3.65000 CONSTS(380) = 1.00000*10.0000 ** - 2.32000 CONSTS(381) = (CONSTS(60)+CONSTS(62)) - 2.00000*CONSTS(61) CONSTS(382) = EXP(- CONSTS(381)/CONSTS(275)) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(11) = CONSTS(382) RATES(57) = CONSTS(383) RATES(20) = CONSTS(384) RATES(6) = CONSTS(385) RATES(7) = CONSTS(386) RATES(8) = CONSTS(387) ALGBRC(17) = CONSTS(6)*CONSTS(11)*(STATES(49) - STATES(50)) RATES(50) = ALGBRC(17)/CONSTS(5) ALGBRC(16) = CONSTS(6)*CONSTS(247)*(STATES(47) - STATES(48)) RATES(48) = ALGBRC(16)/CONSTS(5) RATES(49) = ( - CONSTS(331)*ALGBRC(17))/CONSTS(287) RATES(47) = ( - CONSTS(330)*ALGBRC(16))/CONSTS(287) ALGBRC(1) = 1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326) ALGBRC(30) = ALGBRC(1) ALGBRC(22) = 1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341) ALGBRC(32) = ALGBRC(22) ALGBRC(34) = 1.00000+STATES(2)/CONSTS(352)+STATES(3)/CONSTS(354)+STATES(4)/CONSTS(356) ALGBRC(36) = 1.00000+STATES(2)/CONSTS(358) ALGBRC(42) = 1.00000+STATES(2)/CONSTS(366) ALGBRC(43) = 1.00000+STATES(2)/CONSTS(367)+STATES(3)/CONSTS(368)+STATES(4)/CONSTS(369) ALGBRC(65) = ( CONSTS(320)*ALGBRC(36)*ALGBRC(43)*ALGBRC(30))/( STATES(2)*ALGBRC(32)*ALGBRC(42)*ALGBRC(34)) ALGBRC(66) = ( CONSTS(29)*CONSTS(187)*CONSTS(182)*CONSTS(183))/( ALGBRC(65)*CONSTS(178)*CONSTS(179)*CONSTS(186)) ALGBRC(67) = ( CONSTS(29)*ALGBRC(66)*( STATES(13)*STATES(24)*STATES(14) - ( STATES(19)*STATES(27)*STATES(12))/ALGBRC(65)))/( ALGBRC(66)*CONSTS(178)*CONSTS(179)*CONSTS(186)+ ALGBRC(66)*CONSTS(179)*CONSTS(186)*STATES(13)+ ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)+ ALGBRC(66)*CONSTS(186)*STATES(13)*STATES(24)+ ALGBRC(66)*CONSTS(185)*STATES(13)*STATES(14)+ ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)+ ALGBRC(66)*STATES(13)*STATES(24)*STATES(14)+( CONSTS(29)*CONSTS(183)*CONSTS(188)*STATES(19))/ALGBRC(65)+( CONSTS(29)*CONSTS(182)*CONSTS(187)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(189)*STATES(19)*STATES(27))/ALGBRC(65)+( CONSTS(29)*CONSTS(188)*STATES(19)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(187)*STATES(27)*STATES(12))/ALGBRC(65)+( CONSTS(29)*STATES(19)*STATES(27)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(188)*CONSTS(183)*STATES(13)*STATES(19))/( CONSTS(178)*ALGBRC(65))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(12))/CONSTS(183)+( CONSTS(29)*CONSTS(188)*CONSTS(183)*STATES(13)*STATES(24)*STATES(19))/( CONSTS(178)*CONSTS(179)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(12))/CONSTS(183)+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(19)*STATES(27))/( CONSTS(178)*ALGBRC(65))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(27)*STATES(12))/( CONSTS(182)*CONSTS(183))+( CONSTS(29)*CONSTS(183)*CONSTS(188)*STATES(13)*STATES(24)*STATES(14)*STATES(19))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( CONSTS(29)*CONSTS(181)*CONSTS(189)*STATES(13)*STATES(24)*STATES(14)*STATES(27))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(24)*STATES(19)*STATES(27))/( CONSTS(178)*CONSTS(179)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(27)*STATES(12))/( CONSTS(182)*CONSTS(183))+( ALGBRC(66)*CONSTS(184)*CONSTS(180)*STATES(24)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(24)*STATES(14)*STATES(19)*STATES(27))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))) ALGBRC(78) = ( CONSTS(313)*ALGBRC(32)*ALGBRC(1))/( ALGBRC(30)*ALGBRC(22)) ALGBRC(79) = ( CONSTS(33)*CONSTS(224)*CONSTS(219))/( ALGBRC(78)*CONSTS(217)*CONSTS(222)) ALGBRC(27) = 1.00000+STATES(2)/CONSTS(345)+STATES(3)/CONSTS(348)+STATES(4)/CONSTS(350) ALGBRC(29) = ( STATES(11)*(1.00000+STATES(2)/CONSTS(345)))/ALGBRC(27) ALGBRC(80) = 1.00000+ALGBRC(29)/CONSTS(225) ALGBRC(81) = TERNRY(STATES(12).GT.CONSTS(226).AND.STATES(9).GT.CONSTS(226), ( (( CONSTS(33)*ALGBRC(79))/ALGBRC(80))*STATES(12)*STATES(9))/( ALGBRC(79)*CONSTS(222)*STATES(12)+ ALGBRC(79)*CONSTS(221)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(13))/ALGBRC(78)+( CONSTS(33)*CONSTS(223)*STATES(1))/ALGBRC(78)+ ALGBRC(79)*STATES(12)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(12)*STATES(13))/( ALGBRC(78)*CONSTS(217))+( CONSTS(33)*STATES(13)*STATES(1))/ALGBRC(78)+( ALGBRC(79)*CONSTS(221)*STATES(9)*STATES(1))/CONSTS(220)), 0.00000) ALGBRC(82) = TERNRY(STATES(13).GT.CONSTS(226).AND.STATES(1).GT.CONSTS(226), ALGBRC(81) - (( (( CONSTS(33)*ALGBRC(79))/ALGBRC(80))*STATES(13)*STATES(1))/ALGBRC(78))/( ALGBRC(79)*CONSTS(222)*STATES(12)+ ALGBRC(79)*CONSTS(221)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(13))/ALGBRC(78)+( CONSTS(33)*CONSTS(223)*STATES(1))/ALGBRC(78)+ ALGBRC(79)*STATES(12)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(12)*STATES(13))/( ALGBRC(78)*CONSTS(217))+( CONSTS(33)*STATES(13)*STATES(1))/ALGBRC(78)+( ALGBRC(79)*CONSTS(221)*STATES(9)*STATES(1))/CONSTS(220)), ALGBRC(81)) RATES(12) = (ALGBRC(67) - ALGBRC(82))/CONSTS(4) RATES(13) = (- ALGBRC(67)+ALGBRC(82))/CONSTS(4) ALGBRC(88) = ( CONSTS(38)*( STATES(36)*STATES(30) - STATES(26)*STATES(35)))/( CONSTS(237)*CONSTS(240)*(2.00000+STATES(35)/CONSTS(239)+STATES(30)/CONSTS(240)+STATES(36)/CONSTS(237)+STATES(26)/CONSTS(238)+( STATES(35)*STATES(26))/( CONSTS(239)*CONSTS(238))+( STATES(30)*STATES(36))/( CONSTS(240)*CONSTS(237)))) ALGBRC(90) = CONSTS(6)*CONSTS(11)*(STATES(43) - STATES(36)) RATES(36) = (- ALGBRC(88)+ALGBRC(90))/CONSTS(5) ALGBRC(89) = CONSTS(6)*CONSTS(11)*(STATES(40) - STATES(33)) RATES(40) = ( - CONSTS(332)*ALGBRC(89))/CONSTS(287) RATES(43) = ( - CONSTS(340)*ALGBRC(90))/CONSTS(287) ALGBRC(91) = CONSTS(6)*CONSTS(11)*(STATES(41) - STATES(17)) RATES(41) = ( - CONSTS(336)*ALGBRC(91))/CONSTS(287) ALGBRC(92) = CONSTS(6)*CONSTS(11)*(STATES(42) - STATES(35)) RATES(42) = ( - CONSTS(347)*ALGBRC(92))/CONSTS(287) ALGBRC(93) = CONSTS(6)*CONSTS(11)*(STATES(44) - STATES(39)) RATES(44) = ( - CONSTS(344)*ALGBRC(93))/CONSTS(287) ALGBRC(94) = CONSTS(6)*CONSTS(11)*(STATES(45) - STATES(34)) RATES(45) = ( - CONSTS(349)*ALGBRC(94))/CONSTS(287) ALGBRC(95) = CONSTS(6)*CONSTS(11)*(STATES(46) - STATES(37)) RATES(46) = ( - CONSTS(351)*ALGBRC(95))/CONSTS(287) ALGBRC(99) = CONSTS(6)*CONSTS(9)*(STATES(55) - STATES(15)) RATES(55) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, ( - ALGBRC(99)*CONSTS(266))/CONSTS(287), ( CONSTS(353)*ALGBRC(99)+CONSTS(53))/CONSTS(287)) ALGBRC(18) = 1.00000+STATES(6)/CONSTS(297)+STATES(7)/CONSTS(312)+STATES(8)/CONSTS(326) ALGBRC(21) = (( STATES(5)*STATES(7))/CONSTS(312))/ALGBRC(18) ALGBRC(23) = 1.00000+STATES(6)/CONSTS(333)+STATES(7)/CONSTS(337)+STATES(8)/CONSTS(341) ALGBRC(26) = (( STATES(10)*STATES(7))/CONSTS(337))/ALGBRC(23) ALGBRC(48) = 1.00000+STATES(6)/CONSTS(376) ALGBRC(100) = ( CONSTS(325)*ALGBRC(48)*ALGBRC(23))/( STATES(6)*ALGBRC(18)*CONSTS(276)) ALGBRC(102) = ( ALGBRC(100)*CONSTS(312)*ALGBRC(18))/( CONSTS(337)*ALGBRC(23)) ALGBRC(104) = ( (CONSTS(274)/( CONSTS(253)*CONSTS(256)))*( ALGBRC(21)*STATES(58) - ( STATES(59)*ALGBRC(26))/ALGBRC(102)))/(1.00000+ALGBRC(21)/CONSTS(252)+STATES(58)/CONSTS(253)+( ALGBRC(21)*STATES(58))/( CONSTS(253)*CONSTS(256))+STATES(59)/CONSTS(254)+ALGBRC(26)/CONSTS(255)+( STATES(59)*ALGBRC(26))/( CONSTS(255)*CONSTS(258))+( STATES(59)*STATES(58))/( CONSTS(260)*CONSTS(253))) RATES(58) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, ( - (1.00000+CONSTS(287)/1.00000)*ALGBRC(104))/CONSTS(287), 0.00000) RATES(59) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, ( (1.00000+CONSTS(287)/1.00000)*ALGBRC(104))/CONSTS(287), 0.00000) ALGBRC(106) = CONSTS(19) - STATES(60) ALGBRC(108) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(261)*( (( CONSTS(310)*STATES(10))/ALGBRC(23))*STATES(60)*STATES(6) - (STATES(5)/ALGBRC(18))*ALGBRC(106))) RATES(60) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, 0.00000, - ALGBRC(108)/CONSTS(287)) ALGBRC(97) = CONSTS(6)*CONSTS(10)*(STATES(5) - STATES(52)) ALGBRC(28) = 1.00000+STATES(6)/CONSTS(345)+STATES(7)/CONSTS(348)+STATES(8)/CONSTS(350) ALGBRC(110) = ( CONSTS(382)*ALGBRC(18)*ALGBRC(28))/( ALGBRC(23)*ALGBRC(23)) ALGBRC(112) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(262)*( ALGBRC(110)*STATES(10)*STATES(10) - STATES(53)*STATES(5))) RATES(5) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, ( - (1.00000+CONSTS(287)/1.00000)*ALGBRC(104) - CONSTS(267)*ALGBRC(97))/CONSTS(287), (( - CONSTS(355)*ALGBRC(97)+ALGBRC(108)+ CONSTS(267)*ALGBRC(112)) - CONSTS(53))/CONSTS(287)) ALGBRC(96) = CONSTS(6)*CONSTS(10)*(STATES(10) - STATES(51)) RATES(10) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, ( (1.00000+CONSTS(287)/1.00000)*ALGBRC(104) - CONSTS(268)*ALGBRC(96))/CONSTS(287), ((( - CONSTS(357)*ALGBRC(96) - ALGBRC(108)) - 2.00000*CONSTS(268)*ALGBRC(112))+CONSTS(53))/CONSTS(287)) ALGBRC(98) = CONSTS(6)*CONSTS(10)*(STATES(53) - STATES(54)) RATES(53) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000.OR.CONSTS(1).EQ.6.00000, 0.00000, ( - CONSTS(359)*ALGBRC(98)+ CONSTS(269)*ALGBRC(112))/CONSTS(287)) ALGBRC(114) = CONSTS(15) - STATES(22) ALGBRC(46) = 1.00000+STATES(2)/CONSTS(374) ALGBRC(47) = 1.00000+STATES(3)/CONSTS(375) ALGBRC(51) = ( CONSTS(315)*ALGBRC(46)*CONSTS(279)*CONSTS(278))/( STATES(2)*ALGBRC(47)*ALGBRC(36)*CONSTS(277)) ALGBRC(52) = 1.00000+STATES(21)/CONSTS(147) ALGBRC(53) = 1.00000+STATES(22)/CONSTS(148) ALGBRC(115) = TERNRY(STATES(18).GT.CONSTS(149).AND.STATES(19).GT.CONSTS(149).AND.ALGBRC(114).GT.CONSTS(149), ( CONSTS(23)*( STATES(18)*STATES(19)*ALGBRC(114) - ( STATES(20)*STATES(21)*STATES(22))/ALGBRC(51)))/( CONSTS(146)*ALGBRC(53)*STATES(18)*STATES(19)+ CONSTS(145)*ALGBRC(52)*STATES(18)*ALGBRC(114)+ CONSTS(144)*STATES(19)*ALGBRC(114)+ STATES(18)*STATES(19)*ALGBRC(114)), 0.00000) ALGBRC(37) = 1.00000+STATES(3)/CONSTS(360) ALGBRC(38) = 1.00000+STATES(2)/CONSTS(361)+STATES(3)/CONSTS(362)+STATES(4)/CONSTS(363) ALGBRC(54) = ( CONSTS(316)*ALGBRC(38)*ALGBRC(36))/( STATES(2) ** 2.00000*CONSTS(279)*ALGBRC(37)) ALGBRC(39) = ( STATES(16)*(1.00000+STATES(2)/CONSTS(361)))/ALGBRC(38) ALGBRC(55) = 1.00000+ALGBRC(39)/CONSTS(153) ALGBRC(19) = ( STATES(1)*(1.00000+STATES(2)/CONSTS(297)))/ALGBRC(1) ALGBRC(24) = ( STATES(9)*(1.00000+STATES(2)/CONSTS(333)))/ALGBRC(22) ALGBRC(56) = 1.00000+ALGBRC(19)/CONSTS(154)+ALGBRC(24)/CONSTS(155)+ALGBRC(29)/CONSTS(156)+STATES(19)/CONSTS(157)+STATES(24)/CONSTS(158) ALGBRC(57) = ( CONSTS(24)*( STATES(21)*STATES(23) - ( STATES(19)*STATES(16))/ALGBRC(54)))/( CONSTS(152)*CONSTS(151)*ALGBRC(55)+ CONSTS(150)*ALGBRC(55)*STATES(21)+ CONSTS(151)*ALGBRC(56)*STATES(23)+ STATES(23)*STATES(21)) RATES(21) = (ALGBRC(115) - ALGBRC(57))/CONSTS(4) ALGBRC(41) = 1.00000+STATES(2)/CONSTS(364)+STATES(3)/CONSTS(365) ALGBRC(58) = ( CONSTS(317)*ALGBRC(41))/ALGBRC(38) ALGBRC(59) = ( CONSTS(25)*CONSTS(160))/( CONSTS(159)*ALGBRC(58)) ALGBRC(60) = ( CONSTS(25)*ALGBRC(59)*(STATES(16) - STATES(25)/ALGBRC(58)))/( CONSTS(159)*ALGBRC(59)+ ALGBRC(59)*STATES(16)+( CONSTS(25)*STATES(25))/ALGBRC(58)) ALGBRC(61) = ( CONSTS(318)*ALGBRC(46)*CONSTS(307)*CONSTS(278))/( STATES(2) ** 2.00000*CONSTS(277)*ALGBRC(41)) ALGBRC(62) = 1.00000+ (CONSTS(167)/ALGBRC(24))*(1.00000+ALGBRC(19)/CONSTS(166)) ALGBRC(116) = TERNRY(ALGBRC(114).GT.CONSTS(168).AND.STATES(25).GT.CONSTS(168), ( CONSTS(26)*(1.00000 - ( STATES(20)*STATES(26)*STATES(22))/( ALGBRC(61)*ALGBRC(114)*STATES(25))))/(1.00000+ CONSTS(162)/STATES(25) ** CONSTS(165)*ALGBRC(62)+ (CONSTS(161)/ALGBRC(114))*(1.00000+ CONSTS(163)/STATES(25) ** CONSTS(165)*ALGBRC(62)+( STATES(22)*ALGBRC(62))/CONSTS(164))), 0.00000) RATES(25) = (ALGBRC(60) - ALGBRC(116))/CONSTS(4) ALGBRC(63) = ( CONSTS(319)*ALGBRC(46)*ALGBRC(42)*CONSTS(278))/( STATES(2)*ALGBRC(36)*CONSTS(307)*CONSTS(277)) ALGBRC(64) = 1.00000+ (CONSTS(176)/ALGBRC(24))*(1.00000+ALGBRC(19)/CONSTS(175)) ALGBRC(117) = TERNRY(STATES(26).GT.CONSTS(177).AND.STATES(19).GT.CONSTS(177).AND.ALGBRC(114).GT.CONSTS(177), (( CONSTS(27)*(1.00000 - ( STATES(20)*STATES(24)*STATES(22))/( ALGBRC(63)*STATES(26)*STATES(19)*ALGBRC(114))))/(1.00000+ (CONSTS(169)/STATES(26))*ALGBRC(64)+ (CONSTS(170)/STATES(19))*(1.00000+STATES(24)/CONSTS(172))+ (CONSTS(171)/ALGBRC(114))*(1.00000+STATES(22)/CONSTS(173))))/(1.00000+STATES(22)/CONSTS(174)), 0.00000) ALGBRC(49) = 1.00000+STATES(2)/CONSTS(377)+STATES(3)/CONSTS(378) ALGBRC(50) = 1.00000+STATES(2)/CONSTS(379)+STATES(3)/CONSTS(380) ALGBRC(83) = ( CONSTS(329)*ALGBRC(49)*ALGBRC(37))/( ALGBRC(50)*CONSTS(307)) ALGBRC(84) = ( CONSTS(34)*CONSTS(234)*CONSTS(229))/( ALGBRC(83)*CONSTS(227)*CONSTS(232)) ALGBRC(85) = 1.00000+STATES(26)/CONSTS(235) ALGBRC(86) = TERNRY(STATES(31).GT.CONSTS(236).AND.STATES(26).GT.CONSTS(236), ( CONSTS(34)*ALGBRC(84)*STATES(31)*STATES(26))/( ALGBRC(84)*CONSTS(232)*STATES(31)+ ALGBRC(84)*CONSTS(231)*ALGBRC(85)*STATES(26)+( CONSTS(34)*CONSTS(234)*ALGBRC(85)*STATES(23))/ALGBRC(83)+( CONSTS(34)*CONSTS(233)*STATES(32))/ALGBRC(83)+ ALGBRC(84)*STATES(31)*STATES(26)+( CONSTS(34)*CONSTS(234)*STATES(31)*STATES(23))/( ALGBRC(83)*CONSTS(227))+( CONSTS(34)*STATES(23)*STATES(32))/ALGBRC(83)+( ALGBRC(84)*CONSTS(231)*STATES(26)*STATES(32))/CONSTS(230)), 0.00000) ALGBRC(87) = TERNRY(STATES(23).GT.CONSTS(236).AND.STATES(32).GT.CONSTS(236), ALGBRC(86) - (( CONSTS(34)*ALGBRC(84)*STATES(23)*STATES(32))/ALGBRC(83))/( ALGBRC(84)*CONSTS(232)*STATES(31)+ ALGBRC(84)*CONSTS(231)*ALGBRC(85)*STATES(26)+( CONSTS(34)*CONSTS(234)*ALGBRC(85)*STATES(23))/ALGBRC(83)+( CONSTS(34)*CONSTS(233)*STATES(32))/ALGBRC(83)+ ALGBRC(84)*STATES(31)*STATES(26)+( CONSTS(34)*CONSTS(234)*STATES(31)*STATES(23))/( ALGBRC(83)*CONSTS(227))+( CONSTS(34)*STATES(23)*STATES(32))/ALGBRC(83)+( ALGBRC(84)*CONSTS(231)*STATES(26)*STATES(32))/CONSTS(230)), ALGBRC(86)) RATES(26) = (((ALGBRC(116) - ALGBRC(117)) - ALGBRC(87))+ALGBRC(88))/CONSTS(4) RATES(24) = (ALGBRC(117) - ALGBRC(67))/CONSTS(4) RATES(19) = ((- ALGBRC(115) - ALGBRC(117))+ALGBRC(67)+ALGBRC(57))/CONSTS(4) ALGBRC(45) = 1.00000+STATES(2)/CONSTS(371)+STATES(3)/CONSTS(372)+STATES(4)/CONSTS(373) ALGBRC(75) = ( CONSTS(322)*ALGBRC(37)*CONSTS(278))/( STATES(2)*ALGBRC(45)*CONSTS(277)) ALGBRC(76) = ( CONSTS(32)*CONSTS(209)*CONSTS(212))/( ALGBRC(75)*CONSTS(206)*CONSTS(211)) ALGBRC(77) = 1.00000+ALGBRC(19)/CONSTS(214)+ALGBRC(24)/CONSTS(215)+ALGBRC(29)/CONSTS(216) ALGBRC(118) = ( CONSTS(32)*ALGBRC(76)*( ALGBRC(114)*STATES(30) - ( STATES(23)*STATES(22))/ALGBRC(75)))/( ALGBRC(76)*CONSTS(206)*CONSTS(211)*ALGBRC(77)+ ALGBRC(76)*CONSTS(211)*ALGBRC(114)+ ALGBRC(76)*CONSTS(210)*ALGBRC(77)*STATES(30)+( CONSTS(32)*CONSTS(213)*ALGBRC(77)*STATES(23))/ALGBRC(75)+( CONSTS(32)*CONSTS(212)*STATES(22))/ALGBRC(75)+ ALGBRC(76)*ALGBRC(114)*STATES(30)+( CONSTS(32)*CONSTS(213)*ALGBRC(114)*STATES(23))/( ALGBRC(75)*CONSTS(206))+( CONSTS(32)*STATES(23)*STATES(22))/ALGBRC(75)+( ALGBRC(76)*CONSTS(210)*STATES(30)*STATES(22))/CONSTS(209)+( ALGBRC(76)*ALGBRC(114)*STATES(30)*STATES(23))/CONSTS(208)+( CONSTS(32)*STATES(30)*STATES(23)*STATES(22))/( CONSTS(207)*ALGBRC(75))) RATES(23) = (- ALGBRC(57)+ALGBRC(118)+ALGBRC(87))/CONSTS(4) ALGBRC(119) = CONSTS(16) - STATES(28) ALGBRC(44) = 1.00000+STATES(2)/CONSTS(370) ALGBRC(68) = ( CONSTS(335)*ALGBRC(44)*CONSTS(305))/( ALGBRC(43)*CONSTS(321)) ALGBRC(69) = ( CONSTS(30)*CONSTS(194)*CONSTS(191))/( ALGBRC(68)*CONSTS(190)*CONSTS(193)) ALGBRC(70) = (1.00000+STATES(23)/CONSTS(196)+STATES(27)/CONSTS(197)+STATES(29)/CONSTS(198))/(1.00000+STATES(27)/CONSTS(197)+STATES(29)/CONSTS(198)) ALGBRC(120) = ( CONSTS(30)*ALGBRC(69)*( STATES(27)*ALGBRC(119) - ( STATES(28)*STATES(29))/ALGBRC(68)))/( ALGBRC(69)*CONSTS(190)*CONSTS(193)*ALGBRC(70)+ ALGBRC(69)*CONSTS(193)*STATES(27)+ ALGBRC(69)*CONSTS(192)*ALGBRC(70)*ALGBRC(119)+( CONSTS(30)*CONSTS(195)*ALGBRC(70)*STATES(28))/ALGBRC(68)+( CONSTS(30)*CONSTS(194)*STATES(29))/ALGBRC(68)+ ALGBRC(69)*STATES(27)*ALGBRC(119)+( CONSTS(30)*CONSTS(195)*STATES(27)*STATES(28))/( ALGBRC(68)*CONSTS(190))+( ALGBRC(69)*CONSTS(192)*ALGBRC(119)*STATES(29))/CONSTS(191)+( CONSTS(30)*STATES(28)*STATES(29))/ALGBRC(68)) ALGBRC(71) = ( CONSTS(343)*ALGBRC(45))/ALGBRC(44) ALGBRC(72) = ( CONSTS(31)*CONSTS(200))/( ALGBRC(71)*CONSTS(199)) ALGBRC(31) = ( STATES(12)*(1.00000+STATES(2)/CONSTS(314)))/ALGBRC(30) ALGBRC(33) = ( STATES(13)*(1.00000+STATES(2)/CONSTS(338)))/ALGBRC(32) ALGBRC(73) = 1.00000+STATES(16)/CONSTS(201)+ALGBRC(19)/CONSTS(202)+ALGBRC(24)/CONSTS(203)+ALGBRC(31)/CONSTS(204)+ALGBRC(33)/CONSTS(205) ALGBRC(74) = ( ALGBRC(72)*CONSTS(31)*(STATES(29) - STATES(30)/ALGBRC(71)))/( CONSTS(199)*ALGBRC(72)*ALGBRC(73)+ ALGBRC(72)*STATES(29)+( CONSTS(31)*STATES(30))/ALGBRC(71)) RATES(29) = (ALGBRC(120) - ALGBRC(74))/CONSTS(4) ALGBRC(133) = STATES(6) ALGBRC(101) = 1.00000*EXP(- (CONSTS(288)+ 4.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(139) = ( ALGBRC(101)*STATES(2) ** 5.00000)/ALGBRC(133) ** 4.00000 ALGBRC(140) = CONSTS(42)*( ALGBRC(139)*STATES(22)*ALGBRC(119) - ALGBRC(114)*STATES(28)) RATES(22) = ((ALGBRC(115)+ALGBRC(116)+ALGBRC(117)+ALGBRC(118)) - ALGBRC(140))/CONSTS(4) ALGBRC(103) = 1.00000*EXP(- (CONSTS(289)+ 2.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(105) = ( ALGBRC(103)*STATES(2) ** 2.00000)/STATES(2) ** 4.00000 ALGBRC(107) = TERNRY(STATES(28).GT.CONSTS(248), STATES(28), CONSTS(248)) ALGBRC(121) = TERNRY(ALGBRC(119).GT.CONSTS(248), ALGBRC(119), CONSTS(248)) ALGBRC(122) = CONSTS(17) - STATES(56) ALGBRC(123) = (( CONSTS(43)*(1.00000+STATES(14)/CONSTS(44)))/(1.00000+STATES(14)/CONSTS(45)))*( ALGBRC(105) ** (1.0 / 2)*ALGBRC(122)* ALGBRC(107) ** (1.0 / 2) - STATES(56)* ALGBRC(121) ** (1.0 / 2)) RATES(28) = ((ALGBRC(120)+ALGBRC(140)) - ALGBRC(123))/CONSTS(4) ALGBRC(109) = EXP(- (CONSTS(290)+ 4.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(141) = ( ALGBRC(109)*STATES(2) ** 4.00000)/( ALGBRC(133) ** 2.00000*1.00000 ** 2.00000) ALGBRC(111) = TERNRY(STATES(57).GT.CONSTS(249), STATES(57), CONSTS(249)) ALGBRC(142) = CONSTS(46)*( (ALGBRC(111)/(ALGBRC(111)+CONSTS(14)))*EXP(( CONSTS(3)*STATES(38))/CONSTS(275))*(STATES(56)/CONSTS(17)))*( ALGBRC(141) ** 1.00000/2.00000*STATES(56)*ALGBRC(111)/1.00000 ** 1.00000/4.00000 - ALGBRC(122)) RATES(56) = ( 2.00000*ALGBRC(123) - 2.00000*ALGBRC(142))/CONSTS(5) ALGBRC(150) = STATES(7) ALGBRC(153) = 1.00000+ALGBRC(150)/CONSTS(375) ALGBRC(156) = CONSTS(35)*( (STATES(33)/ALGBRC(153))*ALGBRC(133) - (STATES(18)/ALGBRC(47))*STATES(2)) RATES(18) = (- ALGBRC(115)+ALGBRC(156))/CONSTS(4) RATES(33) = (- ALGBRC(156)+ALGBRC(89))/CONSTS(5) ALGBRC(154) = 1.00000+ALGBRC(133)/CONSTS(377)+ALGBRC(150)/CONSTS(378) ALGBRC(157) = CONSTS(36)*( (STATES(34)/ALGBRC(154))*ALGBRC(133) - (STATES(32)/ALGBRC(49))*STATES(2)) ALGBRC(155) = 1.00000+ALGBRC(133)/CONSTS(379)+ALGBRC(150)/CONSTS(380) ALGBRC(158) = ( EXP(( - CONSTS(3)*STATES(38))/CONSTS(275))*ALGBRC(50)*ALGBRC(154))/( ALGBRC(155)*ALGBRC(49)) ALGBRC(159) = ( CONSTS(41)*( ALGBRC(158)*STATES(37)*STATES(32)*STATES(2) - STATES(31)*STATES(34)*ALGBRC(133)))/( ALGBRC(158)*CONSTS(241)*CONSTS(244)*CONSTS(286)*( 2.00000*CONSTS(245)+( CONSTS(245)*STATES(37))/CONSTS(241)+( STATES(37)*STATES(32)*STATES(2))/( CONSTS(241)*CONSTS(244)*CONSTS(286))+( CONSTS(245)*STATES(31)*ALGBRC(133))/( CONSTS(242)*CONSTS(286))+( STATES(31)*STATES(34)*ALGBRC(133))/( CONSTS(242)*CONSTS(243)*CONSTS(286))+( CONSTS(245)*STATES(31))/CONSTS(242)+( CONSTS(245)*STATES(37)*STATES(2))/( CONSTS(241)*CONSTS(286))+( CONSTS(245)*STATES(2))/CONSTS(286)+( CONSTS(245)*STATES(34)*ALGBRC(133))/( CONSTS(243)*CONSTS(286))+( CONSTS(245)*ALGBRC(133))/CONSTS(286)+( CONSTS(245)*STATES(32)*STATES(2))/( CONSTS(244)*CONSTS(286)))) RATES(32) = ((ALGBRC(87)+ALGBRC(157)) - ALGBRC(159))/CONSTS(4) RATES(31) = (- ALGBRC(87)+ALGBRC(159))/CONSTS(4) RATES(34) = (- ALGBRC(157)+ALGBRC(159)+ALGBRC(94))/CONSTS(5) RATES(37) = (- ALGBRC(159)+ALGBRC(95))/CONSTS(5) ALGBRC(40) = (( STATES(16)*STATES(2))/CONSTS(361))/ALGBRC(38) ALGBRC(163) = STATES(8) ALGBRC(171) = 1.00000+ALGBRC(133)/CONSTS(361)+ALGBRC(150)/CONSTS(362)+ALGBRC(163)/CONSTS(363) ALGBRC(173) = (( STATES(17)*ALGBRC(133))/CONSTS(361))/ALGBRC(171) ALGBRC(175) = 1.00000+ALGBRC(133)/CONSTS(371)+ALGBRC(150)/CONSTS(372)+ALGBRC(163)/CONSTS(373) ALGBRC(176) = CONSTS(37)*(( ALGBRC(173)*STATES(30))/ALGBRC(45) - ( ALGBRC(40)*STATES(35))/ALGBRC(175)) RATES(16) = (- ALGBRC(60)+ALGBRC(176)+ALGBRC(57))/CONSTS(4) RATES(17) = (- ALGBRC(176)+ALGBRC(91))/CONSTS(5) ALGBRC(174) = 1.00000+ALGBRC(133)/CONSTS(367)+ALGBRC(150)/CONSTS(368)+ALGBRC(163)/CONSTS(369) ALGBRC(178) = CONSTS(39)*(( (STATES(39)/ALGBRC(174))*STATES(30))/ALGBRC(45) - ( (STATES(27)/ALGBRC(43))*STATES(35))/ALGBRC(175)) RATES(27) = ((ALGBRC(67) - ALGBRC(120))+ALGBRC(178))/CONSTS(4) ALGBRC(169) = 1.00000+ALGBRC(133)/CONSTS(352)+ALGBRC(150)/CONSTS(354)+ALGBRC(163)/CONSTS(356) ALGBRC(177) = CONSTS(40)*(( (STATES(35)/ALGBRC(175))*STATES(14))/ALGBRC(34) - ( (STATES(30)/ALGBRC(45))*STATES(15))/ALGBRC(169)) RATES(30) = (((((ALGBRC(74) - ALGBRC(118))+ALGBRC(177)) - ALGBRC(88)) - ALGBRC(176)) - ALGBRC(178))/CONSTS(4) RATES(39) = (ALGBRC(93) - ALGBRC(178))/CONSTS(5) RATES(35) = (- ALGBRC(177)+ALGBRC(92)+ALGBRC(88)+ALGBRC(176)+ALGBRC(178))/CONSTS(5) ALGBRC(113) = EXP(- (CONSTS(291) - CONSTS(8)*CONSTS(3)*STATES(38))/CONSTS(275))/1.00000 ALGBRC(143) = ( (( ALGBRC(113)*ALGBRC(133) ** CONSTS(8))/STATES(2) ** CONSTS(8) - 1.00000)*ALGBRC(1))/( ALGBRC(22)*ALGBRC(34)) ALGBRC(144) = CONSTS(47)*( ALGBRC(143)*STATES(9)*STATES(14) - STATES(1)) ALGBRC(164) = 1.00000+ALGBRC(133)/CONSTS(297)+ALGBRC(150)/CONSTS(312)+ALGBRC(163)/CONSTS(326) ALGBRC(166) = 1.00000+ALGBRC(133)/CONSTS(333)+ALGBRC(150)/CONSTS(337)+ALGBRC(163)/CONSTS(341) ALGBRC(179) = TERNRY(STATES(51).GT.CONSTS(250).AND.STATES(1).GT.CONSTS(250), CONSTS(48)*((STATES(51)/ALGBRC(166))/(STATES(51)/ALGBRC(166)+ (STATES(52)/ALGBRC(164))*EXP(( CONSTS(13)*CONSTS(3)*STATES(38))/CONSTS(275))) - (STATES(9)/ALGBRC(22))/(STATES(9)/ALGBRC(22)+ (STATES(1)/ALGBRC(1))*EXP(( CONSTS(13)*(CONSTS(3) - 1.00000)*STATES(38))/CONSTS(275))))*((STATES(51)/ALGBRC(166))/(STATES(51)/ALGBRC(166)+CONSTS(12))), 0.00000) ALGBRC(145) = TERNRY(ABS(STATES(38)).GT.CONSTS(251), ( CONSTS(52)*STATES(38)*( ALGBRC(133)*EXP(( CONSTS(3)*STATES(38))/CONSTS(275)) - STATES(2)))/(EXP(( CONSTS(3)*STATES(38))/CONSTS(275)) - 1.00000), ( CONSTS(52)*CONSTS(275)*(ALGBRC(133) - STATES(2)))/CONSTS(3)) RATES(38) = ((((( 4.00000*ALGBRC(140)+ 2.00000*ALGBRC(123)+ 4.00000*ALGBRC(142)) - CONSTS(8)*ALGBRC(144)) - ALGBRC(179)) - ALGBRC(145))+ALGBRC(159))/CONSTS(22) RATES(1) = ((ALGBRC(82)+ALGBRC(144)) - ALGBRC(179))/CONSTS(4) RATES(9) = ((- ALGBRC(82) - ALGBRC(144))+ALGBRC(179))/CONSTS(4) ALGBRC(35) = (( STATES(14)*STATES(2))/CONSTS(352))/ALGBRC(34) ALGBRC(170) = (( STATES(15)*ALGBRC(133))/CONSTS(352))/ALGBRC(169) ALGBRC(180) = ( CONSTS(49)*( ALGBRC(170)*ALGBRC(133) - ALGBRC(35)*STATES(2)))/( CONSTS(50)*(1.00000+ALGBRC(170)/CONSTS(50))*(1.00000+ALGBRC(35)/CONSTS(50))) RATES(14) = (((- ALGBRC(67) - ALGBRC(144))+ALGBRC(180)) - ALGBRC(177))/CONSTS(4) RATES(15) = (- ALGBRC(180)+ALGBRC(99)+ALGBRC(177))/CONSTS(5) ALGBRC(168) = 1.00000+ALGBRC(133)/CONSTS(345)+ALGBRC(150)/CONSTS(348)+ALGBRC(163)/CONSTS(350) ALGBRC(181) = ( CONSTS(382)*ALGBRC(164)*ALGBRC(168))/( ALGBRC(166)*ALGBRC(166)) ALGBRC(183) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(262)*( ALGBRC(181)*STATES(51)*STATES(51) - STATES(54)*STATES(52))) RATES(52) = (ALGBRC(97)+ALGBRC(179)+ALGBRC(183))/CONSTS(5) RATES(51) = ((ALGBRC(96) - ALGBRC(179)) - 2.00000*ALGBRC(183))/CONSTS(5) RATES(54) = (ALGBRC(98)+ALGBRC(183))/CONSTS(5) ALGBRC(190) = - (( STATES(2)*STATES(2))/( CONSTS(135)*CONSTS(136)*CONSTS(283) ** 2.00000)+( STATES(2)*STATES(1))/( CONSTS(297)*CONSTS(312)*ALGBRC(1) ** 2.00000)+( STATES(2)*STATES(9))/( CONSTS(333)*CONSTS(337)*ALGBRC(22) ** 2.00000)+( STATES(2)*STATES(11))/( CONSTS(345)*CONSTS(348)*ALGBRC(27) ** 2.00000)+( STATES(2)*STATES(12))/( CONSTS(297)*CONSTS(312)*ALGBRC(30) ** 2.00000)+( STATES(2)*STATES(13))/( CONSTS(333)*CONSTS(337)*ALGBRC(32) ** 2.00000)+( STATES(2)*STATES(14))/( CONSTS(352)*CONSTS(354)*ALGBRC(34) ** 2.00000)+( STATES(2)*STATES(22))/( CONSTS(111)*CONSTS(112)*CONSTS(278) ** 2.00000)+( STATES(2)*ALGBRC(114))/( CONSTS(108)*CONSTS(109)*CONSTS(277) ** 2.00000)+( STATES(2)*STATES(28))/( CONSTS(114)*CONSTS(115)*CONSTS(305) ** 2.00000)+( STATES(2)*ALGBRC(119))/( CONSTS(117)*CONSTS(118)*CONSTS(321) ** 2.00000)+( STATES(2)*STATES(18))/( CONSTS(99)*CONSTS(375)*ALGBRC(47) ** 2.00000)+( STATES(2)*STATES(23))/( CONSTS(90)*CONSTS(360)*ALGBRC(37) ** 2.00000)+( STATES(2)*STATES(21))/( CONSTS(120)*CONSTS(121)*CONSTS(279) ** 2.00000)+( STATES(2)*STATES(16))/( CONSTS(361)*CONSTS(362)*ALGBRC(38) ** 2.00000)+( STATES(2)*STATES(25))/( CONSTS(364)*CONSTS(365)*ALGBRC(41) ** 2.00000)+( STATES(2)*STATES(26))/( CONSTS(123)*CONSTS(124)*CONSTS(307) ** 2.00000)+( STATES(2)*STATES(24))/( CONSTS(366)*CONSTS(93)*ALGBRC(42) ** 2.00000)+( STATES(2)*STATES(19))/( CONSTS(358)*CONSTS(88)*ALGBRC(36) ** 2.00000)+( STATES(2)*STATES(27))/( CONSTS(367)*CONSTS(368)*ALGBRC(43) ** 2.00000)+( STATES(2)*STATES(29))/( CONSTS(370)*CONSTS(95)*ALGBRC(44) ** 2.00000)+( STATES(2)*STATES(30))/( CONSTS(371)*CONSTS(372)*ALGBRC(45) ** 2.00000)+( STATES(2)*STATES(32))/( CONSTS(377)*CONSTS(378)*ALGBRC(49) ** 2.00000)+( STATES(2)*STATES(31))/( CONSTS(379)*CONSTS(380)*ALGBRC(50) ** 2.00000)+( STATES(2)*STATES(4))/( CONSTS(141)*CONSTS(142)*CONSTS(285) ** 2.00000)+( STATES(2)*STATES(3))/( CONSTS(138)*CONSTS(139)*CONSTS(284) ** 2.00000)+( STATES(2)*STATES(57))/( CONSTS(126)*CONSTS(127)*CONSTS(280) ** 2.00000)+( STATES(2)*CONSTS(292))/( CONSTS(129)*CONSTS(130)*CONSTS(281) ** 2.00000)+( STATES(2)*CONSTS(311))/( CONSTS(132)*CONSTS(133)*CONSTS(282) ** 2.00000)+( STATES(2)*STATES(20))/( CONSTS(374)*CONSTS(97)*ALGBRC(46) ** 2.00000)) ALGBRC(191) = - (( STATES(2)*STATES(2))/( CONSTS(135)*CONSTS(137)*CONSTS(283) ** 2.00000)+( STATES(2)*STATES(1))/( CONSTS(297)*CONSTS(326)*ALGBRC(1) ** 2.00000)+( STATES(2)*STATES(9))/( CONSTS(333)*CONSTS(341)*ALGBRC(22) ** 2.00000)+( STATES(2)*STATES(11))/( CONSTS(345)*CONSTS(350)*ALGBRC(27) ** 2.00000)+( STATES(2)*STATES(12))/( CONSTS(297)*CONSTS(326)*ALGBRC(30) ** 2.00000)+( STATES(2)*STATES(13))/( CONSTS(333)*CONSTS(341)*ALGBRC(32) ** 2.00000)+( STATES(2)*STATES(14))/( CONSTS(352)*CONSTS(356)*ALGBRC(34) ** 2.00000)+( STATES(2)*STATES(22))/( CONSTS(111)*CONSTS(113)*CONSTS(278) ** 2.00000)+( STATES(2)*ALGBRC(114))/( CONSTS(108)*CONSTS(110)*CONSTS(277) ** 2.00000)+( STATES(2)*STATES(28))/( CONSTS(114)*CONSTS(116)*CONSTS(305) ** 2.00000)+( STATES(2)*ALGBRC(119))/( CONSTS(117)*CONSTS(119)*CONSTS(321) ** 2.00000)+( STATES(2)*STATES(18))/( CONSTS(99)*CONSTS(100)*ALGBRC(47) ** 2.00000)+( STATES(2)*STATES(23))/( CONSTS(90)*CONSTS(91)*ALGBRC(37) ** 2.00000)+( STATES(2)*STATES(21))/( CONSTS(120)*CONSTS(122)*CONSTS(279) ** 2.00000)+( STATES(2)*STATES(16))/( CONSTS(361)*CONSTS(363)*ALGBRC(38) ** 2.00000)+( STATES(2)*STATES(25))/( CONSTS(364)*CONSTS(92)*ALGBRC(41) ** 2.00000)+( STATES(2)*STATES(26))/( CONSTS(123)*CONSTS(125)*CONSTS(307) ** 2.00000)+( STATES(2)*STATES(24))/( CONSTS(366)*CONSTS(94)*ALGBRC(42) ** 2.00000)+( STATES(2)*STATES(19))/( CONSTS(358)*CONSTS(89)*ALGBRC(36) ** 2.00000)+( STATES(2)*STATES(27))/( CONSTS(367)*CONSTS(369)*ALGBRC(43) ** 2.00000)+( STATES(2)*STATES(29))/( CONSTS(370)*CONSTS(96)*ALGBRC(44) ** 2.00000)+( STATES(2)*STATES(30))/( CONSTS(371)*CONSTS(373)*ALGBRC(45) ** 2.00000)+( STATES(2)*STATES(32))/( CONSTS(377)*CONSTS(103)*ALGBRC(49) ** 2.00000)+( STATES(2)*STATES(31))/( CONSTS(379)*CONSTS(104)*ALGBRC(50) ** 2.00000)+( STATES(2)*STATES(4))/( CONSTS(141)*CONSTS(143)*CONSTS(285) ** 2.00000)+( STATES(2)*STATES(3))/( CONSTS(138)*CONSTS(140)*CONSTS(284) ** 2.00000)+( STATES(2)*STATES(57))/( CONSTS(126)*CONSTS(128)*CONSTS(280) ** 2.00000)+( STATES(2)*CONSTS(292))/( CONSTS(129)*CONSTS(131)*CONSTS(281) ** 2.00000)+( STATES(2)*CONSTS(311))/( CONSTS(132)*CONSTS(134)*CONSTS(282) ** 2.00000)+( STATES(2)*STATES(20))/( CONSTS(374)*CONSTS(98)*ALGBRC(46) ** 2.00000)) ALGBRC(194) = - (( STATES(3)*STATES(2))/( CONSTS(136)*CONSTS(137)*CONSTS(283) ** 2.00000)+( STATES(3)*STATES(1))/( CONSTS(312)*CONSTS(326)*ALGBRC(1) ** 2.00000)+( STATES(3)*STATES(9))/( CONSTS(337)*CONSTS(341)*ALGBRC(22) ** 2.00000)+( STATES(3)*STATES(11))/( CONSTS(348)*CONSTS(350)*ALGBRC(27) ** 2.00000)+( STATES(3)*STATES(12))/( CONSTS(312)*CONSTS(326)*ALGBRC(30) ** 2.00000)+( STATES(3)*STATES(13))/( CONSTS(337)*CONSTS(341)*ALGBRC(32) ** 2.00000)+( STATES(3)*STATES(14))/( CONSTS(354)*CONSTS(356)*ALGBRC(34) ** 2.00000)+( STATES(3)*STATES(22))/( CONSTS(112)*CONSTS(113)*CONSTS(278) ** 2.00000)+( STATES(3)*ALGBRC(114))/( CONSTS(109)*CONSTS(110)*CONSTS(277) ** 2.00000)+( STATES(3)*STATES(28))/( CONSTS(115)*CONSTS(116)*CONSTS(305) ** 2.00000)+( STATES(3)*ALGBRC(119))/( CONSTS(118)*CONSTS(119)*CONSTS(321) ** 2.00000)+( STATES(3)*STATES(18))/( CONSTS(375)*CONSTS(100)*ALGBRC(47) ** 2.00000)+( STATES(3)*STATES(23))/( CONSTS(360)*CONSTS(91)*ALGBRC(37) ** 2.00000)+( STATES(3)*STATES(21))/( CONSTS(121)*CONSTS(122)*CONSTS(279) ** 2.00000)+( STATES(3)*STATES(16))/( CONSTS(362)*CONSTS(363)*ALGBRC(38) ** 2.00000)+( STATES(3)*STATES(25))/( CONSTS(365)*CONSTS(92)*ALGBRC(41) ** 2.00000)+( STATES(3)*STATES(26))/( CONSTS(124)*CONSTS(125)*CONSTS(307) ** 2.00000)+( STATES(3)*STATES(24))/( CONSTS(93)*CONSTS(94)*ALGBRC(42) ** 2.00000)+( STATES(3)*STATES(19))/( CONSTS(88)*CONSTS(89)*ALGBRC(36) ** 2.00000)+( STATES(3)*STATES(27))/( CONSTS(368)*CONSTS(369)*ALGBRC(43) ** 2.00000)+( STATES(3)*STATES(29))/( CONSTS(95)*CONSTS(96)*ALGBRC(44) ** 2.00000)+( STATES(3)*STATES(30))/( CONSTS(372)*CONSTS(373)*ALGBRC(45) ** 2.00000)+( STATES(3)*STATES(32))/( CONSTS(378)*CONSTS(103)*ALGBRC(49) ** 2.00000)+( STATES(3)*STATES(31))/( CONSTS(380)*CONSTS(104)*ALGBRC(50) ** 2.00000)+( STATES(3)*STATES(4))/( CONSTS(142)*CONSTS(143)*CONSTS(285) ** 2.00000)+( STATES(3)*STATES(3))/( CONSTS(139)*CONSTS(140)*CONSTS(284) ** 2.00000)+( STATES(3)*STATES(57))/( CONSTS(127)*CONSTS(128)*CONSTS(280) ** 2.00000)+( STATES(3)*CONSTS(292))/( CONSTS(130)*CONSTS(131)*CONSTS(281) ** 2.00000)+( STATES(3)*CONSTS(311))/( CONSTS(133)*CONSTS(134)*CONSTS(282) ** 2.00000)+( STATES(3)*STATES(20))/( CONSTS(97)*CONSTS(98)*ALGBRC(46) ** 2.00000)) ALGBRC(197) = - (( STATES(4)*STATES(2))/( CONSTS(137)*CONSTS(136)*CONSTS(283) ** 2.00000)+( STATES(4)*STATES(1))/( CONSTS(326)*CONSTS(312)*ALGBRC(1) ** 2.00000)+( STATES(4)*STATES(9))/( CONSTS(341)*CONSTS(337)*ALGBRC(22) ** 2.00000)+( STATES(4)*STATES(11))/( CONSTS(350)*CONSTS(348)*ALGBRC(27) ** 2.00000)+( STATES(4)*STATES(12))/( CONSTS(326)*CONSTS(312)*ALGBRC(30) ** 2.00000)+( STATES(4)*STATES(13))/( CONSTS(341)*CONSTS(337)*ALGBRC(32) ** 2.00000)+( STATES(4)*STATES(14))/( CONSTS(356)*CONSTS(354)*ALGBRC(34) ** 2.00000)+( STATES(4)*STATES(22))/( CONSTS(113)*CONSTS(112)*CONSTS(278) ** 2.00000)+( STATES(4)*ALGBRC(114))/( CONSTS(110)*CONSTS(109)*CONSTS(277) ** 2.00000)+( STATES(4)*STATES(28))/( CONSTS(116)*CONSTS(115)*CONSTS(305) ** 2.00000)+( STATES(4)*ALGBRC(119))/( CONSTS(119)*CONSTS(118)*CONSTS(321) ** 2.00000)+( STATES(4)*STATES(18))/( CONSTS(100)*CONSTS(375)*ALGBRC(47) ** 2.00000)+( STATES(4)*STATES(23))/( CONSTS(91)*CONSTS(360)*ALGBRC(37) ** 2.00000)+( STATES(4)*STATES(21))/( CONSTS(122)*CONSTS(121)*CONSTS(279) ** 2.00000)+( STATES(4)*STATES(16))/( CONSTS(363)*CONSTS(362)*ALGBRC(38) ** 2.00000)+( STATES(4)*STATES(25))/( CONSTS(92)*CONSTS(365)*ALGBRC(41) ** 2.00000)+( STATES(4)*STATES(26))/( CONSTS(125)*CONSTS(124)*CONSTS(307) ** 2.00000)+( STATES(4)*STATES(24))/( CONSTS(94)*CONSTS(93)*ALGBRC(42) ** 2.00000)+( STATES(4)*STATES(19))/( CONSTS(89)*CONSTS(88)*ALGBRC(36) ** 2.00000)+( STATES(4)*STATES(27))/( CONSTS(369)*CONSTS(368)*ALGBRC(43) ** 2.00000)+( STATES(4)*STATES(29))/( CONSTS(96)*CONSTS(95)*ALGBRC(44) ** 2.00000)+( STATES(4)*STATES(30))/( CONSTS(373)*CONSTS(372)*ALGBRC(45) ** 2.00000)+( STATES(4)*STATES(32))/( CONSTS(103)*CONSTS(378)*ALGBRC(49) ** 2.00000)+( STATES(4)*STATES(31))/( CONSTS(104)*CONSTS(380)*ALGBRC(50) ** 2.00000)+( STATES(4)*STATES(4))/( CONSTS(143)*CONSTS(142)*CONSTS(285) ** 2.00000)+( STATES(4)*STATES(3))/( CONSTS(140)*CONSTS(139)*CONSTS(284) ** 2.00000)+( STATES(4)*STATES(57))/( CONSTS(128)*CONSTS(127)*CONSTS(280) ** 2.00000)+( STATES(4)*CONSTS(292))/( CONSTS(131)*CONSTS(130)*CONSTS(281) ** 2.00000)+( STATES(4)*CONSTS(311))/( CONSTS(134)*CONSTS(133)*CONSTS(282) ** 2.00000)+( STATES(4)*STATES(20))/( CONSTS(98)*CONSTS(97)*ALGBRC(46) ** 2.00000)) ALGBRC(198) = ( (1.00000+STATES(2)/CONSTS(135)+STATES(3)/CONSTS(136))*STATES(2))/( CONSTS(137)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312))*STATES(1))/( CONSTS(326)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337))*STATES(9))/( CONSTS(341)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(345)+STATES(3)/CONSTS(348))*STATES(11))/( CONSTS(350)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312))*STATES(12))/( CONSTS(326)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337))*STATES(13))/( CONSTS(341)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(352)+STATES(3)/CONSTS(354))*STATES(14))/( CONSTS(356)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(111)+STATES(3)/CONSTS(112))*STATES(22))/( CONSTS(113)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(108)+STATES(3)/CONSTS(109))*ALGBRC(114))/( CONSTS(110)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(114)+STATES(3)/CONSTS(115))*STATES(28))/( CONSTS(116)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(117)+STATES(3)/CONSTS(118))*ALGBRC(119))/( CONSTS(119)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(99)+STATES(3)/CONSTS(375))*STATES(18))/( CONSTS(100)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(90)+STATES(3)/CONSTS(360))*STATES(23))/( CONSTS(91)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(120)+STATES(3)/CONSTS(121))*STATES(21))/( CONSTS(122)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(361)+STATES(3)/CONSTS(362))*STATES(16))/( CONSTS(363)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(364)+STATES(3)/CONSTS(365))*STATES(25))/( CONSTS(92)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(123)+STATES(3)/CONSTS(124))*STATES(26))/( CONSTS(125)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(366)+STATES(3)/CONSTS(93))*STATES(24))/( CONSTS(94)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(358)+STATES(3)/CONSTS(88))*STATES(19))/( CONSTS(89)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(367)+STATES(3)/CONSTS(368))*STATES(27))/( CONSTS(369)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(370)+STATES(3)/CONSTS(95))*STATES(29))/( CONSTS(96)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(371)+STATES(3)/CONSTS(372))*STATES(30))/( CONSTS(373)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(377)+STATES(3)/CONSTS(378))*STATES(32))/( CONSTS(103)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(379)+STATES(3)/CONSTS(380))*STATES(31))/( CONSTS(104)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(141)+STATES(3)/CONSTS(142))*STATES(4))/( CONSTS(143)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(138)+STATES(3)/CONSTS(139))*STATES(3))/( CONSTS(140)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(126)+STATES(3)/CONSTS(127))*STATES(57))/( CONSTS(128)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(129)+STATES(3)/CONSTS(130))*CONSTS(292))/( CONSTS(131)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(132)+STATES(3)/CONSTS(133))*CONSTS(311))/( CONSTS(134)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(374)+STATES(3)/CONSTS(97))*STATES(20))/( CONSTS(98)*ALGBRC(46) ** 2.00000) ALGBRC(206) = 1.00000+ALGBRC(198) ALGBRC(195) = ( (1.00000+STATES(2)/CONSTS(135)+STATES(4)/CONSTS(137))*STATES(2))/( CONSTS(136)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(4)/CONSTS(326))*STATES(1))/( CONSTS(312)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(4)/CONSTS(341))*STATES(9))/( CONSTS(337)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(345)+STATES(4)/CONSTS(350))*STATES(11))/( CONSTS(348)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(4)/CONSTS(326))*STATES(12))/( CONSTS(312)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(4)/CONSTS(341))*STATES(13))/( CONSTS(337)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(352)+STATES(4)/CONSTS(356))*STATES(14))/( CONSTS(354)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(111)+STATES(4)/CONSTS(113))*STATES(22))/( CONSTS(112)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(108)+STATES(4)/CONSTS(110))*ALGBRC(114))/( CONSTS(109)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(114)+STATES(4)/CONSTS(116))*STATES(28))/( CONSTS(115)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(117)+STATES(4)/CONSTS(119))*ALGBRC(119))/( CONSTS(118)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(99)+STATES(4)/CONSTS(100))*STATES(18))/( CONSTS(375)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(90)+STATES(4)/CONSTS(91))*STATES(23))/( CONSTS(360)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(120)+STATES(4)/CONSTS(122))*STATES(21))/( CONSTS(121)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(361)+STATES(4)/CONSTS(363))*STATES(16))/( CONSTS(362)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(364)+STATES(4)/CONSTS(92))*STATES(25))/( CONSTS(365)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(123)+STATES(4)/CONSTS(125))*STATES(26))/( CONSTS(124)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(366)+STATES(4)/CONSTS(94))*STATES(24))/( CONSTS(93)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(358)+STATES(4)/CONSTS(89))*STATES(19))/( CONSTS(88)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(367)+STATES(4)/CONSTS(369))*STATES(27))/( CONSTS(368)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(370)+STATES(4)/CONSTS(96))*STATES(29))/( CONSTS(95)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(371)+STATES(4)/CONSTS(373))*STATES(30))/( CONSTS(372)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(377)+STATES(4)/CONSTS(103))*STATES(32))/( CONSTS(378)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(379)+STATES(4)/CONSTS(104))*STATES(31))/( CONSTS(380)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(141)+STATES(4)/CONSTS(143))*STATES(4))/( CONSTS(142)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(138)+STATES(4)/CONSTS(140))*STATES(3))/( CONSTS(139)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(126)+STATES(4)/CONSTS(128))*STATES(57))/( CONSTS(127)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(129)+STATES(4)/CONSTS(131))*CONSTS(292))/( CONSTS(130)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(132)+STATES(4)/CONSTS(134))*CONSTS(311))/( CONSTS(133)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(374)+STATES(4)/CONSTS(98))*STATES(20))/( CONSTS(97)*ALGBRC(46) ** 2.00000) ALGBRC(205) = 1.00000+ALGBRC(195) ALGBRC(192) = ( (1.00000+STATES(3)/CONSTS(136)+STATES(4)/CONSTS(137))*STATES(2))/( CONSTS(135)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326))*STATES(1))/( CONSTS(297)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341))*STATES(9))/( CONSTS(333)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(348)+STATES(4)/CONSTS(350))*STATES(11))/( CONSTS(345)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326))*STATES(12))/( CONSTS(297)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341))*STATES(13))/( CONSTS(333)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(354)+STATES(4)/CONSTS(356))*STATES(14))/( CONSTS(352)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(112)+STATES(4)/CONSTS(113))*STATES(22))/( CONSTS(111)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(109)+STATES(4)/CONSTS(110))*ALGBRC(114))/( CONSTS(108)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(115)+STATES(4)/CONSTS(116))*STATES(28))/( CONSTS(114)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(118)+STATES(4)/CONSTS(119))*ALGBRC(119))/( CONSTS(117)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(375)+STATES(4)/CONSTS(100))*STATES(18))/( CONSTS(99)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(360)+STATES(4)/CONSTS(91))*STATES(23))/( CONSTS(90)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(121)+STATES(4)/CONSTS(122))*STATES(21))/( CONSTS(120)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(362)+STATES(4)/CONSTS(363))*STATES(16))/( CONSTS(361)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(365)+STATES(4)/CONSTS(92))*STATES(25))/( CONSTS(364)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(124)+STATES(4)/CONSTS(125))*STATES(26))/( CONSTS(123)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(93)+STATES(4)/CONSTS(94))*STATES(24))/( CONSTS(366)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(88)+STATES(4)/CONSTS(89))*STATES(19))/( CONSTS(358)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(368)+STATES(4)/CONSTS(369))*STATES(27))/( CONSTS(367)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(95)+STATES(4)/CONSTS(96))*STATES(29))/( CONSTS(370)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(372)+STATES(4)/CONSTS(373))*STATES(30))/( CONSTS(371)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(378)+STATES(4)/CONSTS(103))*STATES(32))/( CONSTS(377)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(380)+STATES(4)/CONSTS(104))*STATES(31))/( CONSTS(379)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(142)+STATES(4)/CONSTS(143))*STATES(4))/( CONSTS(141)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(139)+STATES(4)/CONSTS(140))*STATES(3))/( CONSTS(138)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(127)+STATES(4)/CONSTS(128))*STATES(57))/( CONSTS(126)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(130)+STATES(4)/CONSTS(131))*CONSTS(292))/( CONSTS(129)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(133)+STATES(4)/CONSTS(134))*CONSTS(311))/( CONSTS(132)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(97)+STATES(4)/CONSTS(98))*STATES(20))/( CONSTS(374)*ALGBRC(46) ** 2.00000) ALGBRC(204) = 1.00000+ALGBRC(192)+CONSTS(272)/( CONSTS(273)*1.00000+STATES(2)/CONSTS(273) ** 2.00000) ALGBRC(196) = - (( STATES(4)*STATES(2))/( CONSTS(137)*CONSTS(135)*CONSTS(283) ** 2.00000)+( STATES(4)*STATES(1))/( CONSTS(326)*CONSTS(297)*ALGBRC(1) ** 2.00000)+( STATES(4)*STATES(9))/( CONSTS(341)*CONSTS(333)*ALGBRC(22) ** 2.00000)+( STATES(4)*STATES(11))/( CONSTS(350)*CONSTS(345)*ALGBRC(27) ** 2.00000)+( STATES(4)*STATES(12))/( CONSTS(326)*CONSTS(297)*ALGBRC(30) ** 2.00000)+( STATES(4)*STATES(13))/( CONSTS(341)*CONSTS(333)*ALGBRC(32) ** 2.00000)+( STATES(4)*STATES(14))/( CONSTS(356)*CONSTS(352)*ALGBRC(34) ** 2.00000)+( STATES(4)*STATES(22))/( CONSTS(113)*CONSTS(111)*CONSTS(278) ** 2.00000)+( STATES(4)*ALGBRC(114))/( CONSTS(110)*CONSTS(108)*CONSTS(277) ** 2.00000)+( STATES(4)*STATES(28))/( CONSTS(116)*CONSTS(114)*CONSTS(305) ** 2.00000)+( STATES(4)*ALGBRC(119))/( CONSTS(119)*CONSTS(117)*CONSTS(321) ** 2.00000)+( STATES(4)*STATES(18))/( CONSTS(100)*CONSTS(99)*ALGBRC(47) ** 2.00000)+( STATES(4)*STATES(23))/( CONSTS(91)*CONSTS(90)*ALGBRC(37) ** 2.00000)+( STATES(4)*STATES(21))/( CONSTS(122)*CONSTS(120)*CONSTS(279) ** 2.00000)+( STATES(4)*STATES(16))/( CONSTS(363)*CONSTS(361)*ALGBRC(38) ** 2.00000)+( STATES(4)*STATES(25))/( CONSTS(92)*CONSTS(364)*ALGBRC(41) ** 2.00000)+( STATES(4)*STATES(26))/( CONSTS(125)*CONSTS(123)*CONSTS(307) ** 2.00000)+( STATES(4)*STATES(24))/( CONSTS(94)*CONSTS(366)*ALGBRC(42) ** 2.00000)+( STATES(4)*STATES(19))/( CONSTS(89)*CONSTS(358)*ALGBRC(36) ** 2.00000)+( STATES(4)*STATES(27))/( CONSTS(369)*CONSTS(367)*ALGBRC(43) ** 2.00000)+( STATES(4)*STATES(29))/( CONSTS(96)*CONSTS(370)*ALGBRC(44) ** 2.00000)+( STATES(4)*STATES(30))/( CONSTS(373)*CONSTS(371)*ALGBRC(45) ** 2.00000)+( STATES(4)*STATES(32))/( CONSTS(103)*CONSTS(377)*ALGBRC(49) ** 2.00000)+( STATES(4)*STATES(31))/( CONSTS(104)*CONSTS(379)*ALGBRC(50) ** 2.00000)+( STATES(4)*STATES(4))/( CONSTS(143)*CONSTS(141)*CONSTS(285) ** 2.00000)+( STATES(4)*STATES(3))/( CONSTS(140)*CONSTS(138)*CONSTS(284) ** 2.00000)+( STATES(4)*STATES(57))/( CONSTS(128)*CONSTS(126)*CONSTS(280) ** 2.00000)+( STATES(4)*CONSTS(292))/( CONSTS(131)*CONSTS(129)*CONSTS(281) ** 2.00000)+( STATES(4)*CONSTS(311))/( CONSTS(134)*CONSTS(132)*CONSTS(282) ** 2.00000)+( STATES(4)*STATES(20))/( CONSTS(98)*CONSTS(374)*ALGBRC(46) ** 2.00000)) ALGBRC(193) = - (( STATES(3)*STATES(2))/( CONSTS(136)*CONSTS(135)*CONSTS(283) ** 2.00000)+( STATES(3)*STATES(1))/( CONSTS(312)*CONSTS(297)*ALGBRC(1) ** 2.00000)+( STATES(3)*STATES(9))/( CONSTS(337)*CONSTS(333)*ALGBRC(22) ** 2.00000)+( STATES(3)*STATES(11))/( CONSTS(348)*CONSTS(345)*ALGBRC(27) ** 2.00000)+( STATES(3)*STATES(12))/( CONSTS(312)*CONSTS(297)*ALGBRC(30) ** 2.00000)+( STATES(3)*STATES(13))/( CONSTS(337)*CONSTS(333)*ALGBRC(32) ** 2.00000)+( STATES(3)*STATES(14))/( CONSTS(354)*CONSTS(352)*ALGBRC(34) ** 2.00000)+( STATES(3)*STATES(22))/( CONSTS(112)*CONSTS(111)*CONSTS(278) ** 2.00000)+( STATES(3)*ALGBRC(114))/( CONSTS(109)*CONSTS(108)*CONSTS(277) ** 2.00000)+( STATES(3)*STATES(28))/( CONSTS(115)*CONSTS(114)*CONSTS(305) ** 2.00000)+( STATES(3)*ALGBRC(119))/( CONSTS(118)*CONSTS(117)*CONSTS(321) ** 2.00000)+( STATES(3)*STATES(18))/( CONSTS(375)*CONSTS(99)*ALGBRC(47) ** 2.00000)+( STATES(3)*STATES(23))/( CONSTS(360)*CONSTS(90)*ALGBRC(37) ** 2.00000)+( STATES(3)*STATES(21))/( CONSTS(121)*CONSTS(120)*CONSTS(279) ** 2.00000)+( STATES(3)*STATES(16))/( CONSTS(362)*CONSTS(361)*ALGBRC(38) ** 2.00000)+( STATES(3)*STATES(25))/( CONSTS(365)*CONSTS(364)*ALGBRC(41) ** 2.00000)+( STATES(3)*STATES(26))/( CONSTS(124)*CONSTS(123)*CONSTS(307) ** 2.00000)+( STATES(3)*STATES(24))/( CONSTS(93)*CONSTS(366)*ALGBRC(42) ** 2.00000)+( STATES(3)*STATES(19))/( CONSTS(88)*CONSTS(358)*ALGBRC(36) ** 2.00000)+( STATES(3)*STATES(27))/( CONSTS(368)*CONSTS(367)*ALGBRC(43) ** 2.00000)+( STATES(3)*STATES(29))/( CONSTS(95)*CONSTS(370)*ALGBRC(44) ** 2.00000)+( STATES(3)*STATES(30))/( CONSTS(372)*CONSTS(371)*ALGBRC(45) ** 2.00000)+( STATES(3)*STATES(32))/( CONSTS(378)*CONSTS(377)*ALGBRC(49) ** 2.00000)+( STATES(3)*STATES(31))/( CONSTS(380)*CONSTS(379)*ALGBRC(50) ** 2.00000)+( STATES(3)*STATES(4))/( CONSTS(142)*CONSTS(141)*CONSTS(285) ** 2.00000)+( STATES(3)*STATES(3))/( CONSTS(139)*CONSTS(138)*CONSTS(284) ** 2.00000)+( STATES(3)*STATES(57))/( CONSTS(127)*CONSTS(126)*CONSTS(280) ** 2.00000)+( STATES(3)*CONSTS(292))/( CONSTS(130)*CONSTS(129)*CONSTS(281) ** 2.00000)+( STATES(3)*CONSTS(311))/( CONSTS(133)*CONSTS(132)*CONSTS(282) ** 2.00000)+( STATES(3)*STATES(20))/( CONSTS(97)*CONSTS(374)*ALGBRC(46) ** 2.00000)) ALGBRC(207) = ((( ALGBRC(204)*ALGBRC(197)*ALGBRC(194)+ ALGBRC(206)*ALGBRC(190)*ALGBRC(193)+ ALGBRC(205)*ALGBRC(191)*ALGBRC(196)) - ALGBRC(205)*ALGBRC(206)*ALGBRC(204)) - ALGBRC(191)*ALGBRC(197)*ALGBRC(193)) - ALGBRC(190)*ALGBRC(194)*ALGBRC(196) ALGBRC(146) = ((ALGBRC(115)+ALGBRC(116)+ALGBRC(117)+ALGBRC(118)) - ALGBRC(140))/CONSTS(4) ALGBRC(147) = - ALGBRC(146) ALGBRC(148) = ((ALGBRC(120)+ALGBRC(140)) - ALGBRC(123))/CONSTS(4) ALGBRC(149) = - ALGBRC(148) ALGBRC(184) = ((ALGBRC(82)+ALGBRC(144)) - ALGBRC(179))/CONSTS(4) ALGBRC(185) = ((- ALGBRC(82) - ALGBRC(144))+ALGBRC(179))/CONSTS(4) ALGBRC(124) = (ALGBRC(67) - ALGBRC(82))/CONSTS(4) ALGBRC(125) = (- ALGBRC(67)+ALGBRC(82))/CONSTS(4) ALGBRC(186) = (((- ALGBRC(67) - ALGBRC(144))+ALGBRC(180)) - ALGBRC(177))/CONSTS(4) ALGBRC(160) = (- ALGBRC(115)+ALGBRC(156))/CONSTS(4) ALGBRC(126) = (ALGBRC(115) - ALGBRC(57))/CONSTS(4) ALGBRC(187) = (- ALGBRC(60)+ALGBRC(176)+ALGBRC(57))/CONSTS(4) ALGBRC(127) = (ALGBRC(60) - ALGBRC(116))/CONSTS(4) ALGBRC(128) = (((ALGBRC(116) - ALGBRC(117)) - ALGBRC(87))+ALGBRC(88))/CONSTS(4) ALGBRC(129) = (ALGBRC(117) - ALGBRC(67))/CONSTS(4) ALGBRC(130) = ((- ALGBRC(115) - ALGBRC(117))+ALGBRC(67)+ALGBRC(57))/CONSTS(4) ALGBRC(188) = ((ALGBRC(67) - ALGBRC(120))+ALGBRC(178))/CONSTS(4) ALGBRC(131) = (ALGBRC(120) - ALGBRC(74))/CONSTS(4) ALGBRC(189) = (((((ALGBRC(74) - ALGBRC(118))+ALGBRC(177)) - ALGBRC(88)) - ALGBRC(176)) - ALGBRC(178))/CONSTS(4) ALGBRC(132) = (- ALGBRC(57)+ALGBRC(118)+ALGBRC(87))/CONSTS(4) ALGBRC(161) = ((ALGBRC(87)+ALGBRC(157)) - ALGBRC(159))/CONSTS(4) ALGBRC(162) = (- ALGBRC(87)+ALGBRC(159))/CONSTS(4) ALGBRC(182) = CONSTS(51)*( ALGBRC(163)*STATES(2) - STATES(4)*ALGBRC(133)) ALGBRC(199) = (((((((((- ALGBRC(115)+ 2.00000*ALGBRC(57)) - ALGBRC(117))+ALGBRC(67)+ALGBRC(118)+ALGBRC(156)+ALGBRC(157)+ALGBRC(176)) - ALGBRC(159)) - 5.00000*ALGBRC(140)) - 2.00000*ALGBRC(123)) - 4.00000*ALGBRC(142))+ (CONSTS(8) - 1.00000)*ALGBRC(144)+ 2.00000*ALGBRC(180)+ALGBRC(145)) - ALGBRC(182))/CONSTS(4) ALGBRC(201) = - (( STATES(2)*CONSTS(293))/( CONSTS(135)*CONSTS(283))+( STATES(2)*ALGBRC(184))/( CONSTS(297)*ALGBRC(1))+( STATES(2)*ALGBRC(185))/( CONSTS(333)*ALGBRC(22))+( STATES(2)*CONSTS(263))/( CONSTS(345)*ALGBRC(27))+( STATES(2)*ALGBRC(124))/( CONSTS(297)*ALGBRC(30))+( STATES(2)*ALGBRC(125))/( CONSTS(333)*ALGBRC(32))+( STATES(2)*ALGBRC(186))/( CONSTS(352)*ALGBRC(34))+( STATES(2)*ALGBRC(146))/( CONSTS(111)*CONSTS(278))+( STATES(2)*ALGBRC(147))/( CONSTS(108)*CONSTS(277))+( STATES(2)*ALGBRC(148))/( CONSTS(114)*CONSTS(305))+( STATES(2)*ALGBRC(149))/( CONSTS(117)*CONSTS(321))+( STATES(2)*ALGBRC(160))/( CONSTS(99)*ALGBRC(47))+( STATES(2)*ALGBRC(132))/( CONSTS(90)*ALGBRC(37))+( STATES(2)*ALGBRC(126))/( CONSTS(120)*CONSTS(279))+( STATES(2)*ALGBRC(187))/( CONSTS(361)*ALGBRC(38))+( STATES(2)*ALGBRC(127))/( CONSTS(364)*ALGBRC(41))+( STATES(2)*ALGBRC(128))/( CONSTS(123)*CONSTS(307))+( STATES(2)*ALGBRC(129))/( CONSTS(366)*ALGBRC(42))+( STATES(2)*ALGBRC(130))/( CONSTS(358)*ALGBRC(36))+( STATES(2)*ALGBRC(188))/( CONSTS(367)*ALGBRC(43))+( STATES(2)*ALGBRC(131))/( CONSTS(370)*ALGBRC(44))+( STATES(2)*ALGBRC(189))/( CONSTS(371)*ALGBRC(45))+( STATES(2)*ALGBRC(161))/( CONSTS(377)*ALGBRC(49))+( STATES(2)*ALGBRC(162))/( CONSTS(379)*ALGBRC(50))+( STATES(2)*CONSTS(295))/( CONSTS(141)*CONSTS(285))+( STATES(2)*CONSTS(294))/( CONSTS(138)*CONSTS(284))+( STATES(2)*CONSTS(264))/( CONSTS(126)*CONSTS(280))+( STATES(2)*CONSTS(271))/( CONSTS(129)*CONSTS(281))+( STATES(2)*CONSTS(270))/( CONSTS(132)*CONSTS(282))+( STATES(2)*CONSTS(265))/( CONSTS(374)*ALGBRC(46)))+ALGBRC(199) ALGBRC(202) = - (( STATES(3)*CONSTS(293))/( CONSTS(136)*CONSTS(283))+( STATES(3)*ALGBRC(184))/( CONSTS(312)*ALGBRC(1))+( STATES(3)*ALGBRC(185))/( CONSTS(337)*ALGBRC(22))+( STATES(3)*CONSTS(263))/( CONSTS(348)*ALGBRC(27))+( STATES(3)*ALGBRC(124))/( CONSTS(312)*ALGBRC(30))+( STATES(3)*ALGBRC(125))/( CONSTS(337)*ALGBRC(32))+( STATES(3)*ALGBRC(186))/( CONSTS(354)*ALGBRC(34))+( STATES(3)*ALGBRC(146))/( CONSTS(112)*CONSTS(278))+( STATES(3)*ALGBRC(147))/( CONSTS(109)*CONSTS(277))+( STATES(3)*ALGBRC(148))/( CONSTS(115)*CONSTS(305))+( STATES(3)*ALGBRC(149))/( CONSTS(118)*CONSTS(321))+( STATES(3)*ALGBRC(160))/( CONSTS(375)*ALGBRC(47))+( STATES(3)*ALGBRC(132))/( CONSTS(360)*ALGBRC(37))+( STATES(3)*ALGBRC(126))/( CONSTS(121)*CONSTS(279))+( STATES(3)*ALGBRC(187))/( CONSTS(362)*ALGBRC(38))+( STATES(3)*ALGBRC(127))/( CONSTS(365)*ALGBRC(41))+( STATES(3)*ALGBRC(128))/( CONSTS(124)*CONSTS(307))+( STATES(3)*ALGBRC(129))/( CONSTS(93)*ALGBRC(42))+( STATES(3)*ALGBRC(130))/( CONSTS(88)*ALGBRC(36))+( STATES(3)*ALGBRC(188))/( CONSTS(368)*ALGBRC(43))+( STATES(3)*ALGBRC(131))/( CONSTS(95)*ALGBRC(44))+( STATES(3)*ALGBRC(189))/( CONSTS(372)*ALGBRC(45))+( STATES(3)*ALGBRC(161))/( CONSTS(378)*ALGBRC(49))+( STATES(3)*ALGBRC(162))/( CONSTS(380)*ALGBRC(50))+( STATES(3)*CONSTS(295))/( CONSTS(142)*CONSTS(285))+( STATES(3)*CONSTS(294))/( CONSTS(139)*CONSTS(284))+( STATES(3)*CONSTS(264))/( CONSTS(127)*CONSTS(280))+( STATES(3)*CONSTS(271))/( CONSTS(130)*CONSTS(281))+( STATES(3)*CONSTS(270))/( CONSTS(133)*CONSTS(282))+( STATES(3)*CONSTS(265))/( CONSTS(97)*ALGBRC(46)))+CONSTS(296) ALGBRC(200) = ALGBRC(182)/CONSTS(4) ALGBRC(203) = - (( STATES(4)*CONSTS(293))/( CONSTS(137)*CONSTS(283))+( STATES(4)*ALGBRC(184))/( CONSTS(326)*ALGBRC(1))+( STATES(4)*ALGBRC(185))/( CONSTS(341)*ALGBRC(22))+( STATES(4)*CONSTS(263))/( CONSTS(350)*ALGBRC(27))+( STATES(4)*ALGBRC(124))/( CONSTS(326)*ALGBRC(30))+( STATES(4)*ALGBRC(125))/( CONSTS(341)*ALGBRC(32))+( STATES(4)*ALGBRC(186))/( CONSTS(356)*ALGBRC(34))+( STATES(4)*ALGBRC(146))/( CONSTS(113)*CONSTS(278))+( STATES(4)*ALGBRC(147))/( CONSTS(110)*CONSTS(277))+( STATES(4)*ALGBRC(148))/( CONSTS(116)*CONSTS(305))+( STATES(4)*ALGBRC(149))/( CONSTS(119)*CONSTS(321))+( STATES(4)*ALGBRC(160))/( CONSTS(100)*ALGBRC(47))+( STATES(4)*ALGBRC(132))/( CONSTS(91)*ALGBRC(37))+( STATES(4)*ALGBRC(126))/( CONSTS(122)*CONSTS(279))+( STATES(4)*ALGBRC(187))/( CONSTS(363)*ALGBRC(38))+( STATES(4)*ALGBRC(127))/( CONSTS(92)*ALGBRC(41))+( STATES(4)*ALGBRC(128))/( CONSTS(125)*CONSTS(307))+( STATES(4)*ALGBRC(129))/( CONSTS(94)*ALGBRC(42))+( STATES(4)*ALGBRC(130))/( CONSTS(89)*ALGBRC(36))+( STATES(4)*ALGBRC(188))/( CONSTS(369)*ALGBRC(43))+( STATES(4)*ALGBRC(131))/( CONSTS(96)*ALGBRC(44))+( STATES(4)*ALGBRC(189))/( CONSTS(373)*ALGBRC(45))+( STATES(4)*ALGBRC(161))/( CONSTS(103)*ALGBRC(49))+( STATES(4)*ALGBRC(162))/( CONSTS(104)*ALGBRC(50))+( STATES(4)*CONSTS(295))/( CONSTS(143)*CONSTS(285))+( STATES(4)*CONSTS(294))/( CONSTS(140)*CONSTS(284))+( STATES(4)*CONSTS(264))/( CONSTS(128)*CONSTS(280))+( STATES(4)*CONSTS(271))/( CONSTS(131)*CONSTS(281))+( STATES(4)*CONSTS(270))/( CONSTS(134)*CONSTS(282))+( STATES(4)*CONSTS(265))/( CONSTS(98)*ALGBRC(46)))+ALGBRC(200) RATES(2) = ( ( ALGBRC(197)*ALGBRC(194) - ALGBRC(205)*ALGBRC(206))*ALGBRC(201)+ ( ALGBRC(206)*ALGBRC(190) - ALGBRC(191)*ALGBRC(197))*ALGBRC(202)+ ( ALGBRC(205)*ALGBRC(191) - ALGBRC(190)*ALGBRC(194))*ALGBRC(203))/ALGBRC(207) RATES(3) = ( ( ALGBRC(206)*ALGBRC(193) - ALGBRC(196)*ALGBRC(194))*ALGBRC(201)+ ( ALGBRC(196)*ALGBRC(191) - ALGBRC(204)*ALGBRC(206))*ALGBRC(202)+ ( ALGBRC(204)*ALGBRC(194) - ALGBRC(191)*ALGBRC(193))*ALGBRC(203))/ALGBRC(207) RATES(4) = ( ( ALGBRC(205)*ALGBRC(196) - ALGBRC(197)*ALGBRC(193))*ALGBRC(201)+ ( ALGBRC(204)*ALGBRC(197) - ALGBRC(196)*ALGBRC(190))*ALGBRC(202)+ ( ALGBRC(193)*ALGBRC(190) - ALGBRC(204)*ALGBRC(205))*ALGBRC(203))/ALGBRC(207) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(17) = CONSTS(6)*CONSTS(11)*(STATES(49) - STATES(50)) ALGBRC(16) = CONSTS(6)*CONSTS(247)*(STATES(47) - STATES(48)) ALGBRC(1) = 1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326) ALGBRC(30) = ALGBRC(1) ALGBRC(22) = 1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341) ALGBRC(32) = ALGBRC(22) ALGBRC(34) = 1.00000+STATES(2)/CONSTS(352)+STATES(3)/CONSTS(354)+STATES(4)/CONSTS(356) ALGBRC(36) = 1.00000+STATES(2)/CONSTS(358) ALGBRC(42) = 1.00000+STATES(2)/CONSTS(366) ALGBRC(43) = 1.00000+STATES(2)/CONSTS(367)+STATES(3)/CONSTS(368)+STATES(4)/CONSTS(369) ALGBRC(65) = ( CONSTS(320)*ALGBRC(36)*ALGBRC(43)*ALGBRC(30))/( STATES(2)*ALGBRC(32)*ALGBRC(42)*ALGBRC(34)) ALGBRC(66) = ( CONSTS(29)*CONSTS(187)*CONSTS(182)*CONSTS(183))/( ALGBRC(65)*CONSTS(178)*CONSTS(179)*CONSTS(186)) ALGBRC(67) = ( CONSTS(29)*ALGBRC(66)*( STATES(13)*STATES(24)*STATES(14) - ( STATES(19)*STATES(27)*STATES(12))/ALGBRC(65)))/( ALGBRC(66)*CONSTS(178)*CONSTS(179)*CONSTS(186)+ ALGBRC(66)*CONSTS(179)*CONSTS(186)*STATES(13)+ ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)+ ALGBRC(66)*CONSTS(186)*STATES(13)*STATES(24)+ ALGBRC(66)*CONSTS(185)*STATES(13)*STATES(14)+ ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)+ ALGBRC(66)*STATES(13)*STATES(24)*STATES(14)+( CONSTS(29)*CONSTS(183)*CONSTS(188)*STATES(19))/ALGBRC(65)+( CONSTS(29)*CONSTS(182)*CONSTS(187)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(189)*STATES(19)*STATES(27))/ALGBRC(65)+( CONSTS(29)*CONSTS(188)*STATES(19)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(187)*STATES(27)*STATES(12))/ALGBRC(65)+( CONSTS(29)*STATES(19)*STATES(27)*STATES(12))/ALGBRC(65)+( CONSTS(29)*CONSTS(188)*CONSTS(183)*STATES(13)*STATES(19))/( CONSTS(178)*ALGBRC(65))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(12))/CONSTS(183)+( CONSTS(29)*CONSTS(188)*CONSTS(183)*STATES(13)*STATES(24)*STATES(19))/( CONSTS(178)*CONSTS(179)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(12))/CONSTS(183)+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(19)*STATES(27))/( CONSTS(178)*ALGBRC(65))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(27)*STATES(12))/( CONSTS(182)*CONSTS(183))+( CONSTS(29)*CONSTS(183)*CONSTS(188)*STATES(13)*STATES(24)*STATES(14)*STATES(19))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( CONSTS(29)*CONSTS(181)*CONSTS(189)*STATES(13)*STATES(24)*STATES(14)*STATES(27))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(24)*STATES(19)*STATES(27))/( CONSTS(178)*CONSTS(179)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(27)*STATES(12))/( CONSTS(182)*CONSTS(183))+( ALGBRC(66)*CONSTS(184)*CONSTS(180)*STATES(24)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))+( ALGBRC(66)*CONSTS(178)*CONSTS(185)*STATES(14)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))+( CONSTS(29)*CONSTS(189)*STATES(13)*STATES(24)*STATES(14)*STATES(19)*STATES(27))/( CONSTS(178)*CONSTS(179)*CONSTS(180)*ALGBRC(65))+( ALGBRC(66)*CONSTS(184)*STATES(24)*STATES(14)*STATES(19)*STATES(27)*STATES(12))/( CONSTS(181)*CONSTS(182)*CONSTS(183))) ALGBRC(78) = ( CONSTS(313)*ALGBRC(32)*ALGBRC(1))/( ALGBRC(30)*ALGBRC(22)) ALGBRC(79) = ( CONSTS(33)*CONSTS(224)*CONSTS(219))/( ALGBRC(78)*CONSTS(217)*CONSTS(222)) ALGBRC(27) = 1.00000+STATES(2)/CONSTS(345)+STATES(3)/CONSTS(348)+STATES(4)/CONSTS(350) ALGBRC(29) = ( STATES(11)*(1.00000+STATES(2)/CONSTS(345)))/ALGBRC(27) ALGBRC(80) = 1.00000+ALGBRC(29)/CONSTS(225) ALGBRC(81) = TERNRY(STATES(12).GT.CONSTS(226).AND.STATES(9).GT.CONSTS(226), ( (( CONSTS(33)*ALGBRC(79))/ALGBRC(80))*STATES(12)*STATES(9))/( ALGBRC(79)*CONSTS(222)*STATES(12)+ ALGBRC(79)*CONSTS(221)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(13))/ALGBRC(78)+( CONSTS(33)*CONSTS(223)*STATES(1))/ALGBRC(78)+ ALGBRC(79)*STATES(12)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(12)*STATES(13))/( ALGBRC(78)*CONSTS(217))+( CONSTS(33)*STATES(13)*STATES(1))/ALGBRC(78)+( ALGBRC(79)*CONSTS(221)*STATES(9)*STATES(1))/CONSTS(220)), 0.00000) ALGBRC(82) = TERNRY(STATES(13).GT.CONSTS(226).AND.STATES(1).GT.CONSTS(226), ALGBRC(81) - (( (( CONSTS(33)*ALGBRC(79))/ALGBRC(80))*STATES(13)*STATES(1))/ALGBRC(78))/( ALGBRC(79)*CONSTS(222)*STATES(12)+ ALGBRC(79)*CONSTS(221)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(13))/ALGBRC(78)+( CONSTS(33)*CONSTS(223)*STATES(1))/ALGBRC(78)+ ALGBRC(79)*STATES(12)*STATES(9)+( CONSTS(33)*CONSTS(224)*STATES(12)*STATES(13))/( ALGBRC(78)*CONSTS(217))+( CONSTS(33)*STATES(13)*STATES(1))/ALGBRC(78)+( ALGBRC(79)*CONSTS(221)*STATES(9)*STATES(1))/CONSTS(220)), ALGBRC(81)) ALGBRC(88) = ( CONSTS(38)*( STATES(36)*STATES(30) - STATES(26)*STATES(35)))/( CONSTS(237)*CONSTS(240)*(2.00000+STATES(35)/CONSTS(239)+STATES(30)/CONSTS(240)+STATES(36)/CONSTS(237)+STATES(26)/CONSTS(238)+( STATES(35)*STATES(26))/( CONSTS(239)*CONSTS(238))+( STATES(30)*STATES(36))/( CONSTS(240)*CONSTS(237)))) ALGBRC(90) = CONSTS(6)*CONSTS(11)*(STATES(43) - STATES(36)) ALGBRC(89) = CONSTS(6)*CONSTS(11)*(STATES(40) - STATES(33)) ALGBRC(91) = CONSTS(6)*CONSTS(11)*(STATES(41) - STATES(17)) ALGBRC(92) = CONSTS(6)*CONSTS(11)*(STATES(42) - STATES(35)) ALGBRC(93) = CONSTS(6)*CONSTS(11)*(STATES(44) - STATES(39)) ALGBRC(94) = CONSTS(6)*CONSTS(11)*(STATES(45) - STATES(34)) ALGBRC(95) = CONSTS(6)*CONSTS(11)*(STATES(46) - STATES(37)) ALGBRC(99) = CONSTS(6)*CONSTS(9)*(STATES(55) - STATES(15)) ALGBRC(18) = 1.00000+STATES(6)/CONSTS(297)+STATES(7)/CONSTS(312)+STATES(8)/CONSTS(326) ALGBRC(21) = (( STATES(5)*STATES(7))/CONSTS(312))/ALGBRC(18) ALGBRC(23) = 1.00000+STATES(6)/CONSTS(333)+STATES(7)/CONSTS(337)+STATES(8)/CONSTS(341) ALGBRC(26) = (( STATES(10)*STATES(7))/CONSTS(337))/ALGBRC(23) ALGBRC(48) = 1.00000+STATES(6)/CONSTS(376) ALGBRC(100) = ( CONSTS(325)*ALGBRC(48)*ALGBRC(23))/( STATES(6)*ALGBRC(18)*CONSTS(276)) ALGBRC(102) = ( ALGBRC(100)*CONSTS(312)*ALGBRC(18))/( CONSTS(337)*ALGBRC(23)) ALGBRC(104) = ( (CONSTS(274)/( CONSTS(253)*CONSTS(256)))*( ALGBRC(21)*STATES(58) - ( STATES(59)*ALGBRC(26))/ALGBRC(102)))/(1.00000+ALGBRC(21)/CONSTS(252)+STATES(58)/CONSTS(253)+( ALGBRC(21)*STATES(58))/( CONSTS(253)*CONSTS(256))+STATES(59)/CONSTS(254)+ALGBRC(26)/CONSTS(255)+( STATES(59)*ALGBRC(26))/( CONSTS(255)*CONSTS(258))+( STATES(59)*STATES(58))/( CONSTS(260)*CONSTS(253))) ALGBRC(106) = CONSTS(19) - STATES(60) ALGBRC(108) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(261)*( (( CONSTS(310)*STATES(10))/ALGBRC(23))*STATES(60)*STATES(6) - (STATES(5)/ALGBRC(18))*ALGBRC(106))) ALGBRC(97) = CONSTS(6)*CONSTS(10)*(STATES(5) - STATES(52)) ALGBRC(28) = 1.00000+STATES(6)/CONSTS(345)+STATES(7)/CONSTS(348)+STATES(8)/CONSTS(350) ALGBRC(110) = ( CONSTS(382)*ALGBRC(18)*ALGBRC(28))/( ALGBRC(23)*ALGBRC(23)) ALGBRC(112) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(262)*( ALGBRC(110)*STATES(10)*STATES(10) - STATES(53)*STATES(5))) ALGBRC(96) = CONSTS(6)*CONSTS(10)*(STATES(10) - STATES(51)) ALGBRC(98) = CONSTS(6)*CONSTS(10)*(STATES(53) - STATES(54)) ALGBRC(114) = CONSTS(15) - STATES(22) ALGBRC(46) = 1.00000+STATES(2)/CONSTS(374) ALGBRC(47) = 1.00000+STATES(3)/CONSTS(375) ALGBRC(51) = ( CONSTS(315)*ALGBRC(46)*CONSTS(279)*CONSTS(278))/( STATES(2)*ALGBRC(47)*ALGBRC(36)*CONSTS(277)) ALGBRC(52) = 1.00000+STATES(21)/CONSTS(147) ALGBRC(53) = 1.00000+STATES(22)/CONSTS(148) ALGBRC(115) = TERNRY(STATES(18).GT.CONSTS(149).AND.STATES(19).GT.CONSTS(149).AND.ALGBRC(114).GT.CONSTS(149), ( CONSTS(23)*( STATES(18)*STATES(19)*ALGBRC(114) - ( STATES(20)*STATES(21)*STATES(22))/ALGBRC(51)))/( CONSTS(146)*ALGBRC(53)*STATES(18)*STATES(19)+ CONSTS(145)*ALGBRC(52)*STATES(18)*ALGBRC(114)+ CONSTS(144)*STATES(19)*ALGBRC(114)+ STATES(18)*STATES(19)*ALGBRC(114)), 0.00000) ALGBRC(37) = 1.00000+STATES(3)/CONSTS(360) ALGBRC(38) = 1.00000+STATES(2)/CONSTS(361)+STATES(3)/CONSTS(362)+STATES(4)/CONSTS(363) ALGBRC(54) = ( CONSTS(316)*ALGBRC(38)*ALGBRC(36))/( STATES(2) ** 2.00000*CONSTS(279)*ALGBRC(37)) ALGBRC(39) = ( STATES(16)*(1.00000+STATES(2)/CONSTS(361)))/ALGBRC(38) ALGBRC(55) = 1.00000+ALGBRC(39)/CONSTS(153) ALGBRC(19) = ( STATES(1)*(1.00000+STATES(2)/CONSTS(297)))/ALGBRC(1) ALGBRC(24) = ( STATES(9)*(1.00000+STATES(2)/CONSTS(333)))/ALGBRC(22) ALGBRC(56) = 1.00000+ALGBRC(19)/CONSTS(154)+ALGBRC(24)/CONSTS(155)+ALGBRC(29)/CONSTS(156)+STATES(19)/CONSTS(157)+STATES(24)/CONSTS(158) ALGBRC(57) = ( CONSTS(24)*( STATES(21)*STATES(23) - ( STATES(19)*STATES(16))/ALGBRC(54)))/( CONSTS(152)*CONSTS(151)*ALGBRC(55)+ CONSTS(150)*ALGBRC(55)*STATES(21)+ CONSTS(151)*ALGBRC(56)*STATES(23)+ STATES(23)*STATES(21)) ALGBRC(41) = 1.00000+STATES(2)/CONSTS(364)+STATES(3)/CONSTS(365) ALGBRC(58) = ( CONSTS(317)*ALGBRC(41))/ALGBRC(38) ALGBRC(59) = ( CONSTS(25)*CONSTS(160))/( CONSTS(159)*ALGBRC(58)) ALGBRC(60) = ( CONSTS(25)*ALGBRC(59)*(STATES(16) - STATES(25)/ALGBRC(58)))/( CONSTS(159)*ALGBRC(59)+ ALGBRC(59)*STATES(16)+( CONSTS(25)*STATES(25))/ALGBRC(58)) ALGBRC(61) = ( CONSTS(318)*ALGBRC(46)*CONSTS(307)*CONSTS(278))/( STATES(2) ** 2.00000*CONSTS(277)*ALGBRC(41)) ALGBRC(62) = 1.00000+ (CONSTS(167)/ALGBRC(24))*(1.00000+ALGBRC(19)/CONSTS(166)) ALGBRC(116) = TERNRY(ALGBRC(114).GT.CONSTS(168).AND.STATES(25).GT.CONSTS(168), ( CONSTS(26)*(1.00000 - ( STATES(20)*STATES(26)*STATES(22))/( ALGBRC(61)*ALGBRC(114)*STATES(25))))/(1.00000+ CONSTS(162)/STATES(25) ** CONSTS(165)*ALGBRC(62)+ (CONSTS(161)/ALGBRC(114))*(1.00000+ CONSTS(163)/STATES(25) ** CONSTS(165)*ALGBRC(62)+( STATES(22)*ALGBRC(62))/CONSTS(164))), 0.00000) ALGBRC(63) = ( CONSTS(319)*ALGBRC(46)*ALGBRC(42)*CONSTS(278))/( STATES(2)*ALGBRC(36)*CONSTS(307)*CONSTS(277)) ALGBRC(64) = 1.00000+ (CONSTS(176)/ALGBRC(24))*(1.00000+ALGBRC(19)/CONSTS(175)) ALGBRC(117) = TERNRY(STATES(26).GT.CONSTS(177).AND.STATES(19).GT.CONSTS(177).AND.ALGBRC(114).GT.CONSTS(177), (( CONSTS(27)*(1.00000 - ( STATES(20)*STATES(24)*STATES(22))/( ALGBRC(63)*STATES(26)*STATES(19)*ALGBRC(114))))/(1.00000+ (CONSTS(169)/STATES(26))*ALGBRC(64)+ (CONSTS(170)/STATES(19))*(1.00000+STATES(24)/CONSTS(172))+ (CONSTS(171)/ALGBRC(114))*(1.00000+STATES(22)/CONSTS(173))))/(1.00000+STATES(22)/CONSTS(174)), 0.00000) ALGBRC(49) = 1.00000+STATES(2)/CONSTS(377)+STATES(3)/CONSTS(378) ALGBRC(50) = 1.00000+STATES(2)/CONSTS(379)+STATES(3)/CONSTS(380) ALGBRC(83) = ( CONSTS(329)*ALGBRC(49)*ALGBRC(37))/( ALGBRC(50)*CONSTS(307)) ALGBRC(84) = ( CONSTS(34)*CONSTS(234)*CONSTS(229))/( ALGBRC(83)*CONSTS(227)*CONSTS(232)) ALGBRC(85) = 1.00000+STATES(26)/CONSTS(235) ALGBRC(86) = TERNRY(STATES(31).GT.CONSTS(236).AND.STATES(26).GT.CONSTS(236), ( CONSTS(34)*ALGBRC(84)*STATES(31)*STATES(26))/( ALGBRC(84)*CONSTS(232)*STATES(31)+ ALGBRC(84)*CONSTS(231)*ALGBRC(85)*STATES(26)+( CONSTS(34)*CONSTS(234)*ALGBRC(85)*STATES(23))/ALGBRC(83)+( CONSTS(34)*CONSTS(233)*STATES(32))/ALGBRC(83)+ ALGBRC(84)*STATES(31)*STATES(26)+( CONSTS(34)*CONSTS(234)*STATES(31)*STATES(23))/( ALGBRC(83)*CONSTS(227))+( CONSTS(34)*STATES(23)*STATES(32))/ALGBRC(83)+( ALGBRC(84)*CONSTS(231)*STATES(26)*STATES(32))/CONSTS(230)), 0.00000) ALGBRC(87) = TERNRY(STATES(23).GT.CONSTS(236).AND.STATES(32).GT.CONSTS(236), ALGBRC(86) - (( CONSTS(34)*ALGBRC(84)*STATES(23)*STATES(32))/ALGBRC(83))/( ALGBRC(84)*CONSTS(232)*STATES(31)+ ALGBRC(84)*CONSTS(231)*ALGBRC(85)*STATES(26)+( CONSTS(34)*CONSTS(234)*ALGBRC(85)*STATES(23))/ALGBRC(83)+( CONSTS(34)*CONSTS(233)*STATES(32))/ALGBRC(83)+ ALGBRC(84)*STATES(31)*STATES(26)+( CONSTS(34)*CONSTS(234)*STATES(31)*STATES(23))/( ALGBRC(83)*CONSTS(227))+( CONSTS(34)*STATES(23)*STATES(32))/ALGBRC(83)+( ALGBRC(84)*CONSTS(231)*STATES(26)*STATES(32))/CONSTS(230)), ALGBRC(86)) ALGBRC(45) = 1.00000+STATES(2)/CONSTS(371)+STATES(3)/CONSTS(372)+STATES(4)/CONSTS(373) ALGBRC(75) = ( CONSTS(322)*ALGBRC(37)*CONSTS(278))/( STATES(2)*ALGBRC(45)*CONSTS(277)) ALGBRC(76) = ( CONSTS(32)*CONSTS(209)*CONSTS(212))/( ALGBRC(75)*CONSTS(206)*CONSTS(211)) ALGBRC(77) = 1.00000+ALGBRC(19)/CONSTS(214)+ALGBRC(24)/CONSTS(215)+ALGBRC(29)/CONSTS(216) ALGBRC(118) = ( CONSTS(32)*ALGBRC(76)*( ALGBRC(114)*STATES(30) - ( STATES(23)*STATES(22))/ALGBRC(75)))/( ALGBRC(76)*CONSTS(206)*CONSTS(211)*ALGBRC(77)+ ALGBRC(76)*CONSTS(211)*ALGBRC(114)+ ALGBRC(76)*CONSTS(210)*ALGBRC(77)*STATES(30)+( CONSTS(32)*CONSTS(213)*ALGBRC(77)*STATES(23))/ALGBRC(75)+( CONSTS(32)*CONSTS(212)*STATES(22))/ALGBRC(75)+ ALGBRC(76)*ALGBRC(114)*STATES(30)+( CONSTS(32)*CONSTS(213)*ALGBRC(114)*STATES(23))/( ALGBRC(75)*CONSTS(206))+( CONSTS(32)*STATES(23)*STATES(22))/ALGBRC(75)+( ALGBRC(76)*CONSTS(210)*STATES(30)*STATES(22))/CONSTS(209)+( ALGBRC(76)*ALGBRC(114)*STATES(30)*STATES(23))/CONSTS(208)+( CONSTS(32)*STATES(30)*STATES(23)*STATES(22))/( CONSTS(207)*ALGBRC(75))) ALGBRC(119) = CONSTS(16) - STATES(28) ALGBRC(44) = 1.00000+STATES(2)/CONSTS(370) ALGBRC(68) = ( CONSTS(335)*ALGBRC(44)*CONSTS(305))/( ALGBRC(43)*CONSTS(321)) ALGBRC(69) = ( CONSTS(30)*CONSTS(194)*CONSTS(191))/( ALGBRC(68)*CONSTS(190)*CONSTS(193)) ALGBRC(70) = (1.00000+STATES(23)/CONSTS(196)+STATES(27)/CONSTS(197)+STATES(29)/CONSTS(198))/(1.00000+STATES(27)/CONSTS(197)+STATES(29)/CONSTS(198)) ALGBRC(120) = ( CONSTS(30)*ALGBRC(69)*( STATES(27)*ALGBRC(119) - ( STATES(28)*STATES(29))/ALGBRC(68)))/( ALGBRC(69)*CONSTS(190)*CONSTS(193)*ALGBRC(70)+ ALGBRC(69)*CONSTS(193)*STATES(27)+ ALGBRC(69)*CONSTS(192)*ALGBRC(70)*ALGBRC(119)+( CONSTS(30)*CONSTS(195)*ALGBRC(70)*STATES(28))/ALGBRC(68)+( CONSTS(30)*CONSTS(194)*STATES(29))/ALGBRC(68)+ ALGBRC(69)*STATES(27)*ALGBRC(119)+( CONSTS(30)*CONSTS(195)*STATES(27)*STATES(28))/( ALGBRC(68)*CONSTS(190))+( ALGBRC(69)*CONSTS(192)*ALGBRC(119)*STATES(29))/CONSTS(191)+( CONSTS(30)*STATES(28)*STATES(29))/ALGBRC(68)) ALGBRC(71) = ( CONSTS(343)*ALGBRC(45))/ALGBRC(44) ALGBRC(72) = ( CONSTS(31)*CONSTS(200))/( ALGBRC(71)*CONSTS(199)) ALGBRC(31) = ( STATES(12)*(1.00000+STATES(2)/CONSTS(314)))/ALGBRC(30) ALGBRC(33) = ( STATES(13)*(1.00000+STATES(2)/CONSTS(338)))/ALGBRC(32) ALGBRC(73) = 1.00000+STATES(16)/CONSTS(201)+ALGBRC(19)/CONSTS(202)+ALGBRC(24)/CONSTS(203)+ALGBRC(31)/CONSTS(204)+ALGBRC(33)/CONSTS(205) ALGBRC(74) = ( ALGBRC(72)*CONSTS(31)*(STATES(29) - STATES(30)/ALGBRC(71)))/( CONSTS(199)*ALGBRC(72)*ALGBRC(73)+ ALGBRC(72)*STATES(29)+( CONSTS(31)*STATES(30))/ALGBRC(71)) ALGBRC(133) = STATES(6) ALGBRC(101) = 1.00000*EXP(- (CONSTS(288)+ 4.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(139) = ( ALGBRC(101)*STATES(2) ** 5.00000)/ALGBRC(133) ** 4.00000 ALGBRC(140) = CONSTS(42)*( ALGBRC(139)*STATES(22)*ALGBRC(119) - ALGBRC(114)*STATES(28)) ALGBRC(103) = 1.00000*EXP(- (CONSTS(289)+ 2.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(105) = ( ALGBRC(103)*STATES(2) ** 2.00000)/STATES(2) ** 4.00000 ALGBRC(107) = TERNRY(STATES(28).GT.CONSTS(248), STATES(28), CONSTS(248)) ALGBRC(121) = TERNRY(ALGBRC(119).GT.CONSTS(248), ALGBRC(119), CONSTS(248)) ALGBRC(122) = CONSTS(17) - STATES(56) ALGBRC(123) = (( CONSTS(43)*(1.00000+STATES(14)/CONSTS(44)))/(1.00000+STATES(14)/CONSTS(45)))*( ALGBRC(105) ** (1.0 / 2)*ALGBRC(122)* ALGBRC(107) ** (1.0 / 2) - STATES(56)* ALGBRC(121) ** (1.0 / 2)) ALGBRC(109) = EXP(- (CONSTS(290)+ 4.00000*CONSTS(3)*STATES(38))/CONSTS(275)) ALGBRC(141) = ( ALGBRC(109)*STATES(2) ** 4.00000)/( ALGBRC(133) ** 2.00000*1.00000 ** 2.00000) ALGBRC(111) = TERNRY(STATES(57).GT.CONSTS(249), STATES(57), CONSTS(249)) ALGBRC(142) = CONSTS(46)*( (ALGBRC(111)/(ALGBRC(111)+CONSTS(14)))*EXP(( CONSTS(3)*STATES(38))/CONSTS(275))*(STATES(56)/CONSTS(17)))*( ALGBRC(141) ** 1.00000/2.00000*STATES(56)*ALGBRC(111)/1.00000 ** 1.00000/4.00000 - ALGBRC(122)) ALGBRC(150) = STATES(7) ALGBRC(153) = 1.00000+ALGBRC(150)/CONSTS(375) ALGBRC(156) = CONSTS(35)*( (STATES(33)/ALGBRC(153))*ALGBRC(133) - (STATES(18)/ALGBRC(47))*STATES(2)) ALGBRC(154) = 1.00000+ALGBRC(133)/CONSTS(377)+ALGBRC(150)/CONSTS(378) ALGBRC(157) = CONSTS(36)*( (STATES(34)/ALGBRC(154))*ALGBRC(133) - (STATES(32)/ALGBRC(49))*STATES(2)) ALGBRC(155) = 1.00000+ALGBRC(133)/CONSTS(379)+ALGBRC(150)/CONSTS(380) ALGBRC(158) = ( EXP(( - CONSTS(3)*STATES(38))/CONSTS(275))*ALGBRC(50)*ALGBRC(154))/( ALGBRC(155)*ALGBRC(49)) ALGBRC(159) = ( CONSTS(41)*( ALGBRC(158)*STATES(37)*STATES(32)*STATES(2) - STATES(31)*STATES(34)*ALGBRC(133)))/( ALGBRC(158)*CONSTS(241)*CONSTS(244)*CONSTS(286)*( 2.00000*CONSTS(245)+( CONSTS(245)*STATES(37))/CONSTS(241)+( STATES(37)*STATES(32)*STATES(2))/( CONSTS(241)*CONSTS(244)*CONSTS(286))+( CONSTS(245)*STATES(31)*ALGBRC(133))/( CONSTS(242)*CONSTS(286))+( STATES(31)*STATES(34)*ALGBRC(133))/( CONSTS(242)*CONSTS(243)*CONSTS(286))+( CONSTS(245)*STATES(31))/CONSTS(242)+( CONSTS(245)*STATES(37)*STATES(2))/( CONSTS(241)*CONSTS(286))+( CONSTS(245)*STATES(2))/CONSTS(286)+( CONSTS(245)*STATES(34)*ALGBRC(133))/( CONSTS(243)*CONSTS(286))+( CONSTS(245)*ALGBRC(133))/CONSTS(286)+( CONSTS(245)*STATES(32)*STATES(2))/( CONSTS(244)*CONSTS(286)))) ALGBRC(40) = (( STATES(16)*STATES(2))/CONSTS(361))/ALGBRC(38) ALGBRC(163) = STATES(8) ALGBRC(171) = 1.00000+ALGBRC(133)/CONSTS(361)+ALGBRC(150)/CONSTS(362)+ALGBRC(163)/CONSTS(363) ALGBRC(173) = (( STATES(17)*ALGBRC(133))/CONSTS(361))/ALGBRC(171) ALGBRC(175) = 1.00000+ALGBRC(133)/CONSTS(371)+ALGBRC(150)/CONSTS(372)+ALGBRC(163)/CONSTS(373) ALGBRC(176) = CONSTS(37)*(( ALGBRC(173)*STATES(30))/ALGBRC(45) - ( ALGBRC(40)*STATES(35))/ALGBRC(175)) ALGBRC(174) = 1.00000+ALGBRC(133)/CONSTS(367)+ALGBRC(150)/CONSTS(368)+ALGBRC(163)/CONSTS(369) ALGBRC(178) = CONSTS(39)*(( (STATES(39)/ALGBRC(174))*STATES(30))/ALGBRC(45) - ( (STATES(27)/ALGBRC(43))*STATES(35))/ALGBRC(175)) ALGBRC(169) = 1.00000+ALGBRC(133)/CONSTS(352)+ALGBRC(150)/CONSTS(354)+ALGBRC(163)/CONSTS(356) ALGBRC(177) = CONSTS(40)*(( (STATES(35)/ALGBRC(175))*STATES(14))/ALGBRC(34) - ( (STATES(30)/ALGBRC(45))*STATES(15))/ALGBRC(169)) ALGBRC(113) = EXP(- (CONSTS(291) - CONSTS(8)*CONSTS(3)*STATES(38))/CONSTS(275))/1.00000 ALGBRC(143) = ( (( ALGBRC(113)*ALGBRC(133) ** CONSTS(8))/STATES(2) ** CONSTS(8) - 1.00000)*ALGBRC(1))/( ALGBRC(22)*ALGBRC(34)) ALGBRC(144) = CONSTS(47)*( ALGBRC(143)*STATES(9)*STATES(14) - STATES(1)) ALGBRC(164) = 1.00000+ALGBRC(133)/CONSTS(297)+ALGBRC(150)/CONSTS(312)+ALGBRC(163)/CONSTS(326) ALGBRC(166) = 1.00000+ALGBRC(133)/CONSTS(333)+ALGBRC(150)/CONSTS(337)+ALGBRC(163)/CONSTS(341) ALGBRC(179) = TERNRY(STATES(51).GT.CONSTS(250).AND.STATES(1).GT.CONSTS(250), CONSTS(48)*((STATES(51)/ALGBRC(166))/(STATES(51)/ALGBRC(166)+ (STATES(52)/ALGBRC(164))*EXP(( CONSTS(13)*CONSTS(3)*STATES(38))/CONSTS(275))) - (STATES(9)/ALGBRC(22))/(STATES(9)/ALGBRC(22)+ (STATES(1)/ALGBRC(1))*EXP(( CONSTS(13)*(CONSTS(3) - 1.00000)*STATES(38))/CONSTS(275))))*((STATES(51)/ALGBRC(166))/(STATES(51)/ALGBRC(166)+CONSTS(12))), 0.00000) ALGBRC(145) = TERNRY(ABS(STATES(38)).GT.CONSTS(251), ( CONSTS(52)*STATES(38)*( ALGBRC(133)*EXP(( CONSTS(3)*STATES(38))/CONSTS(275)) - STATES(2)))/(EXP(( CONSTS(3)*STATES(38))/CONSTS(275)) - 1.00000), ( CONSTS(52)*CONSTS(275)*(ALGBRC(133) - STATES(2)))/CONSTS(3)) ALGBRC(35) = (( STATES(14)*STATES(2))/CONSTS(352))/ALGBRC(34) ALGBRC(170) = (( STATES(15)*ALGBRC(133))/CONSTS(352))/ALGBRC(169) ALGBRC(180) = ( CONSTS(49)*( ALGBRC(170)*ALGBRC(133) - ALGBRC(35)*STATES(2)))/( CONSTS(50)*(1.00000+ALGBRC(170)/CONSTS(50))*(1.00000+ALGBRC(35)/CONSTS(50))) ALGBRC(168) = 1.00000+ALGBRC(133)/CONSTS(345)+ALGBRC(150)/CONSTS(348)+ALGBRC(163)/CONSTS(350) ALGBRC(181) = ( CONSTS(382)*ALGBRC(164)*ALGBRC(168))/( ALGBRC(166)*ALGBRC(166)) ALGBRC(183) = TERNRY(CONSTS(1).EQ.1.00000.OR.CONSTS(1).EQ.2.00000.OR.CONSTS(1).EQ.5.00000, 0.00000, CONSTS(262)*( ALGBRC(181)*STATES(51)*STATES(51) - STATES(54)*STATES(52))) ALGBRC(190) = - (( STATES(2)*STATES(2))/( CONSTS(135)*CONSTS(136)*CONSTS(283) ** 2.00000)+( STATES(2)*STATES(1))/( CONSTS(297)*CONSTS(312)*ALGBRC(1) ** 2.00000)+( STATES(2)*STATES(9))/( CONSTS(333)*CONSTS(337)*ALGBRC(22) ** 2.00000)+( STATES(2)*STATES(11))/( CONSTS(345)*CONSTS(348)*ALGBRC(27) ** 2.00000)+( STATES(2)*STATES(12))/( CONSTS(297)*CONSTS(312)*ALGBRC(30) ** 2.00000)+( STATES(2)*STATES(13))/( CONSTS(333)*CONSTS(337)*ALGBRC(32) ** 2.00000)+( STATES(2)*STATES(14))/( CONSTS(352)*CONSTS(354)*ALGBRC(34) ** 2.00000)+( STATES(2)*STATES(22))/( CONSTS(111)*CONSTS(112)*CONSTS(278) ** 2.00000)+( STATES(2)*ALGBRC(114))/( CONSTS(108)*CONSTS(109)*CONSTS(277) ** 2.00000)+( STATES(2)*STATES(28))/( CONSTS(114)*CONSTS(115)*CONSTS(305) ** 2.00000)+( STATES(2)*ALGBRC(119))/( CONSTS(117)*CONSTS(118)*CONSTS(321) ** 2.00000)+( STATES(2)*STATES(18))/( CONSTS(99)*CONSTS(375)*ALGBRC(47) ** 2.00000)+( STATES(2)*STATES(23))/( CONSTS(90)*CONSTS(360)*ALGBRC(37) ** 2.00000)+( STATES(2)*STATES(21))/( CONSTS(120)*CONSTS(121)*CONSTS(279) ** 2.00000)+( STATES(2)*STATES(16))/( CONSTS(361)*CONSTS(362)*ALGBRC(38) ** 2.00000)+( STATES(2)*STATES(25))/( CONSTS(364)*CONSTS(365)*ALGBRC(41) ** 2.00000)+( STATES(2)*STATES(26))/( CONSTS(123)*CONSTS(124)*CONSTS(307) ** 2.00000)+( STATES(2)*STATES(24))/( CONSTS(366)*CONSTS(93)*ALGBRC(42) ** 2.00000)+( STATES(2)*STATES(19))/( CONSTS(358)*CONSTS(88)*ALGBRC(36) ** 2.00000)+( STATES(2)*STATES(27))/( CONSTS(367)*CONSTS(368)*ALGBRC(43) ** 2.00000)+( STATES(2)*STATES(29))/( CONSTS(370)*CONSTS(95)*ALGBRC(44) ** 2.00000)+( STATES(2)*STATES(30))/( CONSTS(371)*CONSTS(372)*ALGBRC(45) ** 2.00000)+( STATES(2)*STATES(32))/( CONSTS(377)*CONSTS(378)*ALGBRC(49) ** 2.00000)+( STATES(2)*STATES(31))/( CONSTS(379)*CONSTS(380)*ALGBRC(50) ** 2.00000)+( STATES(2)*STATES(4))/( CONSTS(141)*CONSTS(142)*CONSTS(285) ** 2.00000)+( STATES(2)*STATES(3))/( CONSTS(138)*CONSTS(139)*CONSTS(284) ** 2.00000)+( STATES(2)*STATES(57))/( CONSTS(126)*CONSTS(127)*CONSTS(280) ** 2.00000)+( STATES(2)*CONSTS(292))/( CONSTS(129)*CONSTS(130)*CONSTS(281) ** 2.00000)+( STATES(2)*CONSTS(311))/( CONSTS(132)*CONSTS(133)*CONSTS(282) ** 2.00000)+( STATES(2)*STATES(20))/( CONSTS(374)*CONSTS(97)*ALGBRC(46) ** 2.00000)) ALGBRC(191) = - (( STATES(2)*STATES(2))/( CONSTS(135)*CONSTS(137)*CONSTS(283) ** 2.00000)+( STATES(2)*STATES(1))/( CONSTS(297)*CONSTS(326)*ALGBRC(1) ** 2.00000)+( STATES(2)*STATES(9))/( CONSTS(333)*CONSTS(341)*ALGBRC(22) ** 2.00000)+( STATES(2)*STATES(11))/( CONSTS(345)*CONSTS(350)*ALGBRC(27) ** 2.00000)+( STATES(2)*STATES(12))/( CONSTS(297)*CONSTS(326)*ALGBRC(30) ** 2.00000)+( STATES(2)*STATES(13))/( CONSTS(333)*CONSTS(341)*ALGBRC(32) ** 2.00000)+( STATES(2)*STATES(14))/( CONSTS(352)*CONSTS(356)*ALGBRC(34) ** 2.00000)+( STATES(2)*STATES(22))/( CONSTS(111)*CONSTS(113)*CONSTS(278) ** 2.00000)+( STATES(2)*ALGBRC(114))/( CONSTS(108)*CONSTS(110)*CONSTS(277) ** 2.00000)+( STATES(2)*STATES(28))/( CONSTS(114)*CONSTS(116)*CONSTS(305) ** 2.00000)+( STATES(2)*ALGBRC(119))/( CONSTS(117)*CONSTS(119)*CONSTS(321) ** 2.00000)+( STATES(2)*STATES(18))/( CONSTS(99)*CONSTS(100)*ALGBRC(47) ** 2.00000)+( STATES(2)*STATES(23))/( CONSTS(90)*CONSTS(91)*ALGBRC(37) ** 2.00000)+( STATES(2)*STATES(21))/( CONSTS(120)*CONSTS(122)*CONSTS(279) ** 2.00000)+( STATES(2)*STATES(16))/( CONSTS(361)*CONSTS(363)*ALGBRC(38) ** 2.00000)+( STATES(2)*STATES(25))/( CONSTS(364)*CONSTS(92)*ALGBRC(41) ** 2.00000)+( STATES(2)*STATES(26))/( CONSTS(123)*CONSTS(125)*CONSTS(307) ** 2.00000)+( STATES(2)*STATES(24))/( CONSTS(366)*CONSTS(94)*ALGBRC(42) ** 2.00000)+( STATES(2)*STATES(19))/( CONSTS(358)*CONSTS(89)*ALGBRC(36) ** 2.00000)+( STATES(2)*STATES(27))/( CONSTS(367)*CONSTS(369)*ALGBRC(43) ** 2.00000)+( STATES(2)*STATES(29))/( CONSTS(370)*CONSTS(96)*ALGBRC(44) ** 2.00000)+( STATES(2)*STATES(30))/( CONSTS(371)*CONSTS(373)*ALGBRC(45) ** 2.00000)+( STATES(2)*STATES(32))/( CONSTS(377)*CONSTS(103)*ALGBRC(49) ** 2.00000)+( STATES(2)*STATES(31))/( CONSTS(379)*CONSTS(104)*ALGBRC(50) ** 2.00000)+( STATES(2)*STATES(4))/( CONSTS(141)*CONSTS(143)*CONSTS(285) ** 2.00000)+( STATES(2)*STATES(3))/( CONSTS(138)*CONSTS(140)*CONSTS(284) ** 2.00000)+( STATES(2)*STATES(57))/( CONSTS(126)*CONSTS(128)*CONSTS(280) ** 2.00000)+( STATES(2)*CONSTS(292))/( CONSTS(129)*CONSTS(131)*CONSTS(281) ** 2.00000)+( STATES(2)*CONSTS(311))/( CONSTS(132)*CONSTS(134)*CONSTS(282) ** 2.00000)+( STATES(2)*STATES(20))/( CONSTS(374)*CONSTS(98)*ALGBRC(46) ** 2.00000)) ALGBRC(194) = - (( STATES(3)*STATES(2))/( CONSTS(136)*CONSTS(137)*CONSTS(283) ** 2.00000)+( STATES(3)*STATES(1))/( CONSTS(312)*CONSTS(326)*ALGBRC(1) ** 2.00000)+( STATES(3)*STATES(9))/( CONSTS(337)*CONSTS(341)*ALGBRC(22) ** 2.00000)+( STATES(3)*STATES(11))/( CONSTS(348)*CONSTS(350)*ALGBRC(27) ** 2.00000)+( STATES(3)*STATES(12))/( CONSTS(312)*CONSTS(326)*ALGBRC(30) ** 2.00000)+( STATES(3)*STATES(13))/( CONSTS(337)*CONSTS(341)*ALGBRC(32) ** 2.00000)+( STATES(3)*STATES(14))/( CONSTS(354)*CONSTS(356)*ALGBRC(34) ** 2.00000)+( STATES(3)*STATES(22))/( CONSTS(112)*CONSTS(113)*CONSTS(278) ** 2.00000)+( STATES(3)*ALGBRC(114))/( CONSTS(109)*CONSTS(110)*CONSTS(277) ** 2.00000)+( STATES(3)*STATES(28))/( CONSTS(115)*CONSTS(116)*CONSTS(305) ** 2.00000)+( STATES(3)*ALGBRC(119))/( CONSTS(118)*CONSTS(119)*CONSTS(321) ** 2.00000)+( STATES(3)*STATES(18))/( CONSTS(375)*CONSTS(100)*ALGBRC(47) ** 2.00000)+( STATES(3)*STATES(23))/( CONSTS(360)*CONSTS(91)*ALGBRC(37) ** 2.00000)+( STATES(3)*STATES(21))/( CONSTS(121)*CONSTS(122)*CONSTS(279) ** 2.00000)+( STATES(3)*STATES(16))/( CONSTS(362)*CONSTS(363)*ALGBRC(38) ** 2.00000)+( STATES(3)*STATES(25))/( CONSTS(365)*CONSTS(92)*ALGBRC(41) ** 2.00000)+( STATES(3)*STATES(26))/( CONSTS(124)*CONSTS(125)*CONSTS(307) ** 2.00000)+( STATES(3)*STATES(24))/( CONSTS(93)*CONSTS(94)*ALGBRC(42) ** 2.00000)+( STATES(3)*STATES(19))/( CONSTS(88)*CONSTS(89)*ALGBRC(36) ** 2.00000)+( STATES(3)*STATES(27))/( CONSTS(368)*CONSTS(369)*ALGBRC(43) ** 2.00000)+( STATES(3)*STATES(29))/( CONSTS(95)*CONSTS(96)*ALGBRC(44) ** 2.00000)+( STATES(3)*STATES(30))/( CONSTS(372)*CONSTS(373)*ALGBRC(45) ** 2.00000)+( STATES(3)*STATES(32))/( CONSTS(378)*CONSTS(103)*ALGBRC(49) ** 2.00000)+( STATES(3)*STATES(31))/( CONSTS(380)*CONSTS(104)*ALGBRC(50) ** 2.00000)+( STATES(3)*STATES(4))/( CONSTS(142)*CONSTS(143)*CONSTS(285) ** 2.00000)+( STATES(3)*STATES(3))/( CONSTS(139)*CONSTS(140)*CONSTS(284) ** 2.00000)+( STATES(3)*STATES(57))/( CONSTS(127)*CONSTS(128)*CONSTS(280) ** 2.00000)+( STATES(3)*CONSTS(292))/( CONSTS(130)*CONSTS(131)*CONSTS(281) ** 2.00000)+( STATES(3)*CONSTS(311))/( CONSTS(133)*CONSTS(134)*CONSTS(282) ** 2.00000)+( STATES(3)*STATES(20))/( CONSTS(97)*CONSTS(98)*ALGBRC(46) ** 2.00000)) ALGBRC(197) = - (( STATES(4)*STATES(2))/( CONSTS(137)*CONSTS(136)*CONSTS(283) ** 2.00000)+( STATES(4)*STATES(1))/( CONSTS(326)*CONSTS(312)*ALGBRC(1) ** 2.00000)+( STATES(4)*STATES(9))/( CONSTS(341)*CONSTS(337)*ALGBRC(22) ** 2.00000)+( STATES(4)*STATES(11))/( CONSTS(350)*CONSTS(348)*ALGBRC(27) ** 2.00000)+( STATES(4)*STATES(12))/( CONSTS(326)*CONSTS(312)*ALGBRC(30) ** 2.00000)+( STATES(4)*STATES(13))/( CONSTS(341)*CONSTS(337)*ALGBRC(32) ** 2.00000)+( STATES(4)*STATES(14))/( CONSTS(356)*CONSTS(354)*ALGBRC(34) ** 2.00000)+( STATES(4)*STATES(22))/( CONSTS(113)*CONSTS(112)*CONSTS(278) ** 2.00000)+( STATES(4)*ALGBRC(114))/( CONSTS(110)*CONSTS(109)*CONSTS(277) ** 2.00000)+( STATES(4)*STATES(28))/( CONSTS(116)*CONSTS(115)*CONSTS(305) ** 2.00000)+( STATES(4)*ALGBRC(119))/( CONSTS(119)*CONSTS(118)*CONSTS(321) ** 2.00000)+( STATES(4)*STATES(18))/( CONSTS(100)*CONSTS(375)*ALGBRC(47) ** 2.00000)+( STATES(4)*STATES(23))/( CONSTS(91)*CONSTS(360)*ALGBRC(37) ** 2.00000)+( STATES(4)*STATES(21))/( CONSTS(122)*CONSTS(121)*CONSTS(279) ** 2.00000)+( STATES(4)*STATES(16))/( CONSTS(363)*CONSTS(362)*ALGBRC(38) ** 2.00000)+( STATES(4)*STATES(25))/( CONSTS(92)*CONSTS(365)*ALGBRC(41) ** 2.00000)+( STATES(4)*STATES(26))/( CONSTS(125)*CONSTS(124)*CONSTS(307) ** 2.00000)+( STATES(4)*STATES(24))/( CONSTS(94)*CONSTS(93)*ALGBRC(42) ** 2.00000)+( STATES(4)*STATES(19))/( CONSTS(89)*CONSTS(88)*ALGBRC(36) ** 2.00000)+( STATES(4)*STATES(27))/( CONSTS(369)*CONSTS(368)*ALGBRC(43) ** 2.00000)+( STATES(4)*STATES(29))/( CONSTS(96)*CONSTS(95)*ALGBRC(44) ** 2.00000)+( STATES(4)*STATES(30))/( CONSTS(373)*CONSTS(372)*ALGBRC(45) ** 2.00000)+( STATES(4)*STATES(32))/( CONSTS(103)*CONSTS(378)*ALGBRC(49) ** 2.00000)+( STATES(4)*STATES(31))/( CONSTS(104)*CONSTS(380)*ALGBRC(50) ** 2.00000)+( STATES(4)*STATES(4))/( CONSTS(143)*CONSTS(142)*CONSTS(285) ** 2.00000)+( STATES(4)*STATES(3))/( CONSTS(140)*CONSTS(139)*CONSTS(284) ** 2.00000)+( STATES(4)*STATES(57))/( CONSTS(128)*CONSTS(127)*CONSTS(280) ** 2.00000)+( STATES(4)*CONSTS(292))/( CONSTS(131)*CONSTS(130)*CONSTS(281) ** 2.00000)+( STATES(4)*CONSTS(311))/( CONSTS(134)*CONSTS(133)*CONSTS(282) ** 2.00000)+( STATES(4)*STATES(20))/( CONSTS(98)*CONSTS(97)*ALGBRC(46) ** 2.00000)) ALGBRC(198) = ( (1.00000+STATES(2)/CONSTS(135)+STATES(3)/CONSTS(136))*STATES(2))/( CONSTS(137)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312))*STATES(1))/( CONSTS(326)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337))*STATES(9))/( CONSTS(341)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(345)+STATES(3)/CONSTS(348))*STATES(11))/( CONSTS(350)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(3)/CONSTS(312))*STATES(12))/( CONSTS(326)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(3)/CONSTS(337))*STATES(13))/( CONSTS(341)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(352)+STATES(3)/CONSTS(354))*STATES(14))/( CONSTS(356)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(111)+STATES(3)/CONSTS(112))*STATES(22))/( CONSTS(113)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(108)+STATES(3)/CONSTS(109))*ALGBRC(114))/( CONSTS(110)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(114)+STATES(3)/CONSTS(115))*STATES(28))/( CONSTS(116)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(117)+STATES(3)/CONSTS(118))*ALGBRC(119))/( CONSTS(119)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(99)+STATES(3)/CONSTS(375))*STATES(18))/( CONSTS(100)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(90)+STATES(3)/CONSTS(360))*STATES(23))/( CONSTS(91)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(120)+STATES(3)/CONSTS(121))*STATES(21))/( CONSTS(122)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(361)+STATES(3)/CONSTS(362))*STATES(16))/( CONSTS(363)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(364)+STATES(3)/CONSTS(365))*STATES(25))/( CONSTS(92)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(123)+STATES(3)/CONSTS(124))*STATES(26))/( CONSTS(125)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(366)+STATES(3)/CONSTS(93))*STATES(24))/( CONSTS(94)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(358)+STATES(3)/CONSTS(88))*STATES(19))/( CONSTS(89)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(367)+STATES(3)/CONSTS(368))*STATES(27))/( CONSTS(369)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(370)+STATES(3)/CONSTS(95))*STATES(29))/( CONSTS(96)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(371)+STATES(3)/CONSTS(372))*STATES(30))/( CONSTS(373)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(377)+STATES(3)/CONSTS(378))*STATES(32))/( CONSTS(103)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(379)+STATES(3)/CONSTS(380))*STATES(31))/( CONSTS(104)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(141)+STATES(3)/CONSTS(142))*STATES(4))/( CONSTS(143)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(138)+STATES(3)/CONSTS(139))*STATES(3))/( CONSTS(140)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(126)+STATES(3)/CONSTS(127))*STATES(57))/( CONSTS(128)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(129)+STATES(3)/CONSTS(130))*CONSTS(292))/( CONSTS(131)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(132)+STATES(3)/CONSTS(133))*CONSTS(311))/( CONSTS(134)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(374)+STATES(3)/CONSTS(97))*STATES(20))/( CONSTS(98)*ALGBRC(46) ** 2.00000) ALGBRC(206) = 1.00000+ALGBRC(198) ALGBRC(195) = ( (1.00000+STATES(2)/CONSTS(135)+STATES(4)/CONSTS(137))*STATES(2))/( CONSTS(136)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(4)/CONSTS(326))*STATES(1))/( CONSTS(312)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(4)/CONSTS(341))*STATES(9))/( CONSTS(337)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(345)+STATES(4)/CONSTS(350))*STATES(11))/( CONSTS(348)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(297)+STATES(4)/CONSTS(326))*STATES(12))/( CONSTS(312)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(333)+STATES(4)/CONSTS(341))*STATES(13))/( CONSTS(337)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(352)+STATES(4)/CONSTS(356))*STATES(14))/( CONSTS(354)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(111)+STATES(4)/CONSTS(113))*STATES(22))/( CONSTS(112)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(108)+STATES(4)/CONSTS(110))*ALGBRC(114))/( CONSTS(109)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(114)+STATES(4)/CONSTS(116))*STATES(28))/( CONSTS(115)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(117)+STATES(4)/CONSTS(119))*ALGBRC(119))/( CONSTS(118)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(99)+STATES(4)/CONSTS(100))*STATES(18))/( CONSTS(375)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(90)+STATES(4)/CONSTS(91))*STATES(23))/( CONSTS(360)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(120)+STATES(4)/CONSTS(122))*STATES(21))/( CONSTS(121)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(361)+STATES(4)/CONSTS(363))*STATES(16))/( CONSTS(362)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(364)+STATES(4)/CONSTS(92))*STATES(25))/( CONSTS(365)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(123)+STATES(4)/CONSTS(125))*STATES(26))/( CONSTS(124)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(366)+STATES(4)/CONSTS(94))*STATES(24))/( CONSTS(93)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(358)+STATES(4)/CONSTS(89))*STATES(19))/( CONSTS(88)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(367)+STATES(4)/CONSTS(369))*STATES(27))/( CONSTS(368)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(370)+STATES(4)/CONSTS(96))*STATES(29))/( CONSTS(95)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(371)+STATES(4)/CONSTS(373))*STATES(30))/( CONSTS(372)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(377)+STATES(4)/CONSTS(103))*STATES(32))/( CONSTS(378)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(379)+STATES(4)/CONSTS(104))*STATES(31))/( CONSTS(380)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(141)+STATES(4)/CONSTS(143))*STATES(4))/( CONSTS(142)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(138)+STATES(4)/CONSTS(140))*STATES(3))/( CONSTS(139)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(126)+STATES(4)/CONSTS(128))*STATES(57))/( CONSTS(127)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(129)+STATES(4)/CONSTS(131))*CONSTS(292))/( CONSTS(130)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(132)+STATES(4)/CONSTS(134))*CONSTS(311))/( CONSTS(133)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(2)/CONSTS(374)+STATES(4)/CONSTS(98))*STATES(20))/( CONSTS(97)*ALGBRC(46) ** 2.00000) ALGBRC(205) = 1.00000+ALGBRC(195) ALGBRC(192) = ( (1.00000+STATES(3)/CONSTS(136)+STATES(4)/CONSTS(137))*STATES(2))/( CONSTS(135)*CONSTS(283) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326))*STATES(1))/( CONSTS(297)*ALGBRC(1) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341))*STATES(9))/( CONSTS(333)*ALGBRC(22) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(348)+STATES(4)/CONSTS(350))*STATES(11))/( CONSTS(345)*ALGBRC(27) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(312)+STATES(4)/CONSTS(326))*STATES(12))/( CONSTS(297)*ALGBRC(30) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(337)+STATES(4)/CONSTS(341))*STATES(13))/( CONSTS(333)*ALGBRC(32) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(354)+STATES(4)/CONSTS(356))*STATES(14))/( CONSTS(352)*ALGBRC(34) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(112)+STATES(4)/CONSTS(113))*STATES(22))/( CONSTS(111)*CONSTS(278) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(109)+STATES(4)/CONSTS(110))*ALGBRC(114))/( CONSTS(108)*CONSTS(277) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(115)+STATES(4)/CONSTS(116))*STATES(28))/( CONSTS(114)*CONSTS(305) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(118)+STATES(4)/CONSTS(119))*ALGBRC(119))/( CONSTS(117)*CONSTS(321) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(375)+STATES(4)/CONSTS(100))*STATES(18))/( CONSTS(99)*ALGBRC(47) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(360)+STATES(4)/CONSTS(91))*STATES(23))/( CONSTS(90)*ALGBRC(37) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(121)+STATES(4)/CONSTS(122))*STATES(21))/( CONSTS(120)*CONSTS(279) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(362)+STATES(4)/CONSTS(363))*STATES(16))/( CONSTS(361)*ALGBRC(38) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(365)+STATES(4)/CONSTS(92))*STATES(25))/( CONSTS(364)*ALGBRC(41) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(124)+STATES(4)/CONSTS(125))*STATES(26))/( CONSTS(123)*CONSTS(307) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(93)+STATES(4)/CONSTS(94))*STATES(24))/( CONSTS(366)*ALGBRC(42) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(88)+STATES(4)/CONSTS(89))*STATES(19))/( CONSTS(358)*ALGBRC(36) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(368)+STATES(4)/CONSTS(369))*STATES(27))/( CONSTS(367)*ALGBRC(43) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(95)+STATES(4)/CONSTS(96))*STATES(29))/( CONSTS(370)*ALGBRC(44) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(372)+STATES(4)/CONSTS(373))*STATES(30))/( CONSTS(371)*ALGBRC(45) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(378)+STATES(4)/CONSTS(103))*STATES(32))/( CONSTS(377)*ALGBRC(49) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(380)+STATES(4)/CONSTS(104))*STATES(31))/( CONSTS(379)*ALGBRC(50) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(142)+STATES(4)/CONSTS(143))*STATES(4))/( CONSTS(141)*CONSTS(285) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(139)+STATES(4)/CONSTS(140))*STATES(3))/( CONSTS(138)*CONSTS(284) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(127)+STATES(4)/CONSTS(128))*STATES(57))/( CONSTS(126)*CONSTS(280) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(130)+STATES(4)/CONSTS(131))*CONSTS(292))/( CONSTS(129)*CONSTS(281) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(133)+STATES(4)/CONSTS(134))*CONSTS(311))/( CONSTS(132)*CONSTS(282) ** 2.00000)+( (1.00000+STATES(3)/CONSTS(97)+STATES(4)/CONSTS(98))*STATES(20))/( CONSTS(374)*ALGBRC(46) ** 2.00000) ALGBRC(204) = 1.00000+ALGBRC(192)+CONSTS(272)/( CONSTS(273)*1.00000+STATES(2)/CONSTS(273) ** 2.00000) ALGBRC(196) = - (( STATES(4)*STATES(2))/( CONSTS(137)*CONSTS(135)*CONSTS(283) ** 2.00000)+( STATES(4)*STATES(1))/( CONSTS(326)*CONSTS(297)*ALGBRC(1) ** 2.00000)+( STATES(4)*STATES(9))/( CONSTS(341)*CONSTS(333)*ALGBRC(22) ** 2.00000)+( STATES(4)*STATES(11))/( CONSTS(350)*CONSTS(345)*ALGBRC(27) ** 2.00000)+( STATES(4)*STATES(12))/( CONSTS(326)*CONSTS(297)*ALGBRC(30) ** 2.00000)+( STATES(4)*STATES(13))/( CONSTS(341)*CONSTS(333)*ALGBRC(32) ** 2.00000)+( STATES(4)*STATES(14))/( CONSTS(356)*CONSTS(352)*ALGBRC(34) ** 2.00000)+( STATES(4)*STATES(22))/( CONSTS(113)*CONSTS(111)*CONSTS(278) ** 2.00000)+( STATES(4)*ALGBRC(114))/( CONSTS(110)*CONSTS(108)*CONSTS(277) ** 2.00000)+( STATES(4)*STATES(28))/( CONSTS(116)*CONSTS(114)*CONSTS(305) ** 2.00000)+( STATES(4)*ALGBRC(119))/( CONSTS(119)*CONSTS(117)*CONSTS(321) ** 2.00000)+( STATES(4)*STATES(18))/( CONSTS(100)*CONSTS(99)*ALGBRC(47) ** 2.00000)+( STATES(4)*STATES(23))/( CONSTS(91)*CONSTS(90)*ALGBRC(37) ** 2.00000)+( STATES(4)*STATES(21))/( CONSTS(122)*CONSTS(120)*CONSTS(279) ** 2.00000)+( STATES(4)*STATES(16))/( CONSTS(363)*CONSTS(361)*ALGBRC(38) ** 2.00000)+( STATES(4)*STATES(25))/( CONSTS(92)*CONSTS(364)*ALGBRC(41) ** 2.00000)+( STATES(4)*STATES(26))/( CONSTS(125)*CONSTS(123)*CONSTS(307) ** 2.00000)+( STATES(4)*STATES(24))/( CONSTS(94)*CONSTS(366)*ALGBRC(42) ** 2.00000)+( STATES(4)*STATES(19))/( CONSTS(89)*CONSTS(358)*ALGBRC(36) ** 2.00000)+( STATES(4)*STATES(27))/( CONSTS(369)*CONSTS(367)*ALGBRC(43) ** 2.00000)+( STATES(4)*STATES(29))/( CONSTS(96)*CONSTS(370)*ALGBRC(44) ** 2.00000)+( STATES(4)*STATES(30))/( CONSTS(373)*CONSTS(371)*ALGBRC(45) ** 2.00000)+( STATES(4)*STATES(32))/( CONSTS(103)*CONSTS(377)*ALGBRC(49) ** 2.00000)+( STATES(4)*STATES(31))/( CONSTS(104)*CONSTS(379)*ALGBRC(50) ** 2.00000)+( STATES(4)*STATES(4))/( CONSTS(143)*CONSTS(141)*CONSTS(285) ** 2.00000)+( STATES(4)*STATES(3))/( CONSTS(140)*CONSTS(138)*CONSTS(284) ** 2.00000)+( STATES(4)*STATES(57))/( CONSTS(128)*CONSTS(126)*CONSTS(280) ** 2.00000)+( STATES(4)*CONSTS(292))/( CONSTS(131)*CONSTS(129)*CONSTS(281) ** 2.00000)+( STATES(4)*CONSTS(311))/( CONSTS(134)*CONSTS(132)*CONSTS(282) ** 2.00000)+( STATES(4)*STATES(20))/( CONSTS(98)*CONSTS(374)*ALGBRC(46) ** 2.00000)) ALGBRC(193) = - (( STATES(3)*STATES(2))/( CONSTS(136)*CONSTS(135)*CONSTS(283) ** 2.00000)+( STATES(3)*STATES(1))/( CONSTS(312)*CONSTS(297)*ALGBRC(1) ** 2.00000)+( STATES(3)*STATES(9))/( CONSTS(337)*CONSTS(333)*ALGBRC(22) ** 2.00000)+( STATES(3)*STATES(11))/( CONSTS(348)*CONSTS(345)*ALGBRC(27) ** 2.00000)+( STATES(3)*STATES(12))/( CONSTS(312)*CONSTS(297)*ALGBRC(30) ** 2.00000)+( STATES(3)*STATES(13))/( CONSTS(337)*CONSTS(333)*ALGBRC(32) ** 2.00000)+( STATES(3)*STATES(14))/( CONSTS(354)*CONSTS(352)*ALGBRC(34) ** 2.00000)+( STATES(3)*STATES(22))/( CONSTS(112)*CONSTS(111)*CONSTS(278) ** 2.00000)+( STATES(3)*ALGBRC(114))/( CONSTS(109)*CONSTS(108)*CONSTS(277) ** 2.00000)+( STATES(3)*STATES(28))/( CONSTS(115)*CONSTS(114)*CONSTS(305) ** 2.00000)+( STATES(3)*ALGBRC(119))/( CONSTS(118)*CONSTS(117)*CONSTS(321) ** 2.00000)+( STATES(3)*STATES(18))/( CONSTS(375)*CONSTS(99)*ALGBRC(47) ** 2.00000)+( STATES(3)*STATES(23))/( CONSTS(360)*CONSTS(90)*ALGBRC(37) ** 2.00000)+( STATES(3)*STATES(21))/( CONSTS(121)*CONSTS(120)*CONSTS(279) ** 2.00000)+( STATES(3)*STATES(16))/( CONSTS(362)*CONSTS(361)*ALGBRC(38) ** 2.00000)+( STATES(3)*STATES(25))/( CONSTS(365)*CONSTS(364)*ALGBRC(41) ** 2.00000)+( STATES(3)*STATES(26))/( CONSTS(124)*CONSTS(123)*CONSTS(307) ** 2.00000)+( STATES(3)*STATES(24))/( CONSTS(93)*CONSTS(366)*ALGBRC(42) ** 2.00000)+( STATES(3)*STATES(19))/( CONSTS(88)*CONSTS(358)*ALGBRC(36) ** 2.00000)+( STATES(3)*STATES(27))/( CONSTS(368)*CONSTS(367)*ALGBRC(43) ** 2.00000)+( STATES(3)*STATES(29))/( CONSTS(95)*CONSTS(370)*ALGBRC(44) ** 2.00000)+( STATES(3)*STATES(30))/( CONSTS(372)*CONSTS(371)*ALGBRC(45) ** 2.00000)+( STATES(3)*STATES(32))/( CONSTS(378)*CONSTS(377)*ALGBRC(49) ** 2.00000)+( STATES(3)*STATES(31))/( CONSTS(380)*CONSTS(379)*ALGBRC(50) ** 2.00000)+( STATES(3)*STATES(4))/( CONSTS(142)*CONSTS(141)*CONSTS(285) ** 2.00000)+( STATES(3)*STATES(3))/( CONSTS(139)*CONSTS(138)*CONSTS(284) ** 2.00000)+( STATES(3)*STATES(57))/( CONSTS(127)*CONSTS(126)*CONSTS(280) ** 2.00000)+( STATES(3)*CONSTS(292))/( CONSTS(130)*CONSTS(129)*CONSTS(281) ** 2.00000)+( STATES(3)*CONSTS(311))/( CONSTS(133)*CONSTS(132)*CONSTS(282) ** 2.00000)+( STATES(3)*STATES(20))/( CONSTS(97)*CONSTS(374)*ALGBRC(46) ** 2.00000)) ALGBRC(207) = ((( ALGBRC(204)*ALGBRC(197)*ALGBRC(194)+ ALGBRC(206)*ALGBRC(190)*ALGBRC(193)+ ALGBRC(205)*ALGBRC(191)*ALGBRC(196)) - ALGBRC(205)*ALGBRC(206)*ALGBRC(204)) - ALGBRC(191)*ALGBRC(197)*ALGBRC(193)) - ALGBRC(190)*ALGBRC(194)*ALGBRC(196) ALGBRC(146) = ((ALGBRC(115)+ALGBRC(116)+ALGBRC(117)+ALGBRC(118)) - ALGBRC(140))/CONSTS(4) ALGBRC(147) = - ALGBRC(146) ALGBRC(148) = ((ALGBRC(120)+ALGBRC(140)) - ALGBRC(123))/CONSTS(4) ALGBRC(149) = - ALGBRC(148) ALGBRC(184) = ((ALGBRC(82)+ALGBRC(144)) - ALGBRC(179))/CONSTS(4) ALGBRC(185) = ((- ALGBRC(82) - ALGBRC(144))+ALGBRC(179))/CONSTS(4) ALGBRC(124) = (ALGBRC(67) - ALGBRC(82))/CONSTS(4) ALGBRC(125) = (- ALGBRC(67)+ALGBRC(82))/CONSTS(4) ALGBRC(186) = (((- ALGBRC(67) - ALGBRC(144))+ALGBRC(180)) - ALGBRC(177))/CONSTS(4) ALGBRC(160) = (- ALGBRC(115)+ALGBRC(156))/CONSTS(4) ALGBRC(126) = (ALGBRC(115) - ALGBRC(57))/CONSTS(4) ALGBRC(187) = (- ALGBRC(60)+ALGBRC(176)+ALGBRC(57))/CONSTS(4) ALGBRC(127) = (ALGBRC(60) - ALGBRC(116))/CONSTS(4) ALGBRC(128) = (((ALGBRC(116) - ALGBRC(117)) - ALGBRC(87))+ALGBRC(88))/CONSTS(4) ALGBRC(129) = (ALGBRC(117) - ALGBRC(67))/CONSTS(4) ALGBRC(130) = ((- ALGBRC(115) - ALGBRC(117))+ALGBRC(67)+ALGBRC(57))/CONSTS(4) ALGBRC(188) = ((ALGBRC(67) - ALGBRC(120))+ALGBRC(178))/CONSTS(4) ALGBRC(131) = (ALGBRC(120) - ALGBRC(74))/CONSTS(4) ALGBRC(189) = (((((ALGBRC(74) - ALGBRC(118))+ALGBRC(177)) - ALGBRC(88)) - ALGBRC(176)) - ALGBRC(178))/CONSTS(4) ALGBRC(132) = (- ALGBRC(57)+ALGBRC(118)+ALGBRC(87))/CONSTS(4) ALGBRC(161) = ((ALGBRC(87)+ALGBRC(157)) - ALGBRC(159))/CONSTS(4) ALGBRC(162) = (- ALGBRC(87)+ALGBRC(159))/CONSTS(4) ALGBRC(182) = CONSTS(51)*( ALGBRC(163)*STATES(2) - STATES(4)*ALGBRC(133)) ALGBRC(199) = (((((((((- ALGBRC(115)+ 2.00000*ALGBRC(57)) - ALGBRC(117))+ALGBRC(67)+ALGBRC(118)+ALGBRC(156)+ALGBRC(157)+ALGBRC(176)) - ALGBRC(159)) - 5.00000*ALGBRC(140)) - 2.00000*ALGBRC(123)) - 4.00000*ALGBRC(142))+ (CONSTS(8) - 1.00000)*ALGBRC(144)+ 2.00000*ALGBRC(180)+ALGBRC(145)) - ALGBRC(182))/CONSTS(4) ALGBRC(201) = - (( STATES(2)*CONSTS(293))/( CONSTS(135)*CONSTS(283))+( STATES(2)*ALGBRC(184))/( CONSTS(297)*ALGBRC(1))+( STATES(2)*ALGBRC(185))/( CONSTS(333)*ALGBRC(22))+( STATES(2)*CONSTS(263))/( CONSTS(345)*ALGBRC(27))+( STATES(2)*ALGBRC(124))/( CONSTS(297)*ALGBRC(30))+( STATES(2)*ALGBRC(125))/( CONSTS(333)*ALGBRC(32))+( STATES(2)*ALGBRC(186))/( CONSTS(352)*ALGBRC(34))+( STATES(2)*ALGBRC(146))/( CONSTS(111)*CONSTS(278))+( STATES(2)*ALGBRC(147))/( CONSTS(108)*CONSTS(277))+( STATES(2)*ALGBRC(148))/( CONSTS(114)*CONSTS(305))+( STATES(2)*ALGBRC(149))/( CONSTS(117)*CONSTS(321))+( STATES(2)*ALGBRC(160))/( CONSTS(99)*ALGBRC(47))+( STATES(2)*ALGBRC(132))/( CONSTS(90)*ALGBRC(37))+( STATES(2)*ALGBRC(126))/( CONSTS(120)*CONSTS(279))+( STATES(2)*ALGBRC(187))/( CONSTS(361)*ALGBRC(38))+( STATES(2)*ALGBRC(127))/( CONSTS(364)*ALGBRC(41))+( STATES(2)*ALGBRC(128))/( CONSTS(123)*CONSTS(307))+( STATES(2)*ALGBRC(129))/( CONSTS(366)*ALGBRC(42))+( STATES(2)*ALGBRC(130))/( CONSTS(358)*ALGBRC(36))+( STATES(2)*ALGBRC(188))/( CONSTS(367)*ALGBRC(43))+( STATES(2)*ALGBRC(131))/( CONSTS(370)*ALGBRC(44))+( STATES(2)*ALGBRC(189))/( CONSTS(371)*ALGBRC(45))+( STATES(2)*ALGBRC(161))/( CONSTS(377)*ALGBRC(49))+( STATES(2)*ALGBRC(162))/( CONSTS(379)*ALGBRC(50))+( STATES(2)*CONSTS(295))/( CONSTS(141)*CONSTS(285))+( STATES(2)*CONSTS(294))/( CONSTS(138)*CONSTS(284))+( STATES(2)*CONSTS(264))/( CONSTS(126)*CONSTS(280))+( STATES(2)*CONSTS(271))/( CONSTS(129)*CONSTS(281))+( STATES(2)*CONSTS(270))/( CONSTS(132)*CONSTS(282))+( STATES(2)*CONSTS(265))/( CONSTS(374)*ALGBRC(46)))+ALGBRC(199) ALGBRC(202) = - (( STATES(3)*CONSTS(293))/( CONSTS(136)*CONSTS(283))+( STATES(3)*ALGBRC(184))/( CONSTS(312)*ALGBRC(1))+( STATES(3)*ALGBRC(185))/( CONSTS(337)*ALGBRC(22))+( STATES(3)*CONSTS(263))/( CONSTS(348)*ALGBRC(27))+( STATES(3)*ALGBRC(124))/( CONSTS(312)*ALGBRC(30))+( STATES(3)*ALGBRC(125))/( CONSTS(337)*ALGBRC(32))+( STATES(3)*ALGBRC(186))/( CONSTS(354)*ALGBRC(34))+( STATES(3)*ALGBRC(146))/( CONSTS(112)*CONSTS(278))+( STATES(3)*ALGBRC(147))/( CONSTS(109)*CONSTS(277))+( STATES(3)*ALGBRC(148))/( CONSTS(115)*CONSTS(305))+( STATES(3)*ALGBRC(149))/( CONSTS(118)*CONSTS(321))+( STATES(3)*ALGBRC(160))/( CONSTS(375)*ALGBRC(47))+( STATES(3)*ALGBRC(132))/( CONSTS(360)*ALGBRC(37))+( STATES(3)*ALGBRC(126))/( CONSTS(121)*CONSTS(279))+( STATES(3)*ALGBRC(187))/( CONSTS(362)*ALGBRC(38))+( STATES(3)*ALGBRC(127))/( CONSTS(365)*ALGBRC(41))+( STATES(3)*ALGBRC(128))/( CONSTS(124)*CONSTS(307))+( STATES(3)*ALGBRC(129))/( CONSTS(93)*ALGBRC(42))+( STATES(3)*ALGBRC(130))/( CONSTS(88)*ALGBRC(36))+( STATES(3)*ALGBRC(188))/( CONSTS(368)*ALGBRC(43))+( STATES(3)*ALGBRC(131))/( CONSTS(95)*ALGBRC(44))+( STATES(3)*ALGBRC(189))/( CONSTS(372)*ALGBRC(45))+( STATES(3)*ALGBRC(161))/( CONSTS(378)*ALGBRC(49))+( STATES(3)*ALGBRC(162))/( CONSTS(380)*ALGBRC(50))+( STATES(3)*CONSTS(295))/( CONSTS(142)*CONSTS(285))+( STATES(3)*CONSTS(294))/( CONSTS(139)*CONSTS(284))+( STATES(3)*CONSTS(264))/( CONSTS(127)*CONSTS(280))+( STATES(3)*CONSTS(271))/( CONSTS(130)*CONSTS(281))+( STATES(3)*CONSTS(270))/( CONSTS(133)*CONSTS(282))+( STATES(3)*CONSTS(265))/( CONSTS(97)*ALGBRC(46)))+CONSTS(296) ALGBRC(200) = ALGBRC(182)/CONSTS(4) ALGBRC(203) = - (( STATES(4)*CONSTS(293))/( CONSTS(137)*CONSTS(283))+( STATES(4)*ALGBRC(184))/( CONSTS(326)*ALGBRC(1))+( STATES(4)*ALGBRC(185))/( CONSTS(341)*ALGBRC(22))+( STATES(4)*CONSTS(263))/( CONSTS(350)*ALGBRC(27))+( STATES(4)*ALGBRC(124))/( CONSTS(326)*ALGBRC(30))+( STATES(4)*ALGBRC(125))/( CONSTS(341)*ALGBRC(32))+( STATES(4)*ALGBRC(186))/( CONSTS(356)*ALGBRC(34))+( STATES(4)*ALGBRC(146))/( CONSTS(113)*CONSTS(278))+( STATES(4)*ALGBRC(147))/( CONSTS(110)*CONSTS(277))+( STATES(4)*ALGBRC(148))/( CONSTS(116)*CONSTS(305))+( STATES(4)*ALGBRC(149))/( CONSTS(119)*CONSTS(321))+( STATES(4)*ALGBRC(160))/( CONSTS(100)*ALGBRC(47))+( STATES(4)*ALGBRC(132))/( CONSTS(91)*ALGBRC(37))+( STATES(4)*ALGBRC(126))/( CONSTS(122)*CONSTS(279))+( STATES(4)*ALGBRC(187))/( CONSTS(363)*ALGBRC(38))+( STATES(4)*ALGBRC(127))/( CONSTS(92)*ALGBRC(41))+( STATES(4)*ALGBRC(128))/( CONSTS(125)*CONSTS(307))+( STATES(4)*ALGBRC(129))/( CONSTS(94)*ALGBRC(42))+( STATES(4)*ALGBRC(130))/( CONSTS(89)*ALGBRC(36))+( STATES(4)*ALGBRC(188))/( CONSTS(369)*ALGBRC(43))+( STATES(4)*ALGBRC(131))/( CONSTS(96)*ALGBRC(44))+( STATES(4)*ALGBRC(189))/( CONSTS(373)*ALGBRC(45))+( STATES(4)*ALGBRC(161))/( CONSTS(103)*ALGBRC(49))+( STATES(4)*ALGBRC(162))/( CONSTS(104)*ALGBRC(50))+( STATES(4)*CONSTS(295))/( CONSTS(143)*CONSTS(285))+( STATES(4)*CONSTS(294))/( CONSTS(140)*CONSTS(284))+( STATES(4)*CONSTS(264))/( CONSTS(128)*CONSTS(280))+( STATES(4)*CONSTS(271))/( CONSTS(131)*CONSTS(281))+( STATES(4)*CONSTS(270))/( CONSTS(134)*CONSTS(282))+( STATES(4)*CONSTS(265))/( CONSTS(98)*ALGBRC(46)))+ALGBRC(200) ALGBRC(2) = 1.00000+STATES(6)/CONSTS(352)+STATES(7)/CONSTS(354)+STATES(8)/CONSTS(356) ALGBRC(3) = 1.00000+STATES(6)/CONSTS(358) ALGBRC(4) = 1.00000+STATES(7)/CONSTS(360) ALGBRC(5) = 1.00000+STATES(6)/CONSTS(361)+STATES(7)/CONSTS(362)+STATES(8)/CONSTS(363) ALGBRC(6) = 1.00000+STATES(6)/CONSTS(364)+STATES(7)/CONSTS(365) ALGBRC(7) = 1.00000+STATES(6)/CONSTS(366) ALGBRC(8) = 1.00000+STATES(6)/CONSTS(367)+STATES(7)/CONSTS(368)+STATES(8)/CONSTS(369) ALGBRC(9) = 1.00000+STATES(6)/CONSTS(370) ALGBRC(10) = 1.00000+STATES(6)/CONSTS(371)+STATES(7)/CONSTS(372)+STATES(8)/CONSTS(373) ALGBRC(11) = 1.00000+STATES(6)/CONSTS(374) ALGBRC(12) = 1.00000+STATES(7)/CONSTS(375) ALGBRC(13) = 1.00000+STATES(2)/CONSTS(376) ALGBRC(14) = 1.00000+STATES(6)/CONSTS(377)+STATES(7)/CONSTS(378) ALGBRC(15) = 1.00000+STATES(6)/CONSTS(379)+STATES(7)/CONSTS(380) ALGBRC(20) = ALGBRC(18) ALGBRC(25) = ALGBRC(23) ALGBRC(134) = 1.00000+ALGBRC(133)/CONSTS(358) ALGBRC(135) = 1.00000+ALGBRC(133)/CONSTS(366) ALGBRC(136) = 1.00000+ALGBRC(133)/CONSTS(370) ALGBRC(137) = 1.00000+ALGBRC(133)/CONSTS(374) ALGBRC(138) = 1.00000+ALGBRC(133)/CONSTS(376) ALGBRC(151) = 1.00000+ALGBRC(150)/CONSTS(360) ALGBRC(152) = 1.00000+ALGBRC(133)/CONSTS(364)+ALGBRC(150)/CONSTS(365) ALGBRC(165) = ALGBRC(164) ALGBRC(167) = ALGBRC(166) ALGBRC(172) = ( STATES(17)*(1.00000+ALGBRC(133)/CONSTS(361)))/ALGBRC(171) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END