Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 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