C
C There are a total of 164 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 139 entries in the constant variable array.
C
C
C CONSTS(1) is kappa_1a in component BG_parameters (fmol_per_sec).
C CONSTS(2) is kappa_1b in component BG_parameters (fmol_per_sec).
C CONSTS(3) is kappa_2a in component BG_parameters (fmol_per_sec).
C CONSTS(4) is kappa_2b in component BG_parameters (fmol_per_sec).
C CONSTS(5) is kappa_3a in component BG_parameters (fmol_per_sec).
C CONSTS(6) is kappa_3b in component BG_parameters (fmol_per_sec).
C CONSTS(7) is kappa_4a in component BG_parameters (fmol_per_sec).
C CONSTS(8) is kappa_4b in component BG_parameters (fmol_per_sec).
C CONSTS(9) is kappa_5 in component BG_parameters (fmol_per_sec).
C CONSTS(10) is kappa_6 in component BG_parameters (fmol_per_sec).
C CONSTS(11) is kappa_7 in component BG_parameters (fmol_per_sec).
C CONSTS(12) is kappa_GiAC in component BG_parameters (fmol_per_sec).
C CONSTS(13) is kappa_sig1_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(14) is kappa_sig2_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(15) is kappa_sig3_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(16) is kappa_sig4_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(17) is kappa_sig1_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(18) is kappa_sig2_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(19) is kappa_sig3_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(20) is kappa_sig4_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(21) is kappa_Doff1_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(22) is kappa_Ton1_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(23) is kappa_Act1_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(24) is kappa_Doff2_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(25) is kappa_Ton2_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(26) is kappa_Act2_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(27) is kappa_hyd_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(28) is kappa_reassoc_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(29) is kappa_Doff1_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(30) is kappa_Ton1_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(31) is kappa_Act1_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(32) is kappa_Doff2_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(33) is kappa_Ton2_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(34) is kappa_Act2_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(35) is kappa_hyd_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(36) is kappa_reassoc_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(37) is K_ATP in component BG_parameters (per_fmol).
C CONSTS(38) is K_cAMP in component BG_parameters (per_fmol).
C CONSTS(39) is K_AC in component BG_parameters (per_fmol).
C CONSTS(40) is K_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(41) is K_a_Gs_GTP_AC in component BG_parameters (per_fmol).
C CONSTS(42) is K_a_Gs_GTP_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(43) is K_FSK_AC in component BG_parameters (per_fmol).
C CONSTS(44) is K_FSK_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(45) is K_PDE in component BG_parameters (per_fmol).
C CONSTS(46) is K_PDE_cAMP in component BG_parameters (per_fmol).
C CONSTS(47) is K_five_AMP in component BG_parameters (per_fmol).
C CONSTS(48) is K_IBMX in component BG_parameters (per_fmol).
C CONSTS(49) is K_PDEinh in component BG_parameters (per_fmol).
C CONSTS(50) is K_a_Gs_GTP in component BG_parameters (per_fmol).
C CONSTS(51) is K_FSK in component BG_parameters (per_fmol).
C CONSTS(52) is K_a_Gi_GTP in component BG_parameters (per_fmol).
C CONSTS(53) is K_ACinh in component BG_parameters (per_fmol).
C CONSTS(54) is K_PPi in component BG_parameters (per_fmol).
C CONSTS(55) is K_L_B1 in component BG_parameters (per_fmol).
C CONSTS(56) is K_R_B1 in component BG_parameters (per_fmol).
C CONSTS(57) is K_Gs in component BG_parameters (per_fmol).
C CONSTS(58) is K_LR_B1 in component BG_parameters (per_fmol).
C CONSTS(59) is K_R_B1Gs in component BG_parameters (per_fmol).
C CONSTS(60) is K_LR_B1Gs in component BG_parameters (per_fmol).
C CONSTS(61) is K_L_M2 in component BG_parameters (per_fmol).
C CONSTS(62) is K_R_M2 in component BG_parameters (per_fmol).
C CONSTS(63) is K_Gi in component BG_parameters (per_fmol).
C CONSTS(64) is K_LR_M2 in component BG_parameters (per_fmol).
C CONSTS(65) is K_R_M2Gi in component BG_parameters (per_fmol).
C CONSTS(66) is K_LR_M2Gi in component BG_parameters (per_fmol).
C CONSTS(67) is K_beta_gamma_Gs in component BG_parameters (per_fmol).
C CONSTS(68) is K_a_Gs_GDP in component BG_parameters (per_fmol).
C CONSTS(69) is K_Pi in component BG_parameters (per_fmol).
C CONSTS(70) is K_R_B1_aby in component BG_parameters (per_fmol).
C CONSTS(71) is K_R_B1_aby_T in component BG_parameters (per_fmol).
C CONSTS(72) is K_LR_B1_aby in component BG_parameters (per_fmol).
C CONSTS(73) is K_LR_B1_aby_T in component BG_parameters (per_fmol).
C CONSTS(74) is K_GTP in component BG_parameters (per_fmol).
C CONSTS(75) is K_GDP in component BG_parameters (per_fmol).
C CONSTS(76) is K_beta_gamma_Gi in component BG_parameters (per_fmol).
C CONSTS(77) is K_a_Gi_GDP in component BG_parameters (per_fmol).
C CONSTS(78) is K_R_M2_aby in component BG_parameters (per_fmol).
C CONSTS(79) is K_R_M2_aby_T in component BG_parameters (per_fmol).
C CONSTS(80) is K_LR_M2_aby in component BG_parameters (per_fmol).
C CONSTS(81) is K_LR_M2_aby_T in component BG_parameters (per_fmol).
C VOI is time in component environment (second).
C CONSTS(82) is vol_myo in component environment (pL).
C CONSTS(83) is freq in component environment (dimensionless).
C CONSTS(84) is stimSt in component environment (second).
C CONSTS(85) is stimSt2 in component environment (second).
C CONSTS(86) is stimDur in component environment (second).
C CONSTS(87) is tRamp in component environment (second).
C CONSTS(88) is stimMag in component environment (fmol).
C CONSTS(89) is stimHolding in component environment (fmol).
C CONSTS(139) is m in component environment (fmol_per_sec).
C CONSTS(90) is q_ATP_init in component environment (fmol).
C CONSTS(91) is q_AC_init in component environment (fmol).
C CONSTS(92) is q_cAMP_init in component environment (fmol).
C CONSTS(93) is q_AC_ATP_init in component environment (fmol).
C CONSTS(94) is q_FSK_init in component environment (fmol).
C CONSTS(95) is q_FSK_AC_init in component environment (fmol).
C CONSTS(96) is q_FSK_AC_ATP_init in component environment (fmol).
C CONSTS(97) is q_a_Gs_GTP_init in component environment (fmol).
C CONSTS(98) is q_a_Gs_GTP_AC_init in component environment (fmol).
C CONSTS(99) is q_a_Gs_GTP_AC_ATP_init in component environment (fmol).
C CONSTS(100) is q_PDE_init in component environment (fmol).
C CONSTS(101) is q_PDEinh_init in component environment (fmol).
C CONSTS(102) is q_PDE_cAMP_init in component environment (fmol).
C CONSTS(103) is q_IBMX_init in component environment (fmol).
C CONSTS(104) is q_five_AMP_init in component environment (fmol).
C CONSTS(105) is q_a_Gi_GTP_init in component environment (fmol).
C CONSTS(106) is q_ACinh_init in component environment (fmol).
C CONSTS(107) is q_PPi_init in component environment (fmol).
C CONSTS(108) is q_R_B1_init in component environment (fmol).
C CONSTS(109) is q_Gs_init in component environment (fmol).
C ALGBRC(1) is q_L_B1_init in component environment (fmol).
C CONSTS(110) is q_LR_B1_init in component environment (fmol).
C CONSTS(111) is q_R_B1Gs_init in component environment (fmol).
C CONSTS(112) is q_LR_B1Gs_init in component environment (fmol).
C ALGBRC(2) is q_L_M2_init in component environment (fmol).
C CONSTS(113) is q_R_M2_init in component environment (fmol).
C CONSTS(114) is q_Gi_init in component environment (fmol).
C CONSTS(115) is q_LR_M2_init in component environment (fmol).
C CONSTS(116) is q_R_M2Gi_init in component environment (fmol).
C CONSTS(117) is q_LR_M2Gi_init in component environment (fmol).
C CONSTS(118) is q_beta_gamma_Gs_init in component environment (fmol).
C CONSTS(119) is q_a_Gs_GDP_init in component environment (fmol).
C CONSTS(120) is q_Pi_init in component environment (fmol).
C CONSTS(121) is q_R_B1_aby_init in component environment (fmol).
C CONSTS(122) is q_R_B1_aby_T_init in component environment (fmol).
C CONSTS(123) is q_LR_B1_aby_init in component environment (fmol).
C CONSTS(124) is q_LR_B1_aby_T_init in component environment (fmol).
C CONSTS(125) is q_GTP_init in component environment (fmol).
C CONSTS(126) is q_GDP_init in component environment (fmol).
C CONSTS(127) is q_beta_gamma_Gi_init in component environment (fmol).
C CONSTS(128) is q_a_Gi_GDP_init in component environment (fmol).
C CONSTS(129) is q_R_M2_aby_init in component environment (fmol).
C CONSTS(130) is q_R_M2_aby_T_init in component environment (fmol).
C CONSTS(131) is q_LR_M2_aby_init in component environment (fmol).
C CONSTS(132) is q_LR_M2_aby_T_init in component environment (fmol).
C ALGBRC(85) is L_B1_T in component environment (fmol).
C ALGBRC(116) is L_M2_T in component environment (fmol).
C ALGBRC(86) is R_B1_T in component environment (fmol).
C ALGBRC(117) is R_M2_T in component environment (fmol).
C ALGBRC(87) is Gs_T in component environment (fmol).
C ALGBRC(118) is Gi_T in component environment (fmol).
C ALGBRC(17) is adenosine_T in component environment (fmol).
C ALGBRC(3) is q_ATP in component environment (fmol).
C ALGBRC(4) is q_cAMP in component environment (fmol).
C ALGBRC(5) is q_AC in component environment (fmol).
C ALGBRC(6) is q_AC_ATP in component environment (fmol).
C ALGBRC(7) is q_a_Gs_GTP_AC in component environment (fmol).
C ALGBRC(8) is q_a_Gs_GTP_AC_ATP in component environment (fmol).
C ALGBRC(9) is q_FSK_AC in component environment (fmol).
C ALGBRC(10) is q_FSK_AC_ATP in component environment (fmol).
C ALGBRC(11) is q_PDE in component environment (fmol).
C ALGBRC(14) is q_PDE_cAMP in component environment (fmol).
C ALGBRC(15) is q_five_AMP in component environment (fmol).
C ALGBRC(18) is q_IBMX in component environment (fmol).
C ALGBRC(19) is q_PDEinh in component environment (fmol).
C ALGBRC(20) is q_a_Gs_GTP in component environment (fmol).
C ALGBRC(21) is q_FSK in component environment (fmol).
C ALGBRC(23) is q_a_Gi_GTP in component environment (fmol).
C ALGBRC(24) is q_ACinh in component environment (fmol).
C ALGBRC(26) is q_PPi in component environment (fmol).
C ALGBRC(25) is q_L_B1 in component environment (fmol).
C ALGBRC(27) is q_R_B1 in component environment (fmol).
C ALGBRC(28) is q_Gs in component environment (fmol).
C ALGBRC(30) is q_LR_B1 in component environment (fmol).
C ALGBRC(32) is q_R_B1Gs in component environment (fmol).
C ALGBRC(34) is q_LR_B1Gs in component environment (fmol).
C ALGBRC(36) is q_L_M2 in component environment (fmol).
C ALGBRC(42) is q_R_M2 in component environment (fmol).
C ALGBRC(45) is q_Gi in component environment (fmol).
C ALGBRC(48) is q_LR_M2 in component environment (fmol).
C ALGBRC(51) is q_R_M2Gi in component environment (fmol).
C ALGBRC(54) is q_LR_M2Gi in component environment (fmol).
C ALGBRC(57) is q_beta_gamma_Gs in component environment (fmol).
C ALGBRC(64) is q_a_Gs_GDP in component environment (fmol).
C ALGBRC(68) is q_Pi in component environment (fmol).
C ALGBRC(72) is q_R_B1_aby in component environment (fmol).
C ALGBRC(76) is q_R_B1_aby_T in component environment (fmol).
C ALGBRC(79) is q_LR_B1_aby in component environment (fmol).
C ALGBRC(82) is q_LR_B1_aby_T in component environment (fmol).
C ALGBRC(88) is q_GTP in component environment (fmol).
C ALGBRC(93) is q_GDP in component environment (fmol).
C ALGBRC(96) is q_beta_gamma_Gi in component environment (fmol).
C ALGBRC(100) is q_a_Gi_GDP in component environment (fmol).
C ALGBRC(103) is q_R_M2_aby in component environment (fmol).
C ALGBRC(106) is q_R_M2_aby_T in component environment (fmol).
C ALGBRC(109) is q_LR_M2_aby in component environment (fmol).
C ALGBRC(115) is q_LR_M2_aby_T in component environment (fmol).
C STATES(1) is q_ATP in component cAMP (fmol).
C STATES(2) is q_cAMP in component cAMP (fmol).
C STATES(3) is q_AC in component cAMP (fmol).
C STATES(4) is q_AC_ATP in component cAMP (fmol).
C STATES(5) is q_a_Gs_GTP_AC in component cAMP (fmol).
C STATES(6) is q_a_Gs_GTP_AC_ATP in component cAMP (fmol).
C STATES(7) is q_FSK_AC in component cAMP (fmol).
C STATES(8) is q_FSK_AC_ATP in component cAMP (fmol).
C STATES(9) is q_PDE in component cAMP (fmol).
C STATES(10) is q_PDE_cAMP in component cAMP (fmol).
C STATES(11) is q_five_AMP in component cAMP (fmol).
C STATES(12) is q_IBMX in component cAMP (fmol).
C STATES(13) is q_PDEinh in component cAMP (fmol).
C STATES(14) is q_a_Gs_GTP in component cAMP (fmol).
C STATES(15) is q_FSK in component cAMP (fmol).
C STATES(16) is q_a_Gi_GTP in component cAMP (fmol).
C STATES(17) is q_ACinh in component cAMP (fmol).
C STATES(18) is q_PPi in component cAMP (fmol).
C STATES(19) is q_L_B1 in component LRGbinding_B1AR (fmol).
C STATES(20) is q_R_B1 in component LRGbinding_B1AR (fmol).
C STATES(21) is q_Gs in component LRGbinding_B1AR (fmol).
C STATES(22) is q_LR_B1 in component LRGbinding_B1AR (fmol).
C STATES(23) is q_R_B1Gs in component LRGbinding_B1AR (fmol).
C STATES(24) is q_LR_B1Gs in component LRGbinding_B1AR (fmol).
C STATES(25) is q_L_M2 in component LRGbinding_M2 (fmol).
C STATES(26) is q_R_M2 in component LRGbinding_M2 (fmol).
C STATES(27) is q_Gi in component LRGbinding_M2 (fmol).
C STATES(28) is q_LR_M2 in component LRGbinding_M2 (fmol).
C STATES(29) is q_R_M2Gi in component LRGbinding_M2 (fmol).
C STATES(30) is q_LR_M2Gi in component LRGbinding_M2 (fmol).
C STATES(31) is q_R_B1 in component GsProtein (fmol).
C STATES(32) is q_Gs in component GsProtein (fmol).
C STATES(33) is q_R_B1Gs in component GsProtein (fmol).
C STATES(34) is q_LR_B1 in component GsProtein (fmol).
C STATES(35) is q_LR_B1Gs in component GsProtein (fmol).
C STATES(36) is q_a_Gs_GTP in component GsProtein (fmol).
C STATES(37) is q_beta_gamma_Gs in component GsProtein (fmol).
C STATES(38) is q_a_Gs_GDP in component GsProtein (fmol).
C STATES(39) is q_Pi in component GsProtein (fmol).
C STATES(40) is q_R_B1_aby in component GsProtein (fmol).
C STATES(41) is q_R_B1_aby_T in component GsProtein (fmol).
C STATES(42) is q_LR_B1_aby in component GsProtein (fmol).
C STATES(43) is q_LR_B1_aby_T in component GsProtein (fmol).
C STATES(44) is q_GTP in component GsProtein (fmol).
C STATES(45) is q_GDP in component GsProtein (fmol).
C STATES(46) is q_R_M2 in component GiProtein (fmol).
C STATES(47) is q_Gi in component GiProtein (fmol).
C STATES(48) is q_R_M2Gi in component GiProtein (fmol).
C STATES(49) is q_LR_M2 in component GiProtein (fmol).
C STATES(50) is q_LR_M2Gi in component GiProtein (fmol).
C STATES(51) is q_a_Gi_GTP in component GiProtein (fmol).
C STATES(52) is q_beta_gamma_Gi in component GiProtein (fmol).
C STATES(53) is q_a_Gi_GDP in component GiProtein (fmol).
C STATES(54) is q_Pi in component GiProtein (fmol).
C STATES(55) is q_R_M2_aby in component GiProtein (fmol).
C STATES(56) is q_R_M2_aby_T in component GiProtein (fmol).
C STATES(57) is q_LR_M2_aby in component GiProtein (fmol).
C STATES(58) is q_LR_M2_aby_T in component GiProtein (fmol).
C STATES(59) is q_GTP in component GiProtein (fmol).
C STATES(60) is q_GDP in component GiProtein (fmol).
C CONSTS(133) is R in component constants (J_per_K_per_mol).
C CONSTS(134) is T in component constants (kelvin).
C CONSTS(135) is F in component constants (C_per_mol).
C ALGBRC(94) is v1a in component cAMP (fmol_per_sec).
C ALGBRC(97) is v1b in component cAMP (fmol_per_sec).
C ALGBRC(101) is v2a in component cAMP (fmol_per_sec).
C ALGBRC(104) is v2b in component cAMP (fmol_per_sec).
C ALGBRC(107) is v3a in component cAMP (fmol_per_sec).
C ALGBRC(110) is v3b in component cAMP (fmol_per_sec).
C ALGBRC(112) is v4a in component cAMP (fmol_per_sec).
C ALGBRC(119) is v4b in component cAMP (fmol_per_sec).
C ALGBRC(125) is v5 in component cAMP (fmol_per_sec).
C ALGBRC(113) is v6 in component cAMP (fmol_per_sec).
C ALGBRC(120) is v7 in component cAMP (fmol_per_sec).
C ALGBRC(126) is vGiAC in component cAMP (fmol_per_sec).
C ALGBRC(29) is mu_ATP in component cAMP (J_per_mol).
C ALGBRC(33) is mu_AC in component cAMP (J_per_mol).
C ALGBRC(31) is mu_cAMP in component cAMP (J_per_mol).
C ALGBRC(35) is mu_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(77) is mu_FSK in component cAMP (J_per_mol).
C ALGBRC(46) is mu_FSK_AC in component cAMP (J_per_mol).
C ALGBRC(49) is mu_FSK_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(73) is mu_a_Gs_GTP in component cAMP (J_per_mol).
C ALGBRC(37) is mu_a_Gs_GTP_AC in component cAMP (J_per_mol).
C ALGBRC(43) is mu_a_Gs_GTP_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(52) is mu_PDE in component cAMP (J_per_mol).
C ALGBRC(69) is mu_PDEinh in component cAMP (J_per_mol).
C ALGBRC(55) is mu_PDE_cAMP in component cAMP (J_per_mol).
C ALGBRC(65) is mu_IBMX in component cAMP (J_per_mol).
C ALGBRC(58) is mu_five_AMP in component cAMP (J_per_mol).
C ALGBRC(80) is mu_a_Gi_GTP in component cAMP (J_per_mol).
C ALGBRC(83) is mu_ACinh in component cAMP (J_per_mol).
C ALGBRC(89) is mu_PPi in component cAMP (J_per_mol).
C CONSTS(136) is vol in component cAMP (pL).
C ALGBRC(12) is ATP_T in component cAMP (fmol).
C ALGBRC(13) is AC_T in component cAMP (fmol).
C ALGBRC(22) is Gs_T in component cAMP (fmol).
C ALGBRC(16) is cAMP_T in component cAMP (fmol).
C ALGBRC(38) is mu_L_B1 in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(44) is mu_R_B1 in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(47) is mu_Gs in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(50) is mu_LR_B1 in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(53) is mu_R_B1Gs in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(56) is mu_LR_B1Gs in component LRGbinding_B1AR (J_per_mol).
C ALGBRC(59) is vsig1_B1 in component LRGbinding_B1AR (fmol_per_sec).
C ALGBRC(66) is vsig2_B1 in component LRGbinding_B1AR (fmol_per_sec).
C ALGBRC(70) is vsig3_B1 in component LRGbinding_B1AR (fmol_per_sec).
C ALGBRC(74) is vsig4_B1 in component LRGbinding_B1AR (fmol_per_sec).
C CONSTS(137) is vol in component LRGbinding_B1AR (pL).
C ALGBRC(39) is L_T in component LRGbinding_B1AR (fmol).
C ALGBRC(40) is R_T in component LRGbinding_B1AR (fmol).
C ALGBRC(41) is Gs_T in component LRGbinding_B1AR (fmol).
C ALGBRC(60) is mu_L_M2 in component LRGbinding_M2 (J_per_mol).
C ALGBRC(67) is mu_R_M2 in component LRGbinding_M2 (J_per_mol).
C ALGBRC(71) is mu_Gi in component LRGbinding_M2 (J_per_mol).
C ALGBRC(75) is mu_LR_M2 in component LRGbinding_M2 (J_per_mol).
C ALGBRC(78) is mu_R_M2Gi in component LRGbinding_M2 (J_per_mol).
C ALGBRC(81) is mu_LR_M2Gi in component LRGbinding_M2 (J_per_mol).
C ALGBRC(84) is vsig1_M2 in component LRGbinding_M2 (fmol_per_sec).
C ALGBRC(90) is vsig2_M2 in component LRGbinding_M2 (fmol_per_sec).
C ALGBRC(95) is vsig3_M2 in component LRGbinding_M2 (fmol_per_sec).
C ALGBRC(98) is vsig4_M2 in component LRGbinding_M2 (fmol_per_sec).
C CONSTS(138) is vol in component LRGbinding_M2 (pL).
C ALGBRC(61) is L_T in component LRGbinding_M2 (fmol).
C ALGBRC(62) is R_T in component LRGbinding_M2 (fmol).
C ALGBRC(63) is Gi_T in component LRGbinding_M2 (fmol).
C ALGBRC(143) is vDoff1_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(145) is vTon1_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(147) is vAct1_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(149) is vDoff2_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(151) is vTon2_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(153) is vAct2_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(155) is vhyd_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(157) is vreassoc_Gs in component GsProtein (fmol_per_sec).
C ALGBRC(99) is mu_R_B1 in component GsProtein (J_per_mol).
C ALGBRC(102) is mu_Gs in component GsProtein (J_per_mol).
C ALGBRC(105) is mu_R_B1Gs in component GsProtein (J_per_mol).
C ALGBRC(108) is mu_LR_B1 in component GsProtein (J_per_mol).
C ALGBRC(111) is mu_LR_B1Gs in component GsProtein (J_per_mol).
C ALGBRC(114) is mu_a_Gs_GTP in component GsProtein (J_per_mol).
C ALGBRC(121) is mu_beta_gamma_Gs in component GsProtein (J_per_mol).
C ALGBRC(127) is mu_a_Gs_GDP in component GsProtein (J_per_mol).
C ALGBRC(129) is mu_Pi in component GsProtein (J_per_mol).
C ALGBRC(131) is mu_R_B1_aby in component GsProtein (J_per_mol).
C ALGBRC(133) is mu_R_B1_aby_T in component GsProtein (J_per_mol).
C ALGBRC(135) is mu_LR_B1_aby in component GsProtein (J_per_mol).
C ALGBRC(137) is mu_LR_B1_aby_T in component GsProtein (J_per_mol).
C ALGBRC(139) is mu_GTP in component GsProtein (J_per_mol).
C ALGBRC(141) is mu_GDP in component GsProtein (J_per_mol).
C ALGBRC(91) is R_T in component GsProtein (fmol).
C ALGBRC(92) is Gs_T in component GsProtein (fmol).
C ALGBRC(156) is vDoff1_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(158) is vTon1_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(159) is vAct1_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(160) is vDoff2_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(161) is vTon2_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(162) is vAct2_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(163) is vhyd_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(164) is vreassoc_Gi in component GiProtein (fmol_per_sec).
C ALGBRC(122) is mu_R_M2 in component GiProtein (J_per_mol).
C ALGBRC(128) is mu_Gi in component GiProtein (J_per_mol).
C ALGBRC(130) is mu_R_M2Gi in component GiProtein (J_per_mol).
C ALGBRC(132) is mu_LR_M2 in component GiProtein (J_per_mol).
C ALGBRC(134) is mu_LR_M2Gi in component GiProtein (J_per_mol).
C ALGBRC(136) is mu_a_Gi_GTP in component GiProtein (J_per_mol).
C ALGBRC(138) is mu_beta_gamma_Gi in component GiProtein (J_per_mol).
C ALGBRC(140) is mu_a_Gi_GDP in component GiProtein (J_per_mol).
C ALGBRC(142) is mu_Pi in component GiProtein (J_per_mol).
C ALGBRC(144) is mu_R_M2_aby in component GiProtein (J_per_mol).
C ALGBRC(146) is mu_R_M2_aby_T in component GiProtein (J_per_mol).
C ALGBRC(148) is mu_LR_M2_aby in component GiProtein (J_per_mol).
C ALGBRC(150) is mu_LR_M2_aby_T in component GiProtein (J_per_mol).
C ALGBRC(152) is mu_GTP in component GiProtein (J_per_mol).
C ALGBRC(154) is mu_GDP in component GiProtein (J_per_mol).
C ALGBRC(123) is R_T in component GiProtein (fmol).
C ALGBRC(124) is Gi_T in component GiProtein (fmol).
C RATES(1) is d/dt q_ATP in component cAMP (fmol).
C RATES(3) is d/dt q_AC in component cAMP (fmol).
C RATES(4) is d/dt q_AC_ATP in component cAMP (fmol).
C RATES(2) is d/dt q_cAMP in component cAMP (fmol).
C RATES(15) is d/dt q_FSK in component cAMP (fmol).
C RATES(7) is d/dt q_FSK_AC in component cAMP (fmol).
C RATES(8) is d/dt q_FSK_AC_ATP in component cAMP (fmol).
C RATES(14) is d/dt q_a_Gs_GTP in component cAMP (fmol).
C RATES(5) is d/dt q_a_Gs_GTP_AC in component cAMP (fmol).
C RATES(6) is d/dt q_a_Gs_GTP_AC_ATP in component cAMP (fmol).
C RATES(10) is d/dt q_PDE_cAMP in component cAMP (fmol).
C RATES(9) is d/dt q_PDE in component cAMP (fmol).
C RATES(12) is d/dt q_IBMX in component cAMP (fmol).
C RATES(13) is d/dt q_PDEinh in component cAMP (fmol).
C RATES(11) is d/dt q_five_AMP in component cAMP (fmol).
C RATES(16) is d/dt q_a_Gi_GTP in component cAMP (fmol).
C RATES(17) is d/dt q_ACinh in component cAMP (fmol).
C RATES(18) is d/dt q_PPi in component cAMP (fmol).
C RATES(19) is d/dt q_L_B1 in component LRGbinding_B1AR (fmol).
C RATES(20) is d/dt q_R_B1 in component LRGbinding_B1AR (fmol).
C RATES(21) is d/dt q_Gs in component LRGbinding_B1AR (fmol).
C RATES(22) is d/dt q_LR_B1 in component LRGbinding_B1AR (fmol).
C RATES(23) is d/dt q_R_B1Gs in component LRGbinding_B1AR (fmol).
C RATES(24) is d/dt q_LR_B1Gs in component LRGbinding_B1AR (fmol).
C RATES(25) is d/dt q_L_M2 in component LRGbinding_M2 (fmol).
C RATES(26) is d/dt q_R_M2 in component LRGbinding_M2 (fmol).
C RATES(27) is d/dt q_Gi in component LRGbinding_M2 (fmol).
C RATES(28) is d/dt q_LR_M2 in component LRGbinding_M2 (fmol).
C RATES(29) is d/dt q_R_M2Gi in component LRGbinding_M2 (fmol).
C RATES(30) is d/dt q_LR_M2Gi in component LRGbinding_M2 (fmol).
C RATES(31) is d/dt q_R_B1 in component GsProtein (fmol).
C RATES(33) is d/dt q_R_B1Gs in component GsProtein (fmol).
C RATES(32) is d/dt q_Gs in component GsProtein (fmol).
C RATES(34) is d/dt q_LR_B1 in component GsProtein (fmol).
C RATES(35) is d/dt q_LR_B1Gs in component GsProtein (fmol).
C RATES(36) is d/dt q_a_Gs_GTP in component GsProtein (fmol).
C RATES(37) is d/dt q_beta_gamma_Gs in component GsProtein (fmol).
C RATES(38) is d/dt q_a_Gs_GDP in component GsProtein (fmol).
C RATES(39) is d/dt q_Pi in component GsProtein (fmol).
C RATES(40) is d/dt q_R_B1_aby in component GsProtein (fmol).
C RATES(41) is d/dt q_R_B1_aby_T in component GsProtein (fmol).
C RATES(42) is d/dt q_LR_B1_aby in component GsProtein (fmol).
C RATES(43) is d/dt q_LR_B1_aby_T in component GsProtein (fmol).
C RATES(44) is d/dt q_GTP in component GsProtein (fmol).
C RATES(45) is d/dt q_GDP in component GsProtein (fmol).
C RATES(46) is d/dt q_R_M2 in component GiProtein (fmol).
C RATES(48) is d/dt q_R_M2Gi in component GiProtein (fmol).
C RATES(47) is d/dt q_Gi in component GiProtein (fmol).
C RATES(49) is d/dt q_LR_M2 in component GiProtein (fmol).
C RATES(50) is d/dt q_LR_M2Gi in component GiProtein (fmol).
C RATES(51) is d/dt q_a_Gi_GTP in component GiProtein (fmol).
C RATES(52) is d/dt q_beta_gamma_Gi in component GiProtein (fmol).
C RATES(53) is d/dt q_a_Gi_GDP in component GiProtein (fmol).
C RATES(54) is d/dt q_Pi in component GiProtein (fmol).
C RATES(55) is d/dt q_R_M2_aby in component GiProtein (fmol).
C RATES(56) is d/dt q_R_M2_aby_T in component GiProtein (fmol).
C RATES(57) is d/dt q_LR_M2_aby in component GiProtein (fmol).
C RATES(58) is d/dt q_LR_M2_aby_T in component GiProtein (fmol).
C RATES(59) is d/dt q_GTP in component GiProtein (fmol).
C RATES(60) is d/dt q_GDP in component GiProtein (fmol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1.8998e+06
      CONSTS(2) = 0.00039666
      CONSTS(3) = 585859
      CONSTS(4) = 0.158089
      CONSTS(5) = 2.31489e+08
      CONSTS(6) = 2.69173e-17
      CONSTS(7) = 31490.2
      CONSTS(8) = 0.121117
      CONSTS(9) = 815.554
      CONSTS(10) = 8751.63
      CONSTS(11) = 345800
      CONSTS(12) = 2878.6
      CONSTS(13) = 2687.66
      CONSTS(14) = 63.5391
      CONSTS(15) = 14.2293
      CONSTS(16) = 1.13636e+07
      CONSTS(17) = 4069.25
      CONSTS(18) = 51.1466
      CONSTS(19) = 135.504
      CONSTS(20) = 1.78875e+06
      CONSTS(21) = 0.150576
      CONSTS(22) = 0.00442248
      CONSTS(23) = 163.698
      CONSTS(24) = 0.110539
      CONSTS(25) = 0.000540176
      CONSTS(26) = 3.32675
      CONSTS(27) = 0.86602
      CONSTS(28) = 3.20932e+06
      CONSTS(29) = 0.226169
      CONSTS(30) = 0.00556858
      CONSTS(31) = 26.9986
      CONSTS(32) = 0.209142
      CONSTS(33) = 0.00153751
      CONSTS(34) = 0.0445153
      CONSTS(35) = 7.99723
      CONSTS(36) = 1.84852e+06
      CONSTS(37) = 4.34247e-05
      CONSTS(38) = 0.0150226
      CONSTS(39) = 9.52621
      CONSTS(40) = 14.6573
      CONSTS(41) = 3.32164
      CONSTS(42) = 1.563
      CONSTS(43) = 0.0840652
      CONSTS(44) = 0.107997
      CONSTS(45) = 1.78633
      CONSTS(46) = 1.20007
      CONSTS(47) = 0.0150226
      CONSTS(48) = 0.0193351
      CONSTS(49) = 35.6442
      CONSTS(50) = 0.0253404
      CONSTS(51) = 5.83022e-06
      CONSTS(52) = 0.00308164
      CONSTS(53) = 10.0986
      CONSTS(54) = 8.40298e-05
      CONSTS(55) = 0.0361965
      CONSTS(56) = 0.00257076
      CONSTS(57) = 148.667
      CONSTS(58) = 0.310121
      CONSTS(59) = 293.639
      CONSTS(60) = 163.157
      CONSTS(61) = 0.079362
      CONSTS(62) = 0.00692368
      CONSTS(63) = 35.9123
      CONSTS(64) = 0.153698
      CONSTS(65) = 178.993
      CONSTS(66) = 105.77
      CONSTS(67) = 1.51453e-06
      CONSTS(68) = 0.198513
      CONSTS(69) = 5.20902e-06
      CONSTS(70) = 0.0798558
      CONSTS(71) = 0.00432163
      CONSTS(72) = 0.266681
      CONSTS(73) = 0.0867409
      CONSTS(74) = 3639.47
      CONSTS(75) = 4.62054e-05
      CONSTS(76) = 3.06222e-05
      CONSTS(77) = 0.0191425
      CONSTS(78) = 0.0580669
      CONSTS(79) = 0.0037486
      CONSTS(80) = 0.114918
      CONSTS(81) = 0.0248464
      CONSTS(82) = 34.4
      CONSTS(83) = 500
      CONSTS(84) = 3.5e-4
      CONSTS(85) = 3e-4
      CONSTS(86) = 0.25e-4
      CONSTS(87) = 1.8e-4
      CONSTS(88) = 1e1
      CONSTS(89) = 1e-5
      CONSTS(90) = 190
      CONSTS(91) = 1.889E-03
      CONSTS(92) = 1e-18
      CONSTS(93) = 1e-18
      CONSTS(94) = 3.800E-05
      CONSTS(95) = 1e-18
      CONSTS(96) = 1e-18
      CONSTS(97) = 1e-18
      CONSTS(98) = 1e-18
      CONSTS(99) = 1e-18
      CONSTS(100) = 1.482E-03
      CONSTS(101) = 1e-18
      CONSTS(102) = 1e-18
      CONSTS(103) = 3.80E-02
      CONSTS(104) = 1e-18
      CONSTS(105) = 1e-18
      CONSTS(106) = 1e-18
      CONSTS(107) = 1e-18
      CONSTS(108) = 0.0004579000
      CONSTS(109) = 0.1455400000
      CONSTS(110) = 1e-18
      CONSTS(111) = 1e-18
      CONSTS(112) = 1e-18
      CONSTS(113) = 0.00072
      CONSTS(114) = 0.00836
      CONSTS(115) = 1e-18
      CONSTS(116) = 1e-18
      CONSTS(117) = 1e-18
      CONSTS(118) = 1e-18
      CONSTS(119) = 1e-18
      CONSTS(120) = 570
      CONSTS(121) = 1e-18
      CONSTS(122) = 1e-18
      CONSTS(123) = 1e-18
      CONSTS(124) = 1e-18
      CONSTS(125) = 2.2
      CONSTS(126) = 1.1
      CONSTS(127) = 1e-18
      CONSTS(128) = 1e-18
      CONSTS(129) = 1e-18
      CONSTS(130) = 1e-18
      CONSTS(131) = 1e-18
      CONSTS(132) = 1e-18
      STATES(1) = 1e-16
      STATES(2) = 1e-16
      STATES(3) = 1e-16
      STATES(4) = 1e-16
      STATES(5) = 1e-16
      STATES(6) = 1e-16
      STATES(7) = 1e-16
      STATES(8) = 1e-16
      STATES(9) = 1e-16
      STATES(10) = 1e-16
      STATES(11) = 1e-16
      STATES(12) = 1e-16
      STATES(13) = 1e-16
      STATES(14) = 1e-16
      STATES(15) = 1e-16
      STATES(16) = 1e-16
      STATES(17) = 1e-16
      STATES(18) = 1e-16
      STATES(19) = 1e-18
      STATES(20) = 1e-18
      STATES(21) = 1e-18
      STATES(22) = 1e-18
      STATES(23) = 1e-18
      STATES(24) = 1e-18
      STATES(25) = 1e-16
      STATES(26) = 1e-16
      STATES(27) = 1e-16
      STATES(28) = 1e-16
      STATES(29) = 1e-16
      STATES(30) = 1e-16
      STATES(31) = 1e-16
      STATES(32) = 1e-16
      STATES(33) = 1e-16
      STATES(34) = 1e-16
      STATES(35) = 1e-16
      STATES(36) = 1e-16
      STATES(37) = 1e-16
      STATES(38) = 1e-16
      STATES(39) = 1e-16
      STATES(40) = 1e-16
      STATES(41) = 1e-16
      STATES(42) = 1e-16
      STATES(43) = 1e-16
      STATES(44) = 1e-16
      STATES(45) = 1e-16
      STATES(46) = 1e-16
      STATES(47) = 1e-16
      STATES(48) = 1e-16
      STATES(49) = 1e-16
      STATES(50) = 1e-16
      STATES(51) = 1e-16
      STATES(52) = 1e-16
      STATES(53) = 1e-16
      STATES(54) = 1e-16
      STATES(55) = 1e-16
      STATES(56) = 1e-16
      STATES(57) = 1e-16
      STATES(58) = 1e-16
      STATES(59) = 1e-16
      STATES(60) = 1e-16
      CONSTS(133) = 8.31
      CONSTS(134) = 310
      CONSTS(135) = 96485
      CONSTS(136) = 38.0
      CONSTS(137) = 34.4
      CONSTS(138) = 34.4
      CONSTS(139) = CONSTS(88)/CONSTS(87)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(27) = CONSTS(108)+STATES(20)+STATES(31)
      ALGBRC(44) =  CONSTS(133)*CONSTS(134)*log( CONSTS(56)*ALGBRC(27))
      ALGBRC(28) = CONSTS(109)+STATES(21)+STATES(32)
      ALGBRC(47) =  CONSTS(133)*CONSTS(134)*log( CONSTS(57)*ALGBRC(28))
      ALGBRC(32) = CONSTS(111)+STATES(23)+STATES(33)
      ALGBRC(53) =  CONSTS(133)*CONSTS(134)*log( CONSTS(59)*ALGBRC(32))
      ALGBRC(59) =  CONSTS(13)*EXP((ALGBRC(44)+ALGBRC(47))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(53)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(1) = TERNRY(VOI.LT.CONSTS(84).AND.VOI.GT.CONSTS(84) - CONSTS(87), CONSTS(89)+ CONSTS(139)*((VOI - CONSTS(84))+CONSTS(87)), TERNRY(VOI.GE.CONSTS(84).AND.VOI.LT.CONSTS(84)+CONSTS(86), CONSTS(88)+CONSTS(89), TERNRY(VOI.LT.CONSTS(84)+CONSTS(87)+CONSTS(86).AND.VOI.GE.CONSTS(84)+CONSTS(86), CONSTS(89)+ - CONSTS(139)*(((VOI - CONSTS(84)) - CONSTS(87)) - CONSTS(86)), CONSTS(89))
      ALGBRC(25) = ALGBRC(1)+STATES(19)
      ALGBRC(38) =  CONSTS(133)*CONSTS(134)*log( CONSTS(55)*ALGBRC(25))
      ALGBRC(34) = CONSTS(112)+STATES(24)+STATES(35)
      ALGBRC(56) =  CONSTS(133)*CONSTS(134)*log( CONSTS(60)*ALGBRC(34))
      ALGBRC(66) =  CONSTS(14)*EXP((ALGBRC(53)+ALGBRC(38))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(56)/( CONSTS(133)*CONSTS(134)))
      RATES(23) = ALGBRC(59) - ALGBRC(66)
      ALGBRC(30) = CONSTS(110)+STATES(22)+STATES(34)
      ALGBRC(50) =  CONSTS(133)*CONSTS(134)*log( CONSTS(58)*ALGBRC(30))
      ALGBRC(70) =  CONSTS(15)*EXP((ALGBRC(50)+ALGBRC(47))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(56)/( CONSTS(133)*CONSTS(134)))
      RATES(21) = - ALGBRC(59) - ALGBRC(70)
      RATES(24) = ALGBRC(66)+ALGBRC(70)
      ALGBRC(74) =  CONSTS(16)*EXP((ALGBRC(44)+ALGBRC(38))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(50)/( CONSTS(133)*CONSTS(134)))
      RATES(19) = - ALGBRC(66) - ALGBRC(74)
      RATES(20) = - ALGBRC(59) - ALGBRC(74)
      RATES(22) = - ALGBRC(70)+ALGBRC(74)
      ALGBRC(42) = CONSTS(113)+STATES(26)+STATES(46)
      ALGBRC(67) =  CONSTS(133)*CONSTS(134)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(45) = CONSTS(114)+STATES(27)+STATES(47)
      ALGBRC(71) =  CONSTS(133)*CONSTS(134)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(51) = CONSTS(116)+STATES(29)+STATES(48)
      ALGBRC(78) =  CONSTS(133)*CONSTS(134)*log( CONSTS(65)*ALGBRC(51))
      ALGBRC(84) =  CONSTS(17)*EXP((ALGBRC(67)+ALGBRC(71))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(78)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(85).AND.VOI.GT.CONSTS(85) - CONSTS(87), CONSTS(89)+ CONSTS(139)*((VOI - CONSTS(85))+CONSTS(87)), TERNRY(VOI.GE.CONSTS(85).AND.VOI.LT.CONSTS(85)+CONSTS(86), CONSTS(88)+CONSTS(89), TERNRY(VOI.LT.CONSTS(85)+CONSTS(87)+CONSTS(86).AND.VOI.GE.CONSTS(85)+CONSTS(86), CONSTS(89)+ - CONSTS(139)*(((VOI - CONSTS(85)) - CONSTS(87)) - CONSTS(86)), CONSTS(89))
      ALGBRC(36) = ALGBRC(2)+STATES(25)
      ALGBRC(60) =  CONSTS(133)*CONSTS(134)*log( CONSTS(61)*ALGBRC(36))
      ALGBRC(54) = CONSTS(117)+STATES(30)+STATES(50)
      ALGBRC(81) =  CONSTS(133)*CONSTS(134)*log( CONSTS(66)*ALGBRC(54))
      ALGBRC(90) =  CONSTS(18)*EXP((ALGBRC(78)+ALGBRC(60))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(81)/( CONSTS(133)*CONSTS(134)))
      RATES(29) = ALGBRC(84) - ALGBRC(90)
      ALGBRC(48) = CONSTS(115)+STATES(28)+STATES(49)
      ALGBRC(75) =  CONSTS(133)*CONSTS(134)*log( CONSTS(64)*ALGBRC(48))
      ALGBRC(95) =  CONSTS(19)*EXP((ALGBRC(75)+ALGBRC(71))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(81)/( CONSTS(133)*CONSTS(134)))
      RATES(27) = - ALGBRC(84) - ALGBRC(95)
      RATES(30) = ALGBRC(90)+ALGBRC(95)
      ALGBRC(3) = CONSTS(90)+STATES(1)
      ALGBRC(29) =  CONSTS(133)*CONSTS(134)*log( CONSTS(37)*ALGBRC(3))
      ALGBRC(5) = CONSTS(91)+STATES(3)
      ALGBRC(33) =  CONSTS(133)*CONSTS(134)*log( CONSTS(39)*ALGBRC(5))
      ALGBRC(6) = CONSTS(93)+STATES(4)
      ALGBRC(35) =  CONSTS(133)*CONSTS(134)*log( CONSTS(40)*ALGBRC(6))
      ALGBRC(94) =  CONSTS(1)*(EXP((ALGBRC(33)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(35)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(4) = CONSTS(92)+STATES(2)
      ALGBRC(31) =  CONSTS(133)*CONSTS(134)*log( CONSTS(38)*ALGBRC(4))
      ALGBRC(26) = CONSTS(107)+STATES(18)
      ALGBRC(89) =  CONSTS(133)*CONSTS(134)*log( CONSTS(54)*ALGBRC(26))
      ALGBRC(97) =  CONSTS(2)*(EXP(ALGBRC(35)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(33)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      RATES(4) = ALGBRC(94) - ALGBRC(97)
      ALGBRC(98) =  CONSTS(20)*EXP((ALGBRC(67)+ALGBRC(60))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(75)/( CONSTS(133)*CONSTS(134)))
      RATES(25) = - ALGBRC(90) - ALGBRC(98)
      RATES(26) = - ALGBRC(84) - ALGBRC(98)
      RATES(28) = - ALGBRC(95)+ALGBRC(98)
      ALGBRC(7) = CONSTS(98)+STATES(5)
      ALGBRC(37) =  CONSTS(133)*CONSTS(134)*log( CONSTS(41)*ALGBRC(7))
      ALGBRC(8) = CONSTS(99)+STATES(6)
      ALGBRC(43) =  CONSTS(133)*CONSTS(134)*log( CONSTS(42)*ALGBRC(8))
      ALGBRC(101) =  CONSTS(3)*(EXP((ALGBRC(37)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(43)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(104) =  CONSTS(4)*(EXP(ALGBRC(43)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(37)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      RATES(6) = ALGBRC(101) - ALGBRC(104)
      ALGBRC(9) = CONSTS(95)+STATES(7)
      ALGBRC(46) =  CONSTS(133)*CONSTS(134)*log( CONSTS(43)*ALGBRC(9))
      ALGBRC(10) = CONSTS(96)+STATES(8)
      ALGBRC(49) =  CONSTS(133)*CONSTS(134)*log( CONSTS(44)*ALGBRC(10))
      ALGBRC(107) =  CONSTS(5)*(EXP((ALGBRC(46)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(49)/( CONSTS(133)*CONSTS(134))))
      RATES(1) = (- ALGBRC(94) - ALGBRC(107)) - ALGBRC(101)
      ALGBRC(110) =  CONSTS(6)*(EXP(ALGBRC(49)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(46)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      RATES(8) = ALGBRC(107) - ALGBRC(110)
      ALGBRC(11) = CONSTS(100)+STATES(9)
      ALGBRC(52) =  CONSTS(133)*CONSTS(134)*log( CONSTS(45)*ALGBRC(11))
      ALGBRC(14) = CONSTS(102)+STATES(10)
      ALGBRC(55) =  CONSTS(133)*CONSTS(134)*log( CONSTS(46)*ALGBRC(14))
      ALGBRC(112) =  CONSTS(7)*(EXP((ALGBRC(52)+ALGBRC(31))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(55)/( CONSTS(133)*CONSTS(134))))
      RATES(2) = (ALGBRC(97)+ALGBRC(110)+ALGBRC(104)) - ALGBRC(112)
      ALGBRC(20) = CONSTS(97)+STATES(14)+STATES(36)
      ALGBRC(73) =  CONSTS(133)*CONSTS(134)*log( CONSTS(50)*ALGBRC(20))
      ALGBRC(113) =  CONSTS(10)*(EXP((ALGBRC(33)+ALGBRC(73))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(37)/( CONSTS(133)*CONSTS(134))))
      RATES(14) = - ALGBRC(113)
      RATES(5) = (ALGBRC(113) - ALGBRC(101))+ALGBRC(104)
      ALGBRC(21) = CONSTS(94)+STATES(15)
      ALGBRC(77) =  CONSTS(133)*CONSTS(134)*log( CONSTS(51)*ALGBRC(21))
      ALGBRC(120) =  CONSTS(11)*(EXP((ALGBRC(77)+ALGBRC(33))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(46)/( CONSTS(133)*CONSTS(134))))
      RATES(15) = - ALGBRC(120)
      RATES(7) = (ALGBRC(120)+ALGBRC(110)) - ALGBRC(107)
      ALGBRC(15) = CONSTS(104)+STATES(11)
      ALGBRC(58) =  CONSTS(133)*CONSTS(134)*log( CONSTS(47)*ALGBRC(15))
      ALGBRC(119) =  CONSTS(8)*(EXP(ALGBRC(55)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(52)+ALGBRC(58))/( CONSTS(133)*CONSTS(134))))
      RATES(10) = ALGBRC(112) - ALGBRC(119)
      RATES(11) = ALGBRC(119)
      ALGBRC(23) = CONSTS(105)+STATES(16)+STATES(51)
      ALGBRC(80) =  CONSTS(133)*CONSTS(134)*log( CONSTS(52)*ALGBRC(23))
      ALGBRC(24) = CONSTS(106)+STATES(17)
      ALGBRC(83) =  CONSTS(133)*CONSTS(134)*log( CONSTS(53)*ALGBRC(24))
      ALGBRC(126) =  CONSTS(12)*(EXP((ALGBRC(33)+ALGBRC(80))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(83)/( CONSTS(133)*CONSTS(134))))
      RATES(3) = (((ALGBRC(97) - ALGBRC(94)) - ALGBRC(113)) - ALGBRC(120)) - ALGBRC(126)
      ALGBRC(19) = CONSTS(101)+STATES(13)
      ALGBRC(69) =  CONSTS(133)*CONSTS(134)*log( CONSTS(49)*ALGBRC(19))
      ALGBRC(18) = CONSTS(103)+STATES(12)
      ALGBRC(65) =  CONSTS(133)*CONSTS(134)*log( CONSTS(48)*ALGBRC(18))
      ALGBRC(125) =  CONSTS(9)*(EXP((ALGBRC(52)+ALGBRC(65))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(69)/( CONSTS(133)*CONSTS(134))))
      RATES(9) = (ALGBRC(119) - ALGBRC(112)) - ALGBRC(125)
      RATES(12) = - ALGBRC(125)
      RATES(13) = ALGBRC(125)
      RATES(16) = - ALGBRC(126)
      RATES(17) = ALGBRC(126)
      RATES(18) = ALGBRC(126)
      ALGBRC(105) =  CONSTS(133)*CONSTS(134)*log( CONSTS(59)*ALGBRC(32))
      ALGBRC(72) = CONSTS(121)+STATES(40)
      ALGBRC(131) =  CONSTS(133)*CONSTS(134)*log( CONSTS(70)*ALGBRC(72))
      ALGBRC(93) = CONSTS(126)+STATES(45)+STATES(60)
      ALGBRC(141) =  CONSTS(133)*CONSTS(134)*log( CONSTS(75)*ALGBRC(93))
      ALGBRC(143) =  CONSTS(21)*(EXP(ALGBRC(105)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(131)+ALGBRC(141))/( CONSTS(133)*CONSTS(134))))
      RATES(33) = - ALGBRC(143)
      ALGBRC(76) = CONSTS(122)+STATES(41)
      ALGBRC(133) =  CONSTS(133)*CONSTS(134)*log( CONSTS(71)*ALGBRC(76))
      ALGBRC(88) = CONSTS(125)+STATES(44)+STATES(59)
      ALGBRC(139) =  CONSTS(133)*CONSTS(134)*log( CONSTS(74)*ALGBRC(88))
      ALGBRC(145) =  CONSTS(22)*(EXP((ALGBRC(131)+ALGBRC(139))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(133)/( CONSTS(133)*CONSTS(134))))
      RATES(40) = ALGBRC(143) - ALGBRC(145)
      ALGBRC(99) =  CONSTS(133)*CONSTS(134)*log( CONSTS(56)*ALGBRC(27))
      ALGBRC(114) =  CONSTS(133)*CONSTS(134)*log( CONSTS(50)*ALGBRC(20))
      ALGBRC(57) = CONSTS(118)+STATES(37)
      ALGBRC(121) =  CONSTS(133)*CONSTS(134)*log( CONSTS(67)*ALGBRC(57))
      ALGBRC(147) =  CONSTS(23)*(EXP(ALGBRC(133)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(114)+ALGBRC(121)+ALGBRC(99))/( CONSTS(133)*CONSTS(134))))
      RATES(31) = ALGBRC(147)
      RATES(41) = ALGBRC(145) - ALGBRC(147)
      ALGBRC(111) =  CONSTS(133)*CONSTS(134)*log( CONSTS(60)*ALGBRC(34))
      ALGBRC(79) = CONSTS(123)+STATES(42)
      ALGBRC(135) =  CONSTS(133)*CONSTS(134)*log( CONSTS(72)*ALGBRC(79))
      ALGBRC(149) =  CONSTS(24)*(EXP(ALGBRC(111)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(135)+ALGBRC(141))/( CONSTS(133)*CONSTS(134))))
      RATES(35) = - ALGBRC(149)
      RATES(45) = ALGBRC(143)+ALGBRC(149)
      ALGBRC(82) = CONSTS(124)+STATES(43)
      ALGBRC(137) =  CONSTS(133)*CONSTS(134)*log( CONSTS(73)*ALGBRC(82))
      ALGBRC(151) =  CONSTS(25)*(EXP((ALGBRC(135)+ALGBRC(139))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(137)/( CONSTS(133)*CONSTS(134))))
      RATES(42) = ALGBRC(149) - ALGBRC(151)
      RATES(44) = - ALGBRC(145) - ALGBRC(151)
      ALGBRC(108) =  CONSTS(133)*CONSTS(134)*log( CONSTS(58)*ALGBRC(30))
      ALGBRC(153) =  CONSTS(26)*(EXP(ALGBRC(137)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(114)+ALGBRC(121)+ALGBRC(108))/( CONSTS(133)*CONSTS(134))))
      RATES(34) = ALGBRC(153)
      RATES(43) = ALGBRC(151) - ALGBRC(153)
      ALGBRC(64) = CONSTS(119)+STATES(38)
      ALGBRC(127) =  CONSTS(133)*CONSTS(134)*log( CONSTS(68)*ALGBRC(64))
      ALGBRC(68) = CONSTS(120)+STATES(39)+STATES(54)
      ALGBRC(129) =  CONSTS(133)*CONSTS(134)*log( CONSTS(69)*ALGBRC(68))
      ALGBRC(155) =  CONSTS(27)*(EXP(ALGBRC(114)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(127)+ALGBRC(129))/( CONSTS(133)*CONSTS(134))))
      RATES(36) = (ALGBRC(147)+ALGBRC(153)) - ALGBRC(155)
      RATES(39) = ALGBRC(155)
      ALGBRC(130) =  CONSTS(133)*CONSTS(134)*log( CONSTS(65)*ALGBRC(51))
      ALGBRC(103) = CONSTS(129)+STATES(55)
      ALGBRC(144) =  CONSTS(133)*CONSTS(134)*log( CONSTS(78)*ALGBRC(103))
      ALGBRC(154) =  CONSTS(133)*CONSTS(134)*log( CONSTS(75)*ALGBRC(93))
      ALGBRC(156) =  CONSTS(29)*(EXP(ALGBRC(130)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(144)+ALGBRC(154))/( CONSTS(133)*CONSTS(134))))
      RATES(48) = - ALGBRC(156)
      ALGBRC(102) =  CONSTS(133)*CONSTS(134)*log( CONSTS(57)*ALGBRC(28))
      ALGBRC(157) =  CONSTS(28)*(EXP((ALGBRC(127)+ALGBRC(121))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(102)/( CONSTS(133)*CONSTS(134))))
      RATES(32) = ALGBRC(157)
      RATES(37) = (ALGBRC(147)+ALGBRC(153)) - ALGBRC(157)
      RATES(38) = ALGBRC(155) - ALGBRC(157)
      ALGBRC(106) = CONSTS(130)+STATES(56)
      ALGBRC(146) =  CONSTS(133)*CONSTS(134)*log( CONSTS(79)*ALGBRC(106))
      ALGBRC(152) =  CONSTS(133)*CONSTS(134)*log( CONSTS(74)*ALGBRC(88))
      ALGBRC(158) =  CONSTS(30)*(EXP((ALGBRC(144)+ALGBRC(152))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(146)/( CONSTS(133)*CONSTS(134))))
      RATES(55) = ALGBRC(156) - ALGBRC(158)
      ALGBRC(122) =  CONSTS(133)*CONSTS(134)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(136) =  CONSTS(133)*CONSTS(134)*log( CONSTS(52)*ALGBRC(23))
      ALGBRC(96) = CONSTS(127)+STATES(52)
      ALGBRC(138) =  CONSTS(133)*CONSTS(134)*log( CONSTS(76)*ALGBRC(96))
      ALGBRC(159) =  CONSTS(31)*(EXP(ALGBRC(146)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(136)+ALGBRC(138)+ALGBRC(122))/( CONSTS(133)*CONSTS(134))))
      RATES(46) = ALGBRC(159)
      RATES(56) = ALGBRC(158) - ALGBRC(159)
      ALGBRC(134) =  CONSTS(133)*CONSTS(134)*log( CONSTS(66)*ALGBRC(54))
      ALGBRC(109) = CONSTS(131)+STATES(57)
      ALGBRC(148) =  CONSTS(133)*CONSTS(134)*log( CONSTS(80)*ALGBRC(109))
      ALGBRC(160) =  CONSTS(32)*(EXP(ALGBRC(134)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(148)+ALGBRC(154))/( CONSTS(133)*CONSTS(134))))
      RATES(50) = - ALGBRC(160)
      RATES(60) = ALGBRC(156)+ALGBRC(160)
      ALGBRC(115) = CONSTS(132)+STATES(58)
      ALGBRC(150) =  CONSTS(133)*CONSTS(134)*log( CONSTS(81)*ALGBRC(115))
      ALGBRC(161) =  CONSTS(33)*(EXP((ALGBRC(148)+ALGBRC(152))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(150)/( CONSTS(133)*CONSTS(134))))
      RATES(57) = ALGBRC(160) - ALGBRC(161)
      RATES(59) = - ALGBRC(158) - ALGBRC(161)
      ALGBRC(132) =  CONSTS(133)*CONSTS(134)*log( CONSTS(64)*ALGBRC(48))
      ALGBRC(162) =  CONSTS(34)*(EXP(ALGBRC(150)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(136)+ALGBRC(138)+ALGBRC(132))/( CONSTS(133)*CONSTS(134))))
      RATES(49) = ALGBRC(162)
      RATES(58) = ALGBRC(161) - ALGBRC(162)
      ALGBRC(100) = CONSTS(128)+STATES(53)
      ALGBRC(140) =  CONSTS(133)*CONSTS(134)*log( CONSTS(77)*ALGBRC(100))
      ALGBRC(142) =  CONSTS(133)*CONSTS(134)*log( CONSTS(69)*ALGBRC(68))
      ALGBRC(163) =  CONSTS(35)*(EXP(ALGBRC(136)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(140)+ALGBRC(142))/( CONSTS(133)*CONSTS(134))))
      RATES(51) = (ALGBRC(159)+ALGBRC(162)) - ALGBRC(163)
      RATES(54) = ALGBRC(163)
      ALGBRC(128) =  CONSTS(133)*CONSTS(134)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(164) =  CONSTS(36)*(EXP((ALGBRC(140)+ALGBRC(138))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(128)/( CONSTS(133)*CONSTS(134))))
      RATES(47) = ALGBRC(164)
      RATES(52) = (ALGBRC(159)+ALGBRC(162)) - ALGBRC(164)
      RATES(53) = ALGBRC(163) - ALGBRC(164)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(27) = CONSTS(108)+STATES(20)+STATES(31)
      ALGBRC(44) =  CONSTS(133)*CONSTS(134)*log( CONSTS(56)*ALGBRC(27))
      ALGBRC(28) = CONSTS(109)+STATES(21)+STATES(32)
      ALGBRC(47) =  CONSTS(133)*CONSTS(134)*log( CONSTS(57)*ALGBRC(28))
      ALGBRC(32) = CONSTS(111)+STATES(23)+STATES(33)
      ALGBRC(53) =  CONSTS(133)*CONSTS(134)*log( CONSTS(59)*ALGBRC(32))
      ALGBRC(59) =  CONSTS(13)*EXP((ALGBRC(44)+ALGBRC(47))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(53)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(1) = TERNRY(VOI.LT.CONSTS(84).AND.VOI.GT.CONSTS(84) - CONSTS(87), CONSTS(89)+ CONSTS(139)*((VOI - CONSTS(84))+CONSTS(87)), TERNRY(VOI.GE.CONSTS(84).AND.VOI.LT.CONSTS(84)+CONSTS(86), CONSTS(88)+CONSTS(89), TERNRY(VOI.LT.CONSTS(84)+CONSTS(87)+CONSTS(86).AND.VOI.GE.CONSTS(84)+CONSTS(86), CONSTS(89)+ - CONSTS(139)*(((VOI - CONSTS(84)) - CONSTS(87)) - CONSTS(86)), CONSTS(89))
      ALGBRC(25) = ALGBRC(1)+STATES(19)
      ALGBRC(38) =  CONSTS(133)*CONSTS(134)*log( CONSTS(55)*ALGBRC(25))
      ALGBRC(34) = CONSTS(112)+STATES(24)+STATES(35)
      ALGBRC(56) =  CONSTS(133)*CONSTS(134)*log( CONSTS(60)*ALGBRC(34))
      ALGBRC(66) =  CONSTS(14)*EXP((ALGBRC(53)+ALGBRC(38))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(56)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(30) = CONSTS(110)+STATES(22)+STATES(34)
      ALGBRC(50) =  CONSTS(133)*CONSTS(134)*log( CONSTS(58)*ALGBRC(30))
      ALGBRC(70) =  CONSTS(15)*EXP((ALGBRC(50)+ALGBRC(47))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(56)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(74) =  CONSTS(16)*EXP((ALGBRC(44)+ALGBRC(38))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(50)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(42) = CONSTS(113)+STATES(26)+STATES(46)
      ALGBRC(67) =  CONSTS(133)*CONSTS(134)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(45) = CONSTS(114)+STATES(27)+STATES(47)
      ALGBRC(71) =  CONSTS(133)*CONSTS(134)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(51) = CONSTS(116)+STATES(29)+STATES(48)
      ALGBRC(78) =  CONSTS(133)*CONSTS(134)*log( CONSTS(65)*ALGBRC(51))
      ALGBRC(84) =  CONSTS(17)*EXP((ALGBRC(67)+ALGBRC(71))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(78)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(85).AND.VOI.GT.CONSTS(85) - CONSTS(87), CONSTS(89)+ CONSTS(139)*((VOI - CONSTS(85))+CONSTS(87)), TERNRY(VOI.GE.CONSTS(85).AND.VOI.LT.CONSTS(85)+CONSTS(86), CONSTS(88)+CONSTS(89), TERNRY(VOI.LT.CONSTS(85)+CONSTS(87)+CONSTS(86).AND.VOI.GE.CONSTS(85)+CONSTS(86), CONSTS(89)+ - CONSTS(139)*(((VOI - CONSTS(85)) - CONSTS(87)) - CONSTS(86)), CONSTS(89))
      ALGBRC(36) = ALGBRC(2)+STATES(25)
      ALGBRC(60) =  CONSTS(133)*CONSTS(134)*log( CONSTS(61)*ALGBRC(36))
      ALGBRC(54) = CONSTS(117)+STATES(30)+STATES(50)
      ALGBRC(81) =  CONSTS(133)*CONSTS(134)*log( CONSTS(66)*ALGBRC(54))
      ALGBRC(90) =  CONSTS(18)*EXP((ALGBRC(78)+ALGBRC(60))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(81)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(48) = CONSTS(115)+STATES(28)+STATES(49)
      ALGBRC(75) =  CONSTS(133)*CONSTS(134)*log( CONSTS(64)*ALGBRC(48))
      ALGBRC(95) =  CONSTS(19)*EXP((ALGBRC(75)+ALGBRC(71))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(81)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(3) = CONSTS(90)+STATES(1)
      ALGBRC(29) =  CONSTS(133)*CONSTS(134)*log( CONSTS(37)*ALGBRC(3))
      ALGBRC(5) = CONSTS(91)+STATES(3)
      ALGBRC(33) =  CONSTS(133)*CONSTS(134)*log( CONSTS(39)*ALGBRC(5))
      ALGBRC(6) = CONSTS(93)+STATES(4)
      ALGBRC(35) =  CONSTS(133)*CONSTS(134)*log( CONSTS(40)*ALGBRC(6))
      ALGBRC(94) =  CONSTS(1)*(EXP((ALGBRC(33)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(35)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(4) = CONSTS(92)+STATES(2)
      ALGBRC(31) =  CONSTS(133)*CONSTS(134)*log( CONSTS(38)*ALGBRC(4))
      ALGBRC(26) = CONSTS(107)+STATES(18)
      ALGBRC(89) =  CONSTS(133)*CONSTS(134)*log( CONSTS(54)*ALGBRC(26))
      ALGBRC(97) =  CONSTS(2)*(EXP(ALGBRC(35)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(33)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(98) =  CONSTS(20)*EXP((ALGBRC(67)+ALGBRC(60))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(75)/( CONSTS(133)*CONSTS(134)))
      ALGBRC(7) = CONSTS(98)+STATES(5)
      ALGBRC(37) =  CONSTS(133)*CONSTS(134)*log( CONSTS(41)*ALGBRC(7))
      ALGBRC(8) = CONSTS(99)+STATES(6)
      ALGBRC(43) =  CONSTS(133)*CONSTS(134)*log( CONSTS(42)*ALGBRC(8))
      ALGBRC(101) =  CONSTS(3)*(EXP((ALGBRC(37)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(43)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(104) =  CONSTS(4)*(EXP(ALGBRC(43)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(37)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(9) = CONSTS(95)+STATES(7)
      ALGBRC(46) =  CONSTS(133)*CONSTS(134)*log( CONSTS(43)*ALGBRC(9))
      ALGBRC(10) = CONSTS(96)+STATES(8)
      ALGBRC(49) =  CONSTS(133)*CONSTS(134)*log( CONSTS(44)*ALGBRC(10))
      ALGBRC(107) =  CONSTS(5)*(EXP((ALGBRC(46)+ALGBRC(29))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(49)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(110) =  CONSTS(6)*(EXP(ALGBRC(49)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(46)+ALGBRC(31)+ALGBRC(89))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(11) = CONSTS(100)+STATES(9)
      ALGBRC(52) =  CONSTS(133)*CONSTS(134)*log( CONSTS(45)*ALGBRC(11))
      ALGBRC(14) = CONSTS(102)+STATES(10)
      ALGBRC(55) =  CONSTS(133)*CONSTS(134)*log( CONSTS(46)*ALGBRC(14))
      ALGBRC(112) =  CONSTS(7)*(EXP((ALGBRC(52)+ALGBRC(31))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(55)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(20) = CONSTS(97)+STATES(14)+STATES(36)
      ALGBRC(73) =  CONSTS(133)*CONSTS(134)*log( CONSTS(50)*ALGBRC(20))
      ALGBRC(113) =  CONSTS(10)*(EXP((ALGBRC(33)+ALGBRC(73))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(37)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(21) = CONSTS(94)+STATES(15)
      ALGBRC(77) =  CONSTS(133)*CONSTS(134)*log( CONSTS(51)*ALGBRC(21))
      ALGBRC(120) =  CONSTS(11)*(EXP((ALGBRC(77)+ALGBRC(33))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(46)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(15) = CONSTS(104)+STATES(11)
      ALGBRC(58) =  CONSTS(133)*CONSTS(134)*log( CONSTS(47)*ALGBRC(15))
      ALGBRC(119) =  CONSTS(8)*(EXP(ALGBRC(55)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(52)+ALGBRC(58))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(23) = CONSTS(105)+STATES(16)+STATES(51)
      ALGBRC(80) =  CONSTS(133)*CONSTS(134)*log( CONSTS(52)*ALGBRC(23))
      ALGBRC(24) = CONSTS(106)+STATES(17)
      ALGBRC(83) =  CONSTS(133)*CONSTS(134)*log( CONSTS(53)*ALGBRC(24))
      ALGBRC(126) =  CONSTS(12)*(EXP((ALGBRC(33)+ALGBRC(80))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(83)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(19) = CONSTS(101)+STATES(13)
      ALGBRC(69) =  CONSTS(133)*CONSTS(134)*log( CONSTS(49)*ALGBRC(19))
      ALGBRC(18) = CONSTS(103)+STATES(12)
      ALGBRC(65) =  CONSTS(133)*CONSTS(134)*log( CONSTS(48)*ALGBRC(18))
      ALGBRC(125) =  CONSTS(9)*(EXP((ALGBRC(52)+ALGBRC(65))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(69)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(105) =  CONSTS(133)*CONSTS(134)*log( CONSTS(59)*ALGBRC(32))
      ALGBRC(72) = CONSTS(121)+STATES(40)
      ALGBRC(131) =  CONSTS(133)*CONSTS(134)*log( CONSTS(70)*ALGBRC(72))
      ALGBRC(93) = CONSTS(126)+STATES(45)+STATES(60)
      ALGBRC(141) =  CONSTS(133)*CONSTS(134)*log( CONSTS(75)*ALGBRC(93))
      ALGBRC(143) =  CONSTS(21)*(EXP(ALGBRC(105)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(131)+ALGBRC(141))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(76) = CONSTS(122)+STATES(41)
      ALGBRC(133) =  CONSTS(133)*CONSTS(134)*log( CONSTS(71)*ALGBRC(76))
      ALGBRC(88) = CONSTS(125)+STATES(44)+STATES(59)
      ALGBRC(139) =  CONSTS(133)*CONSTS(134)*log( CONSTS(74)*ALGBRC(88))
      ALGBRC(145) =  CONSTS(22)*(EXP((ALGBRC(131)+ALGBRC(139))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(133)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(99) =  CONSTS(133)*CONSTS(134)*log( CONSTS(56)*ALGBRC(27))
      ALGBRC(114) =  CONSTS(133)*CONSTS(134)*log( CONSTS(50)*ALGBRC(20))
      ALGBRC(57) = CONSTS(118)+STATES(37)
      ALGBRC(121) =  CONSTS(133)*CONSTS(134)*log( CONSTS(67)*ALGBRC(57))
      ALGBRC(147) =  CONSTS(23)*(EXP(ALGBRC(133)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(114)+ALGBRC(121)+ALGBRC(99))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(111) =  CONSTS(133)*CONSTS(134)*log( CONSTS(60)*ALGBRC(34))
      ALGBRC(79) = CONSTS(123)+STATES(42)
      ALGBRC(135) =  CONSTS(133)*CONSTS(134)*log( CONSTS(72)*ALGBRC(79))
      ALGBRC(149) =  CONSTS(24)*(EXP(ALGBRC(111)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(135)+ALGBRC(141))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(82) = CONSTS(124)+STATES(43)
      ALGBRC(137) =  CONSTS(133)*CONSTS(134)*log( CONSTS(73)*ALGBRC(82))
      ALGBRC(151) =  CONSTS(25)*(EXP((ALGBRC(135)+ALGBRC(139))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(137)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(108) =  CONSTS(133)*CONSTS(134)*log( CONSTS(58)*ALGBRC(30))
      ALGBRC(153) =  CONSTS(26)*(EXP(ALGBRC(137)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(114)+ALGBRC(121)+ALGBRC(108))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(64) = CONSTS(119)+STATES(38)
      ALGBRC(127) =  CONSTS(133)*CONSTS(134)*log( CONSTS(68)*ALGBRC(64))
      ALGBRC(68) = CONSTS(120)+STATES(39)+STATES(54)
      ALGBRC(129) =  CONSTS(133)*CONSTS(134)*log( CONSTS(69)*ALGBRC(68))
      ALGBRC(155) =  CONSTS(27)*(EXP(ALGBRC(114)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(127)+ALGBRC(129))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(130) =  CONSTS(133)*CONSTS(134)*log( CONSTS(65)*ALGBRC(51))
      ALGBRC(103) = CONSTS(129)+STATES(55)
      ALGBRC(144) =  CONSTS(133)*CONSTS(134)*log( CONSTS(78)*ALGBRC(103))
      ALGBRC(154) =  CONSTS(133)*CONSTS(134)*log( CONSTS(75)*ALGBRC(93))
      ALGBRC(156) =  CONSTS(29)*(EXP(ALGBRC(130)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(144)+ALGBRC(154))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(102) =  CONSTS(133)*CONSTS(134)*log( CONSTS(57)*ALGBRC(28))
      ALGBRC(157) =  CONSTS(28)*(EXP((ALGBRC(127)+ALGBRC(121))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(102)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(106) = CONSTS(130)+STATES(56)
      ALGBRC(146) =  CONSTS(133)*CONSTS(134)*log( CONSTS(79)*ALGBRC(106))
      ALGBRC(152) =  CONSTS(133)*CONSTS(134)*log( CONSTS(74)*ALGBRC(88))
      ALGBRC(158) =  CONSTS(30)*(EXP((ALGBRC(144)+ALGBRC(152))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(146)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(122) =  CONSTS(133)*CONSTS(134)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(136) =  CONSTS(133)*CONSTS(134)*log( CONSTS(52)*ALGBRC(23))
      ALGBRC(96) = CONSTS(127)+STATES(52)
      ALGBRC(138) =  CONSTS(133)*CONSTS(134)*log( CONSTS(76)*ALGBRC(96))
      ALGBRC(159) =  CONSTS(31)*(EXP(ALGBRC(146)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(136)+ALGBRC(138)+ALGBRC(122))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(134) =  CONSTS(133)*CONSTS(134)*log( CONSTS(66)*ALGBRC(54))
      ALGBRC(109) = CONSTS(131)+STATES(57)
      ALGBRC(148) =  CONSTS(133)*CONSTS(134)*log( CONSTS(80)*ALGBRC(109))
      ALGBRC(160) =  CONSTS(32)*(EXP(ALGBRC(134)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(148)+ALGBRC(154))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(115) = CONSTS(132)+STATES(58)
      ALGBRC(150) =  CONSTS(133)*CONSTS(134)*log( CONSTS(81)*ALGBRC(115))
      ALGBRC(161) =  CONSTS(33)*(EXP((ALGBRC(148)+ALGBRC(152))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(150)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(132) =  CONSTS(133)*CONSTS(134)*log( CONSTS(64)*ALGBRC(48))
      ALGBRC(162) =  CONSTS(34)*(EXP(ALGBRC(150)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(136)+ALGBRC(138)+ALGBRC(132))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(100) = CONSTS(128)+STATES(53)
      ALGBRC(140) =  CONSTS(133)*CONSTS(134)*log( CONSTS(77)*ALGBRC(100))
      ALGBRC(142) =  CONSTS(133)*CONSTS(134)*log( CONSTS(69)*ALGBRC(68))
      ALGBRC(163) =  CONSTS(35)*(EXP(ALGBRC(136)/( CONSTS(133)*CONSTS(134))) - EXP((ALGBRC(140)+ALGBRC(142))/( CONSTS(133)*CONSTS(134))))
      ALGBRC(128) =  CONSTS(133)*CONSTS(134)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(164) =  CONSTS(36)*(EXP((ALGBRC(140)+ALGBRC(138))/( CONSTS(133)*CONSTS(134))) - EXP(ALGBRC(128)/( CONSTS(133)*CONSTS(134))))
      ALGBRC(12) = ALGBRC(3)+ALGBRC(6)+ALGBRC(10)+ALGBRC(8)
      ALGBRC(13) = ALGBRC(5)+ALGBRC(6)+ALGBRC(9)+ALGBRC(10)+ALGBRC(7)+ALGBRC(8)
      ALGBRC(16) = ALGBRC(4)+ALGBRC(14)+STATES(11)
      ALGBRC(17) = ALGBRC(4)+ALGBRC(14)+ALGBRC(15)+ALGBRC(3)+ALGBRC(6)+ALGBRC(8)+ALGBRC(10)
      ALGBRC(22) = ALGBRC(20)+ALGBRC(7)+ALGBRC(8)
      ALGBRC(39) = ALGBRC(25)+ALGBRC(30)+ALGBRC(34)
      ALGBRC(40) = ALGBRC(27)+ALGBRC(30)+ALGBRC(32)+ALGBRC(34)
      ALGBRC(41) = ALGBRC(28)+ALGBRC(32)+ALGBRC(34)
      ALGBRC(61) = ALGBRC(36)+ALGBRC(48)+ALGBRC(54)
      ALGBRC(62) = ALGBRC(42)+ALGBRC(48)+ALGBRC(51)+ALGBRC(54)
      ALGBRC(63) = ALGBRC(45)+ALGBRC(51)+ALGBRC(54)
      ALGBRC(85) = ALGBRC(25)+ALGBRC(34)+ALGBRC(30)+ALGBRC(79)+ALGBRC(82)
      ALGBRC(86) = ALGBRC(27)+ALGBRC(32)+ALGBRC(30)+ALGBRC(34)+ALGBRC(72)+ALGBRC(76)+ALGBRC(79)+ALGBRC(82)
      ALGBRC(87) = ALGBRC(28)+ALGBRC(32)+ALGBRC(34)+ALGBRC(20)+ALGBRC(64)+ALGBRC(72)+ALGBRC(76)+ALGBRC(79)+ALGBRC(82)
      ALGBRC(91) = ALGBRC(27)+ALGBRC(32)+ALGBRC(30)+ALGBRC(34)+ALGBRC(72)+ALGBRC(76)+ALGBRC(79)+ALGBRC(82)
      ALGBRC(92) = ALGBRC(28)+ALGBRC(32)+ALGBRC(34)+ALGBRC(20)+ALGBRC(64)+ALGBRC(72)+ALGBRC(76)+ALGBRC(79)+ALGBRC(82)
      ALGBRC(116) = ALGBRC(36)+ALGBRC(54)+ALGBRC(48)+ALGBRC(109)+ALGBRC(115)
      ALGBRC(117) = ALGBRC(42)+ALGBRC(51)+ALGBRC(48)+ALGBRC(54)+ALGBRC(103)+ALGBRC(106)+ALGBRC(109)+ALGBRC(115)
      ALGBRC(118) = ALGBRC(45)+ALGBRC(51)+ALGBRC(54)+ALGBRC(23)+ALGBRC(100)+ALGBRC(103)+ALGBRC(106)+ALGBRC(109)+ALGBRC(115)
      ALGBRC(123) = ALGBRC(42)+ALGBRC(51)+ALGBRC(48)+ALGBRC(54)+ALGBRC(103)+ALGBRC(106)+ALGBRC(109)+ALGBRC(115)
      ALGBRC(124) = ALGBRC(45)+ALGBRC(51)+ALGBRC(54)+ALGBRC(23)+ALGBRC(100)+ALGBRC(103)+ALGBRC(106)+ALGBRC(109)+ALGBRC(115)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END