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 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