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 148 entries in the algebraic variable array.
C There are a total of 103 entries in each of the rate and state variable arrays.
C There are a total of 96 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is c1 in component c1 (molar).
C ALGBRC(1) is v1 in component v1 (flux).
C STATES(2) is c2 in component c2 (molar).
C ALGBRC(6) is v6 in component v6 (flux).
C ALGBRC(7) is v13 in component v13 (flux).
C STATES(3) is c3 in component c3 (molar).
C ALGBRC(2) is v2 in component v2 (flux).
C STATES(4) is c4 in component c4 (molar).
C ALGBRC(3) is v3 in component v3 (flux).
C STATES(5) is c5 in component c5 (molar).
C ALGBRC(8) is v7 in component v7 (flux).
C ALGBRC(9) is v8 in component v8 (flux).
C STATES(6) is c6 in component c6 (molar).
C ALGBRC(11) is v10 in component v10 (flux).
C ALGBRC(12) is v60 in component v60 (flux).
C STATES(7) is c7 in component c7 (molar).
C ALGBRC(4) is v4 in component v4 (flux).
C ALGBRC(5) is v5 in component v5 (flux).
C STATES(8) is c8 in component c8 (molar).
C ALGBRC(15) is v12 in component v12 (flux).
C ALGBRC(16) is v14 in component v14 (flux).
C ALGBRC(62) is v62 in component v62 (flux).
C STATES(9) is c9 in component c9 (molar).
C ALGBRC(17) is v15 in component v15 (flux).
C ALGBRC(107) is v107 in component v107 (flux).
C ALGBRC(110) is v110 in component v110 (flux).
C ALGBRC(113) is v113 in component v113 (flux).
C ALGBRC(116) is v116 in component v116 (flux).
C ALGBRC(119) is v119 in component v119 (flux).
C ALGBRC(122) is v122 in component v122 (flux).
C ALGBRC(125) is v125 in component v125 (flux).
C STATES(10) is c10 in component c10 (molar).
C ALGBRC(13) is v11 in component v11 (flux).
C STATES(11) is c11 in component c11 (molar).
C STATES(12) is c12 in component c12 (molar).
C ALGBRC(106) is v106 in component v106 (flux).
C ALGBRC(109) is v109 in component v109 (flux).
C ALGBRC(112) is v112 in component v112 (flux).
C ALGBRC(115) is v115 in component v115 (flux).
C ALGBRC(118) is v118 in component v118 (flux).
C ALGBRC(121) is v121 in component v121 (flux).
C ALGBRC(124) is v124 in component v124 (flux).
C STATES(13) is c13 in component c13 (molar).
C ALGBRC(14) is v61 in component v61 (flux).
C STATES(14) is c14 in component c14 (molar).
C STATES(15) is c15 in component c15 (molar).
C ALGBRC(18) is v16 in component v16 (flux).
C ALGBRC(24) is v22 in component v22 (flux).
C ALGBRC(34) is v32 in component v32 (flux).
C ALGBRC(36) is v34 in component v34 (flux).
C ALGBRC(39) is v37 in component v37 (flux).
C ALGBRC(41) is v39 in component v39 (flux).
C ALGBRC(99) is v102 in component v102 (flux).
C STATES(16) is c16 in component c16 (molar).
C STATES(17) is c17 in component c17 (molar).
C ALGBRC(63) is v63 in component v63 (flux).
C ALGBRC(69) is v69 in component v69 (flux).
C ALGBRC(79) is v79 in component v79 (flux).
C ALGBRC(80) is v80 in component v80 (flux).
C ALGBRC(81) is v81 in component v81 (flux).
C ALGBRC(82) is v82 in component v82 (flux).
C ALGBRC(130) is v132 in component v132 (flux).
C STATES(18) is c18 in component c18 (molar).
C ALGBRC(10) is v9 in component v9 (flux).
C ALGBRC(64) is v64 in component v64 (flux).
C ALGBRC(132) is v133 in component v133 (flux).
C STATES(19) is c19 in component c19 (molar).
C ALGBRC(65) is v65 in component v65 (flux).
C ALGBRC(66) is v66 in component v66 (flux).
C ALGBRC(67) is v67 in component v67 (flux).
C ALGBRC(68) is v68 in component v68 (flux).
C ALGBRC(105) is v105 in component v105 (flux).
C ALGBRC(127) is v127 in component v127 (flux).
C ALGBRC(134) is v134 in component v134 (flux).
C STATES(20) is c20 in component c20 (molar).
C ALGBRC(108) is v108 in component v108 (flux).
C ALGBRC(136) is v135 in component v135 (flux).
C STATES(21) is c21 in component c21 (molar).
C ALGBRC(111) is v111 in component v111 (flux).
C ALGBRC(138) is v136 in component v136 (flux).
C STATES(22) is c22 in component c22 (molar).
C ALGBRC(26) is v24 in component v24 (flux).
C ALGBRC(37) is v35 in component v35 (flux).
C ALGBRC(40) is v38 in component v38 (flux).
C ALGBRC(71) is v71 in component v71 (flux).
C STATES(23) is c23 in component c23 (molar).
C ALGBRC(19) is v17 in component v17 (flux).
C STATES(24) is c24 in component c24 (molar).
C ALGBRC(72) is v72 in component v72 (flux).
C ALGBRC(131) is v130 in component v130 (flux).
C ALGBRC(133) is v131 in component v131 (flux).
C ALGBRC(27) is v25 in component v25 (flux).
C ALGBRC(42) is v40 in component v40 (flux).
C STATES(25) is c25 in component c25 (molar).
C ALGBRC(20) is v18 in component v18 (flux).
C ALGBRC(21) is v19 in component v19 (flux).
C ALGBRC(22) is v20 in component v20 (flux).
C ALGBRC(23) is v21 in component v21 (flux).
C ALGBRC(126) is v126 in component v126 (flux).
C STATES(26) is c26 in component c26 (molar).
C ALGBRC(28) is v26 in component v26 (flux).
C ALGBRC(33) is v31 in component v31 (flux).
C ALGBRC(78) is v78 in component v78 (flux).
C ALGBRC(73) is v73 in component v73 (flux).
C STATES(27) is c27 in component c27 (molar).
C STATES(28) is c28 in component c28 (molar).
C ALGBRC(29) is v27 in component v27 (flux).
C ALGBRC(30) is v28 in component v28 (flux).
C STATES(29) is c29 in component c29 (molar).
C STATES(30) is c30 in component c30 (molar).
C ALGBRC(35) is v33 in component v33 (flux).
C ALGBRC(43) is v41 in component v41 (flux).
C ALGBRC(83) is v83 in component v83 (flux).
C STATES(31) is c31 in component c31 (molar).
C ALGBRC(38) is v36 in component v36 (flux).
C STATES(32) is c32 in component c32 (molar).
C ALGBRC(25) is v23 in component v23 (flux).
C ALGBRC(101) is v103 in component v103 (flux).
C STATES(33) is c33 in component c33 (molar).
C ALGBRC(103) is v104 in component v104 (flux).
C STATES(34) is c34 in component c34 (molar).
C ALGBRC(114) is v114 in component v114 (flux).
C STATES(35) is c35 in component c35 (molar).
C ALGBRC(32) is v30 in component v30 (flux).
C ALGBRC(117) is v117 in component v117 (flux).
C ALGBRC(128) is v128 in component v128 (flux).
C STATES(36) is c36 in component c36 (molar).
C ALGBRC(120) is v120 in component v120 (flux).
C STATES(37) is c37 in component c37 (molar).
C ALGBRC(123) is v123 in component v123 (flux).
C STATES(38) is c38 in component c38 (molar).
C STATES(39) is c39 in component c39 (molar).
C STATES(40) is c40 in component c40 (molar).
C STATES(41) is c41 in component c41 (molar).
C ALGBRC(45) is v43 in component v43 (flux).
C ALGBRC(75) is v75 in component v75 (flux).
C ALGBRC(85) is v85 in component v85 (flux).
C STATES(42) is c42 in component c42 (molar).
C ALGBRC(31) is v29 in component v29 (flux).
C STATES(43) is c43 in component c43 (molar).
C STATES(44) is c44 in component c44 (molar).
C ALGBRC(44) is v42 in component v42 (flux).
C ALGBRC(84) is v84 in component v84 (flux).
C STATES(45) is c45 in component c45 (molar).
C ALGBRC(46) is v44 in component v44 (flux).
C ALGBRC(47) is v45 in component v45 (flux).
C ALGBRC(48) is v46 in component v46 (flux).
C ALGBRC(49) is v47 in component v47 (flux).
C STATES(46) is c46 in component c46 (molar).
C STATES(47) is c47 in component c47 (molar).
C ALGBRC(53) is v51 in component v51 (flux).
C ALGBRC(86) is v86 in component v86 (flux).
C ALGBRC(96) is v93 in component v93 (flux).
C STATES(48) is c48 in component c48 (molar).
C STATES(49) is c49 in component c49 (molar).
C ALGBRC(51) is v49 in component v49 (flux).
C ALGBRC(52) is v50 in component v50 (flux).
C STATES(50) is c50 in component c50 (molar).
C STATES(51) is c51 in component c51 (molar).
C ALGBRC(50) is v48 in component v48 (flux).
C ALGBRC(54) is v52 in component v52 (flux).
C ALGBRC(55) is v53 in component v53 (flux).
C ALGBRC(56) is v54 in component v54 (flux).
C ALGBRC(57) is v55 in component v55 (flux).
C STATES(52) is c52 in component c52 (molar).
C STATES(53) is c53 in component c53 (molar).
C ALGBRC(90) is v90 in component v90 (flux).
C ALGBRC(92) is v91 in component v91 (flux).
C ALGBRC(94) is v92 in component v92 (flux).
C STATES(54) is c54 in component c54 (molar).
C STATES(55) is c55 in component c55 (molar).
C ALGBRC(61) is v59 in component v59 (flux).
C ALGBRC(91) is v94 in component v94 (flux).
C ALGBRC(104) is v101 in component v101 (flux).
C STATES(56) is c56 in component c56 (molar).
C STATES(57) is c57 in component c57 (molar).
C ALGBRC(59) is v57 in component v57 (flux).
C ALGBRC(60) is v58 in component v58 (flux).
C STATES(58) is c58 in component c58 (molar).
C STATES(59) is c59 in component c59 (molar).
C ALGBRC(58) is v56 in component v56 (flux).
C ALGBRC(135) is v143 in component v143 (flux).
C ALGBRC(137) is v144 in component v144 (flux).
C ALGBRC(139) is v145 in component v145 (flux).
C STATES(60) is c60 in component c60 (molar).
C ALGBRC(98) is v98 in component v98 (flux).
C ALGBRC(100) is v99 in component v99 (flux).
C ALGBRC(102) is v100 in component v100 (flux).
C STATES(61) is c61 in component c61 (molar).
C STATES(62) is c62 in component c62 (molar).
C STATES(63) is c63 in component c63 (molar).
C ALGBRC(70) is v70 in component v70 (flux).
C ALGBRC(140) is v137 in component v137 (flux).
C STATES(64) is c64 in component c64 (molar).
C ALGBRC(142) is v138 in component v138 (flux).
C STATES(65) is c65 in component c65 (molar).
C ALGBRC(144) is v139 in component v139 (flux).
C STATES(66) is c66 in component c66 (molar).
C ALGBRC(74) is v74 in component v74 (flux).
C ALGBRC(77) is v77 in component v77 (flux).
C ALGBRC(129) is v129 in component v129 (flux).
C ALGBRC(146) is v140 in component v140 (flux).
C STATES(67) is c67 in component c67 (molar).
C ALGBRC(147) is v141 in component v141 (flux).
C STATES(68) is c68 in component c68 (molar).
C ALGBRC(148) is v142 in component v142 (flux).
C STATES(69) is c69 in component c69 (molar).
C STATES(70) is c70 in component c70 (molar).
C ALGBRC(76) is v76 in component v76 (flux).
C STATES(71) is c71 in component c71 (molar).
C STATES(72) is c72 in component c72 (molar).
C ALGBRC(87) is v87 in component v87 (flux).
C ALGBRC(88) is v88 in component v88 (flux).
C ALGBRC(89) is v89 in component v89 (flux).
C STATES(73) is c73 in component c73 (molar).
C STATES(74) is c74 in component c74 (molar).
C STATES(75) is c75 in component c75 (molar).
C STATES(76) is c76 in component c76 (molar).
C STATES(77) is c77 in component c77 (molar).
C ALGBRC(93) is v95 in component v95 (flux).
C ALGBRC(95) is v96 in component v96 (flux).
C ALGBRC(97) is v97 in component v97 (flux).
C STATES(78) is c78 in component c78 (molar).
C STATES(79) is c79 in component c79 (molar).
C STATES(80) is c80 in component c80 (molar).
C STATES(81) is c81 in component c81 (molar).
C STATES(82) is c82 in component c82 (molar).
C STATES(83) is c83 in component c83 (molar).
C ALGBRC(141) is v146 in component v146 (flux).
C ALGBRC(143) is v147 in component v147 (flux).
C ALGBRC(145) is v148 in component v148 (flux).
C STATES(84) is c84 in component c84 (molar).
C STATES(85) is c85 in component c85 (molar).
C STATES(86) is c86 in component c86 (molar).
C STATES(87) is c87 in component c87 (molar).
C STATES(88) is c88 in component c88 (molar).
C STATES(89) is c89 in component c89 (molar).
C STATES(90) is c90 in component c90 (molar).
C STATES(91) is c91 in component c91 (molar).
C STATES(92) is c92 in component c92 (molar).
C STATES(93) is c93 in component c93 (molar).
C STATES(94) is c94 in component c94 (molar).
C STATES(95) is c95 in component c95 (molar).
C STATES(96) is c96 in component c96 (molar).
C STATES(97) is c97 in component c97 (molar).
C STATES(98) is c98 in component c98 (molar).
C STATES(99) is c99 in component c99 (molar).
C STATES(100) is c100 in component c100 (molar).
C STATES(101) is c101 in component c101 (molar).
C STATES(102) is c102 in component c102 (molar).
C STATES(103) is c103 in component c103 (molar).
C CONSTS(1) is k1 in component reaction_rates (second_order_rate_constant).
C CONSTS(2) is kd1 in component reaction_rates (first_order_rate_constant).
C CONSTS(3) is k2 in component reaction_rates (second_order_rate_constant).
C CONSTS(4) is kd2 in component reaction_rates (first_order_rate_constant).
C CONSTS(5) is k3 in component reaction_rates (first_order_rate_constant).
C CONSTS(6) is kd3 in component reaction_rates (first_order_rate_constant).
C CONSTS(7) is k4 in component reaction_rates (second_order_rate_constant).
C CONSTS(8) is kd4 in component reaction_rates (first_order_rate_constant).
C CONSTS(9) is k5 in component reaction_rates (second_order_rate_constant).
C CONSTS(10) is kd5 in component reaction_rates (first_order_rate_constant).
C CONSTS(11) is k6 in component reaction_rates (first_order_rate_constant).
C CONSTS(12) is kd6 in component reaction_rates (first_order_rate_constant).
C CONSTS(13) is k8 in component reaction_rates (second_order_rate_constant).
C CONSTS(14) is kd8 in component reaction_rates (first_order_rate_constant).
C CONSTS(15) is k10b in component reaction_rates (second_order_rate_constant).
C CONSTS(16) is kd10 in component reaction_rates (first_order_rate_constant).
C CONSTS(17) is k13 in component reaction_rates (first_order_rate_constant).
C CONSTS(18) is kd13 in component reaction_rates (first_order_rate_constant).
C CONSTS(19) is k15 in component reaction_rates (first_order_rate_constant).
C CONSTS(20) is kd15 in component reaction_rates (first_order_rate_constant).
C CONSTS(21) is k16 in component reaction_rates (second_order_rate_constant).
C CONSTS(22) is kd63 in component reaction_rates (first_order_rate_constant).
C CONSTS(23) is k17 in component reaction_rates (second_order_rate_constant).
C CONSTS(24) is kd17 in component reaction_rates (first_order_rate_constant).
C CONSTS(25) is k18 in component reaction_rates (second_order_rate_constant).
C CONSTS(26) is kd18 in component reaction_rates (first_order_rate_constant).
C CONSTS(27) is k19 in component reaction_rates (second_order_rate_constant).
C CONSTS(28) is kd19 in component reaction_rates (first_order_rate_constant).
C CONSTS(29) is k20 in component reaction_rates (second_order_rate_constant).
C CONSTS(30) is kd20 in component reaction_rates (first_order_rate_constant).
C CONSTS(31) is k21 in component reaction_rates (second_order_rate_constant).
C CONSTS(32) is kd21 in component reaction_rates (first_order_rate_constant).
C CONSTS(33) is k22 in component reaction_rates (second_order_rate_constant).
C CONSTS(34) is kd22 in component reaction_rates (first_order_rate_constant).
C CONSTS(35) is k23 in component reaction_rates (first_order_rate_constant).
C CONSTS(36) is kd23 in component reaction_rates (first_order_rate_constant).
C CONSTS(37) is kd24 in component reaction_rates (first_order_rate_constant).
C CONSTS(38) is k25 in component reaction_rates (second_order_rate_constant).
C CONSTS(39) is kd25 in component reaction_rates (first_order_rate_constant).
C CONSTS(40) is k28 in component reaction_rates (second_order_rate_constant).
C CONSTS(41) is kd28 in component reaction_rates (first_order_rate_constant).
C CONSTS(42) is k29 in component reaction_rates (second_order_rate_constant).
C CONSTS(43) is kd29 in component reaction_rates (first_order_rate_constant).
C CONSTS(44) is k32 in component reaction_rates (second_order_rate_constant).
C CONSTS(45) is kd32 in component reaction_rates (first_order_rate_constant).
C CONSTS(46) is k33 in component reaction_rates (second_order_rate_constant).
C CONSTS(47) is kd33 in component reaction_rates (first_order_rate_constant).
C CONSTS(48) is k34 in component reaction_rates (second_order_rate_constant).
C CONSTS(49) is kd34 in component reaction_rates (first_order_rate_constant).
C CONSTS(50) is k35 in component reaction_rates (second_order_rate_constant).
C CONSTS(51) is kd35 in component reaction_rates (first_order_rate_constant).
C CONSTS(52) is k36 in component reaction_rates (first_order_rate_constant).
C CONSTS(53) is kd36 in component reaction_rates (first_order_rate_constant).
C CONSTS(54) is k37 in component reaction_rates (second_order_rate_constant).
C CONSTS(55) is kd37 in component reaction_rates (first_order_rate_constant).
C CONSTS(56) is k40 in component reaction_rates (second_order_rate_constant).
C CONSTS(57) is kd40 in component reaction_rates (first_order_rate_constant).
C CONSTS(58) is k41 in component reaction_rates (second_order_rate_constant).
C CONSTS(59) is kd41 in component reaction_rates (first_order_rate_constant).
C CONSTS(60) is k42 in component reaction_rates (second_order_rate_constant).
C CONSTS(61) is kd42 in component reaction_rates (first_order_rate_constant).
C CONSTS(62) is k43 in component reaction_rates (second_order_rate_constant).
C CONSTS(63) is kd43 in component reaction_rates (first_order_rate_constant).
C CONSTS(64) is k44 in component reaction_rates (second_order_rate_constant).
C CONSTS(65) is kd52 in component reaction_rates (first_order_rate_constant).
C CONSTS(66) is k45 in component reaction_rates (second_order_rate_constant).
C CONSTS(67) is kd45 in component reaction_rates (first_order_rate_constant).
C CONSTS(68) is k47 in component reaction_rates (second_order_rate_constant).
C CONSTS(69) is kd47 in component reaction_rates (first_order_rate_constant).
C CONSTS(70) is k48 in component reaction_rates (second_order_rate_constant).
C CONSTS(71) is kd48 in component reaction_rates (first_order_rate_constant).
C CONSTS(72) is k49 in component reaction_rates (second_order_rate_constant).
C CONSTS(73) is kd49 in component reaction_rates (first_order_rate_constant).
C CONSTS(74) is k50 in component reaction_rates (second_order_rate_constant).
C CONSTS(75) is kd50 in component reaction_rates (first_order_rate_constant).
C CONSTS(76) is k52 in component reaction_rates (second_order_rate_constant).
C CONSTS(77) is kd44 in component reaction_rates (first_order_rate_constant).
C CONSTS(78) is k53 in component reaction_rates (second_order_rate_constant).
C CONSTS(79) is kd53 in component reaction_rates (first_order_rate_constant).
C CONSTS(80) is k55 in component reaction_rates (second_order_rate_constant).
C CONSTS(81) is kd55 in component reaction_rates (first_order_rate_constant).
C CONSTS(82) is k56 in component reaction_rates (second_order_rate_constant).
C CONSTS(83) is kd56 in component reaction_rates (first_order_rate_constant).
C CONSTS(84) is k57 in component reaction_rates (second_order_rate_constant).
C CONSTS(85) is kd57 in component reaction_rates (first_order_rate_constant).
C CONSTS(86) is k58 in component reaction_rates (second_order_rate_constant).
C CONSTS(87) is kd58 in component reaction_rates (first_order_rate_constant).
C CONSTS(88) is k60 in component reaction_rates (first_order_rate_constant).
C CONSTS(89) is kd60 in component reaction_rates (first_order_rate_constant).
C CONSTS(90) is k61 in component reaction_rates (first_order_rate_constant).
C CONSTS(91) is kd61 in component reaction_rates (first_order_rate_constant).
C CONSTS(92) is k126 in component reaction_rates (second_order_rate_constant).
C CONSTS(93) is kd126 in component reaction_rates (first_order_rate_constant).
C CONSTS(94) is k127 in component reaction_rates (second_order_rate_constant).
C CONSTS(95) is kd127 in component reaction_rates (first_order_rate_constant).
C RATES(1) is d/dt c1 in component c1 (molar).
C RATES(2) is d/dt c2 in component c2 (molar).
C RATES(3) is d/dt c3 in component c3 (molar).
C RATES(4) is d/dt c4 in component c4 (molar).
C RATES(5) is d/dt c5 in component c5 (molar).
C RATES(6) is d/dt c6 in component c6 (molar).
C RATES(7) is d/dt c7 in component c7 (molar).
C RATES(8) is d/dt c8 in component c8 (molar).
C RATES(9) is d/dt c9 in component c9 (molar).
C RATES(10) is d/dt c10 in component c10 (molar).
C RATES(11) is d/dt c11 in component c11 (molar).
C RATES(12) is d/dt c12 in component c12 (molar).
C RATES(13) is d/dt c13 in component c13 (molar).
C RATES(14) is d/dt c14 in component c14 (molar).
C RATES(15) is d/dt c15 in component c15 (molar).
C RATES(16) is d/dt c16 in component c16 (molar).
C RATES(17) is d/dt c17 in component c17 (molar).
C RATES(18) is d/dt c18 in component c18 (molar).
C RATES(19) is d/dt c19 in component c19 (molar).
C RATES(20) is d/dt c20 in component c20 (molar).
C RATES(21) is d/dt c21 in component c21 (molar).
C RATES(22) is d/dt c22 in component c22 (molar).
C RATES(23) is d/dt c23 in component c23 (molar).
C RATES(24) is d/dt c24 in component c24 (molar).
C RATES(25) is d/dt c25 in component c25 (molar).
C RATES(26) is d/dt c26 in component c26 (molar).
C RATES(27) is d/dt c27 in component c27 (molar).
C RATES(28) is d/dt c28 in component c28 (molar).
C RATES(29) is d/dt c29 in component c29 (molar).
C RATES(30) is d/dt c30 in component c30 (molar).
C RATES(31) is d/dt c31 in component c31 (molar).
C RATES(32) is d/dt c32 in component c32 (molar).
C RATES(33) is d/dt c33 in component c33 (molar).
C RATES(34) is d/dt c34 in component c34 (molar).
C RATES(35) is d/dt c35 in component c35 (molar).
C RATES(36) is d/dt c36 in component c36 (molar).
C RATES(37) is d/dt c37 in component c37 (molar).
C RATES(38) is d/dt c38 in component c38 (molar).
C RATES(39) is d/dt c39 in component c39 (molar).
C RATES(40) is d/dt c40 in component c40 (molar).
C RATES(41) is d/dt c41 in component c41 (molar).
C RATES(42) is d/dt c42 in component c42 (molar).
C RATES(43) is d/dt c43 in component c43 (molar).
C RATES(44) is d/dt c44 in component c44 (molar).
C RATES(45) is d/dt c45 in component c45 (molar).
C RATES(46) is d/dt c46 in component c46 (molar).
C RATES(47) is d/dt c47 in component c47 (molar).
C RATES(48) is d/dt c48 in component c48 (molar).
C RATES(49) is d/dt c49 in component c49 (molar).
C RATES(50) is d/dt c50 in component c50 (molar).
C RATES(51) is d/dt c51 in component c51 (molar).
C RATES(52) is d/dt c52 in component c52 (molar).
C RATES(53) is d/dt c53 in component c53 (molar).
C RATES(54) is d/dt c54 in component c54 (molar).
C RATES(55) is d/dt c55 in component c55 (molar).
C RATES(56) is d/dt c56 in component c56 (molar).
C RATES(57) is d/dt c57 in component c57 (molar).
C RATES(58) is d/dt c58 in component c58 (molar).
C RATES(59) is d/dt c59 in component c59 (molar).
C RATES(60) is d/dt c60 in component c60 (molar).
C RATES(61) is d/dt c61 in component c61 (molar).
C RATES(62) is d/dt c62 in component c62 (molar).
C RATES(63) is d/dt c63 in component c63 (molar).
C RATES(64) is d/dt c64 in component c64 (molar).
C RATES(65) is d/dt c65 in component c65 (molar).
C RATES(66) is d/dt c66 in component c66 (molar).
C RATES(67) is d/dt c67 in component c67 (molar).
C RATES(68) is d/dt c68 in component c68 (molar).
C RATES(69) is d/dt c69 in component c69 (molar).
C RATES(70) is d/dt c70 in component c70 (molar).
C RATES(71) is d/dt c71 in component c71 (molar).
C RATES(72) is d/dt c72 in component c72 (molar).
C RATES(73) is d/dt c73 in component c73 (molar).
C RATES(74) is d/dt c74 in component c74 (molar).
C RATES(75) is d/dt c75 in component c75 (molar).
C RATES(76) is d/dt c76 in component c76 (molar).
C RATES(77) is d/dt c77 in component c77 (molar).
C RATES(78) is d/dt c78 in component c78 (molar).
C RATES(79) is d/dt c79 in component c79 (molar).
C RATES(80) is d/dt c80 in component c80 (molar).
C RATES(81) is d/dt c81 in component c81 (molar).
C RATES(82) is d/dt c82 in component c82 (molar).
C RATES(83) is d/dt c83 in component c83 (molar).
C RATES(84) is d/dt c84 in component c84 (molar).
C RATES(85) is d/dt c85 in component c85 (molar).
C RATES(86) is d/dt c86 in component c86 (molar).
C RATES(87) is d/dt c87 in component c87 (molar).
C RATES(88) is d/dt c88 in component c88 (molar).
C RATES(89) is d/dt c89 in component c89 (molar).
C RATES(90) is d/dt c90 in component c90 (molar).
C RATES(91) is d/dt c91 in component c91 (molar).
C RATES(92) is d/dt c92 in component c92 (molar).
C RATES(93) is d/dt c93 in component c93 (molar).
C RATES(94) is d/dt c94 in component c94 (molar).
C RATES(95) is d/dt c95 in component c95 (molar).
C RATES(96) is d/dt c96 in component c96 (molar).
C RATES(97) is d/dt c97 in component c97 (molar).
C RATES(98) is d/dt c98 in component c98 (molar).
C RATES(99) is d/dt c99 in component c99 (molar).
C RATES(100) is d/dt c100 in component c100 (molar).
C RATES(101) is d/dt c101 in component c101 (molar).
C RATES(102) is d/dt c102 in component c102 (molar).
C RATES(103) is d/dt c103 in component c103 (molar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 5.0E-8
      STATES(2) = 5.0E4
      STATES(3) = 0.0
      STATES(4) = 0.0
      STATES(5) = 0.0
      STATES(6) = 0.0
      STATES(7) = 0.0
      STATES(8) = 0.0
      STATES(9) = 0.0
      STATES(10) = 0.0
      STATES(11) = 0.0
      STATES(12) = 8.1E4
      STATES(13) = 0.0
      STATES(14) = 1.2E4
      STATES(15) = 0.0
      STATES(16) = 0.0
      STATES(17) = 0.0
      STATES(18) = 0.0
      STATES(19) = 0.0
      STATES(20) = 0.0
      STATES(21) = 0.0
      STATES(22) = 1.1E4
      STATES(23) = 0.0
      STATES(24) = 2.63E4
      STATES(25) = 0.0
      STATES(26) = 7.2E4
      STATES(27) = 0.0
      STATES(28) = 0.0
      STATES(29) = 0.0
      STATES(30) = 4.0E4
      STATES(31) = 1.01E5
      STATES(32) = 0.0
      STATES(33) = 0.0
      STATES(34) = 0.0
      STATES(35) = 0.0
      STATES(36) = 0.0
      STATES(37) = 0.0
      STATES(38) = 0.0
      STATES(39) = 0.0
      STATES(40) = 0.0
      STATES(41) = 4.0E4
      STATES(42) = 0.0
      STATES(43) = 0.0
      STATES(44) = 4.0E4
      STATES(45) = 0.0
      STATES(46) = 0.0
      STATES(47) = 2.1E7
      STATES(48) = 0.0
      STATES(49) = 0.0
      STATES(50) = 0.0
      STATES(51) = 0.0
      STATES(52) = 0.0
      STATES(53) = 4.0E4
      STATES(54) = 0.0
      STATES(55) = 2.21E7
      STATES(56) = 0.0
      STATES(57) = 0.0
      STATES(58) = 0.0
      STATES(59) = 0.0
      STATES(60) = 1.0E7
      STATES(61) = 0.0
      STATES(62) = 0.0
      STATES(63) = 0.0
      STATES(64) = 0.0
      STATES(65) = 0.0
      STATES(66) = 0.0
      STATES(67) = 0.0
      STATES(68) = 0.0
      STATES(69) = 0.0
      STATES(70) = 0.0
      STATES(71) = 0.0
      STATES(72) = 0.0
      STATES(73) = 0.0
      STATES(74) = 0.0
      STATES(75) = 0.0
      STATES(76) = 0.0
      STATES(77) = 0.0
      STATES(78) = 0.0
      STATES(79) = 0.0
      STATES(80) = 0.0
      STATES(81) = 0.0
      STATES(82) = 0.0
      STATES(83) = 0.0
      STATES(84) = 0.0
      STATES(85) = 0.0
      STATES(86) = 0.0
      STATES(87) = 0.0
      STATES(88) = 0.0
      STATES(89) = 0.0
      STATES(90) = 0.0
      STATES(91) = 0.0
      STATES(92) = 0.0
      STATES(93) = 0.0
      STATES(94) = 0.0
      STATES(95) = 0.0
      STATES(96) = 0.0
      STATES(97) = 0.0
      STATES(98) = 0.0
      STATES(99) = 0.0
      STATES(100) = 0.0
      STATES(101) = 0.0
      STATES(102) = 0.0
      STATES(103) = 0.0
      CONSTS(1) = 3.0E7
      CONSTS(2) = 3.84E-3
      CONSTS(3) = 1.66E-5
      CONSTS(4) = 1E-1
      CONSTS(5) = 1E0
      CONSTS(6) = 1E-2
      CONSTS(7) = 1.73E-7
      CONSTS(8) = 1.66E-3
      CONSTS(9) = 0E0
      CONSTS(10) = 1.46E-2
      CONSTS(11) = 5E-4
      CONSTS(12) = 5E-3
      CONSTS(13) = 1.66E-6
      CONSTS(14) = 2E-1
      CONSTS(15) = 5.43E-2
      CONSTS(16) = 1.1E-2
      CONSTS(17) = 2.17E0
      CONSTS(18) = 0E0
      CONSTS(19) = 1E4
      CONSTS(20) = 0E0
      CONSTS(21) = 1.66E-5
      CONSTS(22) = 2.75E-1
      CONSTS(23) = 1.66E-5
      CONSTS(24) = 6E-2
      CONSTS(25) = 2.5E-5
      CONSTS(26) = 1.3E0
      CONSTS(27) = 1.66E-7
      CONSTS(28) = 5E-1
      CONSTS(29) = 3.5E-6
      CONSTS(30) = 4E-1
      CONSTS(31) = 3.66E-7
      CONSTS(32) = 2.3E-2
      CONSTS(33) = 3.5E-5
      CONSTS(34) = 1E-1
      CONSTS(35) = 6E0
      CONSTS(36) = 6E-2
      CONSTS(37) = 5.5E-1
      CONSTS(38) = 1.66E-5
      CONSTS(39) = 2.14E-2
      CONSTS(40) = 1.66E-6
      CONSTS(41) = 5.3E-3
      CONSTS(42) = 1.17E-6
      CONSTS(43) = 1E0
      CONSTS(44) = 4E-7
      CONSTS(45) = 1E-1
      CONSTS(46) = 3.5E-5
      CONSTS(47) = 2E-1
      CONSTS(48) = 7.5E-6
      CONSTS(49) = 3E-2
      CONSTS(50) = 7.5E-6
      CONSTS(51) = 1.5E-3
      CONSTS(52) = 5E-3
      CONSTS(53) = 0E0
      CONSTS(54) = 1.5E-6
      CONSTS(55) = 3E-1
      CONSTS(56) = 5E-5
      CONSTS(57) = 6.4E-2
      CONSTS(58) = 5E-5
      CONSTS(59) = 4.29E-2
      CONSTS(60) = 1.18E-4
      CONSTS(61) = 2E-1
      CONSTS(62) = 0E0
      CONSTS(63) = 1E0
      CONSTS(64) = 1.95E-5
      CONSTS(65) = 3.3E-2
      CONSTS(66) = 0E0
      CONSTS(67) = 3.5E0
      CONSTS(68) = 0E0
      CONSTS(69) = 2.9E0
      CONSTS(70) = 2.38E-5
      CONSTS(71) = 8E-1
      CONSTS(72) = 0E0
      CONSTS(73) = 5.68E-2
      CONSTS(74) = 4.5E-7
      CONSTS(75) = 5E-1
      CONSTS(76) = 8.91E-5
      CONSTS(77) = 1.83E-2
      CONSTS(78) = 0E0
      CONSTS(79) = 1.6E1
      CONSTS(80) = 0E0
      CONSTS(81) = 5.7E0
      CONSTS(82) = 2.35E-5
      CONSTS(83) = 6E-1
      CONSTS(84) = 0E0
      CONSTS(85) = 2.46E-1
      CONSTS(86) = 8.33E-6
      CONSTS(87) = 5E-1
      CONSTS(88) = 5.5E-3
      CONSTS(89) = 0E0
      CONSTS(90) = 6.7E-4
      CONSTS(91) = 0E0
      CONSTS(92) = 1.66E-7
      CONSTS(93) = 2E0
      CONSTS(94) = 0E0
      CONSTS(95) = 1E-4
      CONSTS(95) = 0.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = CONSTS(95)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(2) -  CONSTS(2)*STATES(3)
      ALGBRC(2) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      RATES(3) = ALGBRC(1) -  2.00000*ALGBRC(2)
      ALGBRC(3) =  CONSTS(5)*STATES(4)*1.00000 -  CONSTS(6)*STATES(5)
      RATES(4) = ALGBRC(2) - ALGBRC(3)
      ALGBRC(4) =  CONSTS(7)*STATES(23)*STATES(12) -  CONSTS(8)*STATES(7)
      ALGBRC(5) =  CONSTS(9)*STATES(18)*STATES(9) -  CONSTS(10)*STATES(7)
      RATES(7) = ALGBRC(4)+ALGBRC(5)
      ALGBRC(6) =  CONSTS(11)*STATES(2) -  CONSTS(12)*STATES(6)
      ALGBRC(7) =  CONSTS(17)*1.00000 -  CONSTS(18)*STATES(2)
      RATES(2) = ALGBRC(7) - (ALGBRC(1)+ALGBRC(6))
      ALGBRC(8) =  CONSTS(11)*STATES(5) -  CONSTS(12)*STATES(8)
      ALGBRC(9) =  CONSTS(13)*STATES(5)*STATES(14) -  CONSTS(14)*STATES(15)
      RATES(5) = ALGBRC(3) - (ALGBRC(8)+ALGBRC(9))
      ALGBRC(11) =  CONSTS(15)*STATES(6)*STATES(16) -  CONSTS(16)*STATES(10)
      ALGBRC(12) =  CONSTS(88)*STATES(6) -  CONSTS(89)*STATES(86)
      RATES(6) = ALGBRC(6) - (ALGBRC(11)+ALGBRC(12))
      ALGBRC(13) =  CONSTS(3)*STATES(10)*STATES(10) -  CONSTS(4)*STATES(11)
      RATES(10) = ALGBRC(11) -  2.00000*ALGBRC(13)
      ALGBRC(14) =  CONSTS(90)*STATES(16) -  CONSTS(91)*STATES(13)
      RATES(13) = ALGBRC(14)
      RATES(16) = - (ALGBRC(11)+ALGBRC(14))
      RATES(86) = ALGBRC(12)
      ALGBRC(15) =  CONSTS(5)*STATES(11) -  CONSTS(6)*STATES(8)
      RATES(11) = ALGBRC(13) - ALGBRC(15)
      ALGBRC(16) =  CONSTS(13)*STATES(8)*STATES(14) -  CONSTS(14)*STATES(17)
      RATES(14) = - (ALGBRC(9)+ALGBRC(16))
      ALGBRC(18) =  CONSTS(21)*STATES(22)*STATES(15) -  CONSTS(22)*STATES(23)
      ALGBRC(10) =  CONSTS(11)*STATES(23) -  CONSTS(12)*STATES(18)
      ALGBRC(19) =  CONSTS(23)*STATES(24)*STATES(23) -  CONSTS(24)*STATES(25)
      RATES(23) = ALGBRC(18) - (ALGBRC(4)+ALGBRC(10)+ALGBRC(19))
      ALGBRC(21) =  CONSTS(27)*STATES(28)*STATES(25) -  CONSTS(28)*STATES(27)
      ALGBRC(29) =  CONSTS(27)*STATES(35)*STATES(28) -  CONSTS(28)*STATES(36)
      ALGBRC(30) =  CONSTS(40)*STATES(28)*STATES(41) -  CONSTS(41)*STATES(42)
      RATES(28) = - (ALGBRC(21)+ALGBRC(29)+ALGBRC(30))
      ALGBRC(31) =  CONSTS(42)*STATES(43)*STATES(45) -  CONSTS(43)*STATES(42)
      RATES(42) = ALGBRC(30)+ALGBRC(31)
      ALGBRC(22) =  CONSTS(29)*STATES(25)*STATES(43) -  CONSTS(30)*STATES(29)
      ALGBRC(32) =  CONSTS(29)*STATES(35)*STATES(43) -  CONSTS(30)*STATES(37)
      RATES(43) = - (ALGBRC(31)+ALGBRC(22)+ALGBRC(32))
      ALGBRC(45) =  CONSTS(62)*STATES(41)*STATES(44) -  CONSTS(63)*STATES(46)
      ALGBRC(44) =  CONSTS(60)*STATES(44)*STATES(45) -  CONSTS(61)*STATES(46)
      RATES(46) = ALGBRC(44)+ALGBRC(45)
      ALGBRC(46) =  CONSTS(64)*STATES(47)*STATES(45) -  CONSTS(65)*STATES(48)
      ALGBRC(47) =  CONSTS(66)*STATES(49)*STATES(45) -  CONSTS(67)*STATES(48)
      RATES(48) = ALGBRC(46)+ALGBRC(47)
      ALGBRC(48) =  CONSTS(64)*STATES(49)*STATES(45) -  CONSTS(65)*STATES(50)
      ALGBRC(49) =  CONSTS(68)*STATES(51)*STATES(45) -  CONSTS(69)*STATES(50)
      RATES(45) = - (ALGBRC(31)+ALGBRC(44)+ALGBRC(46)+ALGBRC(47)+ALGBRC(48)+ALGBRC(49))
      RATES(50) = ALGBRC(48)+ALGBRC(49)
      ALGBRC(51) =  CONSTS(72)*STATES(49)*STATES(53) -  CONSTS(73)*STATES(52)
      ALGBRC(50) =  CONSTS(70)*STATES(51)*STATES(53) -  CONSTS(71)*STATES(52)
      RATES(52) = ALGBRC(50)+ALGBRC(51)
      ALGBRC(52) =  CONSTS(74)*STATES(53)*STATES(49) -  CONSTS(75)*STATES(54)
      RATES(49) = - (ALGBRC(47)+ALGBRC(48)+ALGBRC(51)+ALGBRC(52))
      ALGBRC(53) =  CONSTS(72)*STATES(47)*STATES(53) -  CONSTS(73)*STATES(54)
      RATES(54) = ALGBRC(52)+ALGBRC(53)
      ALGBRC(54) =  CONSTS(76)*STATES(55)*STATES(51) -  CONSTS(77)*STATES(56)
      ALGBRC(55) =  CONSTS(78)*STATES(51)*STATES(57) -  CONSTS(79)*STATES(56)
      RATES(56) = ALGBRC(54)+ALGBRC(55)
      ALGBRC(56) =  CONSTS(76)*STATES(51)*STATES(57) -  CONSTS(77)*STATES(58)
      ALGBRC(57) =  CONSTS(80)*STATES(59)*STATES(51) -  CONSTS(81)*STATES(58)
      RATES(51) = - (ALGBRC(49)+ALGBRC(55)+ALGBRC(57)+ALGBRC(50)+ALGBRC(54)+ALGBRC(56))
      RATES(58) = ALGBRC(56)+ALGBRC(57)
      ALGBRC(59) =  CONSTS(84)*STATES(57)*STATES(60) -  CONSTS(85)*STATES(61)
      ALGBRC(58) =  CONSTS(82)*STATES(59)*STATES(60) -  CONSTS(83)*STATES(61)
      RATES(61) = ALGBRC(58)+ALGBRC(59)
      ALGBRC(60) =  CONSTS(86)*STATES(60)*STATES(57) -  CONSTS(87)*STATES(62)
      RATES(57) = - (ALGBRC(55)+ALGBRC(56)+ALGBRC(59)+ALGBRC(60))
      ALGBRC(61) =  CONSTS(84)*STATES(55)*STATES(60) -  CONSTS(85)*STATES(62)
      RATES(62) = ALGBRC(60)+ALGBRC(61)
      ALGBRC(62) =  CONSTS(88)*STATES(8) -  CONSTS(89)*STATES(87)
      RATES(8) = (ALGBRC(8)+ALGBRC(15)) - (ALGBRC(16)+ALGBRC(62))
      ALGBRC(24) =  CONSTS(33)*STATES(31)*STATES(15) -  CONSTS(34)*STATES(32)
      ALGBRC(69) =  CONSTS(33)*STATES(31)*STATES(17) -  CONSTS(34)*STATES(63)
      ALGBRC(38) =  CONSTS(52)*STATES(40) -  CONSTS(53)*STATES(31)
      RATES(31) = ALGBRC(38) - (ALGBRC(24)+ALGBRC(69))
      ALGBRC(63) =  CONSTS(21)*STATES(17)*STATES(22) -  CONSTS(22)*STATES(18)
      ALGBRC(26) =  CONSTS(21)*STATES(22)*STATES(33) -  CONSTS(37)*STATES(34)
      ALGBRC(37) =  CONSTS(50)*STATES(24)*STATES(22) -  CONSTS(51)*STATES(30)
      ALGBRC(40) =  CONSTS(21)*STATES(22)*STATES(40) -  CONSTS(37)*STATES(39)
      ALGBRC(71) =  CONSTS(21)*STATES(22)*STATES(64) -  CONSTS(37)*STATES(65)
      RATES(22) = - (ALGBRC(18)+ALGBRC(26)+ALGBRC(37)+ALGBRC(40)+ALGBRC(63)+ALGBRC(71))
      ALGBRC(66) =  CONSTS(27)*STATES(69)*STATES(19) -  CONSTS(28)*STATES(20)
      ALGBRC(75) =  CONSTS(40)*STATES(69)*STATES(41) -  CONSTS(41)*STATES(70)
      ALGBRC(74) =  CONSTS(27)*STATES(66)*STATES(69) -  CONSTS(28)*STATES(67)
      RATES(69) = - (ALGBRC(66)+ALGBRC(74)+ALGBRC(75))
      ALGBRC(76) =  CONSTS(42)*STATES(71)*STATES(72) -  CONSTS(43)*STATES(70)
      RATES(70) = ALGBRC(75)+ALGBRC(76)
      ALGBRC(67) =  CONSTS(29)*STATES(71)*STATES(19) -  CONSTS(30)*STATES(21)
      ALGBRC(77) =  CONSTS(29)*STATES(71)*STATES(66) -  CONSTS(30)*STATES(68)
      RATES(71) = - (ALGBRC(67)+ALGBRC(76)+ALGBRC(77))
      ALGBRC(65) =  CONSTS(25)*STATES(26)*STATES(19) -  CONSTS(26)*STATES(20)
      ALGBRC(68) =  CONSTS(31)*STATES(19)*STATES(26) -  CONSTS(32)*STATES(21)
      ALGBRC(20) =  CONSTS(25)*STATES(26)*STATES(25) -  CONSTS(26)*STATES(27)
      ALGBRC(23) =  CONSTS(31)*STATES(25)*STATES(26) -  CONSTS(32)*STATES(29)
      ALGBRC(28) =  CONSTS(25)*STATES(26)*STATES(35) -  CONSTS(26)*STATES(36)
      ALGBRC(33) =  CONSTS(31)*STATES(35)*STATES(26) -  CONSTS(32)*STATES(37)
      ALGBRC(78) =  CONSTS(31)*STATES(66)*STATES(26) -  CONSTS(32)*STATES(68)
      ALGBRC(73) =  CONSTS(25)*STATES(26)*STATES(66) -  CONSTS(26)*STATES(67)
      RATES(26) = - (ALGBRC(20)+ALGBRC(23)+ALGBRC(28)+ALGBRC(33)+ALGBRC(65)+ALGBRC(68)+ALGBRC(73)+ALGBRC(78))
      ALGBRC(34) =  CONSTS(44)*STATES(38)*STATES(15) -  CONSTS(45)*STATES(35)
      ALGBRC(79) =  CONSTS(44)*STATES(17)*STATES(38) -  CONSTS(45)*STATES(66)
      ALGBRC(42) =  CONSTS(56)*STATES(24)*STATES(39) -  CONSTS(57)*STATES(38)
      ALGBRC(35) =  CONSTS(46)*STATES(40)*STATES(30) -  CONSTS(47)*STATES(38)
      RATES(38) = (ALGBRC(35)+ALGBRC(42)) - (ALGBRC(34)+ALGBRC(79))
      ALGBRC(39) =  CONSTS(54)*STATES(15)*STATES(40) -  CONSTS(55)*STATES(33)
      ALGBRC(81) =  CONSTS(54)*STATES(17)*STATES(40) -  CONSTS(55)*STATES(64)
      RATES(40) = - (ALGBRC(35)+ALGBRC(38)+ALGBRC(39)+ALGBRC(40)+ALGBRC(81))
      ALGBRC(41) =  CONSTS(54)*STATES(15)*STATES(39) -  CONSTS(55)*STATES(34)
      ALGBRC(82) =  CONSTS(54)*STATES(17)*STATES(39) -  CONSTS(55)*STATES(65)
      RATES(39) = ALGBRC(40) - (ALGBRC(41)+ALGBRC(42)+ALGBRC(82))
      ALGBRC(36) =  CONSTS(48)*STATES(15)*STATES(30) -  CONSTS(49)*STATES(25)
      ALGBRC(80) =  CONSTS(48)*STATES(17)*STATES(30) -  CONSTS(49)*STATES(19)
      ALGBRC(43) =  CONSTS(58)*STATES(30)*STATES(33) -  CONSTS(59)*STATES(35)
      ALGBRC(83) =  CONSTS(58)*STATES(30)*STATES(64) -  CONSTS(59)*STATES(66)
      RATES(30) = ALGBRC(37) - (ALGBRC(35)+ALGBRC(36)+ALGBRC(43)+ALGBRC(80)+ALGBRC(83))
      ALGBRC(85) =  CONSTS(62)*STATES(41)*STATES(44) -  CONSTS(63)*STATES(73)
      RATES(41) = - (ALGBRC(30)+ALGBRC(45)+ALGBRC(75)+ALGBRC(85))
      ALGBRC(84) =  CONSTS(60)*STATES(44)*STATES(72) -  CONSTS(61)*STATES(73)
      RATES(44) = - (ALGBRC(44)+ALGBRC(45)+ALGBRC(84)+ALGBRC(85))
      RATES(73) = ALGBRC(84)+ALGBRC(85)
      ALGBRC(86) =  CONSTS(64)*STATES(47)*STATES(72) -  CONSTS(65)*STATES(74)
      ALGBRC(87) =  CONSTS(66)*STATES(75)*STATES(72) -  CONSTS(67)*STATES(74)
      RATES(74) = ALGBRC(86)+ALGBRC(87)
      ALGBRC(88) =  CONSTS(64)*STATES(72)*STATES(75) -  CONSTS(65)*STATES(76)
      ALGBRC(89) =  CONSTS(68)*STATES(72)*STATES(77) -  CONSTS(69)*STATES(76)
      RATES(72) = - (ALGBRC(76)+ALGBRC(84)+ALGBRC(86)+ALGBRC(87)+ALGBRC(88)+ALGBRC(89))
      RATES(76) = ALGBRC(88)+ALGBRC(89)
      ALGBRC(90) =  CONSTS(70)*STATES(77)*STATES(53) -  CONSTS(71)*STATES(78)
      ALGBRC(92) =  CONSTS(72)*STATES(75)*STATES(53) -  CONSTS(73)*STATES(78)
      RATES(78) = ALGBRC(90)+ALGBRC(92)
      ALGBRC(94) =  CONSTS(74)*STATES(53)*STATES(75) -  CONSTS(75)*STATES(79)
      RATES(75) = - (ALGBRC(87)+ALGBRC(88)+ALGBRC(92)+ALGBRC(94))
      ALGBRC(91) =  CONSTS(76)*STATES(55)*STATES(77) -  CONSTS(77)*STATES(80)
      ALGBRC(93) =  CONSTS(78)*STATES(81)*STATES(77) -  CONSTS(79)*STATES(80)
      RATES(80) = ALGBRC(91)+ALGBRC(93)
      ALGBRC(96) =  CONSTS(72)*STATES(47)*STATES(53) -  CONSTS(73)*STATES(79)
      RATES(47) = - (ALGBRC(46)+ALGBRC(53)+ALGBRC(86)+ALGBRC(96))
      RATES(53) = - (ALGBRC(50)+ALGBRC(51)+ALGBRC(52)+ALGBRC(53)+ALGBRC(90)+ALGBRC(92)+ALGBRC(94)+ALGBRC(96))
      RATES(79) = ALGBRC(94)+ALGBRC(96)
      ALGBRC(95) =  CONSTS(76)*STATES(77)*STATES(81) -  CONSTS(77)*STATES(82)
      ALGBRC(97) =  CONSTS(80)*STATES(83)*STATES(77) -  CONSTS(81)*STATES(82)
      RATES(77) = - (ALGBRC(89)+ALGBRC(90)+ALGBRC(91)+ALGBRC(93)+ALGBRC(95)+ALGBRC(97))
      RATES(82) = ALGBRC(95)+ALGBRC(97)
      ALGBRC(99) =  CONSTS(11)*STATES(15) -  CONSTS(12)*STATES(17)
      RATES(15) = ALGBRC(9) - (ALGBRC(18)+ALGBRC(24)+ALGBRC(34)+ALGBRC(36)+ALGBRC(39)+ALGBRC(41)+ALGBRC(99))
      ALGBRC(98) =  CONSTS(82)*STATES(83)*STATES(60) -  CONSTS(83)*STATES(84)
      ALGBRC(100) =  CONSTS(84)*STATES(81)*STATES(60) -  CONSTS(85)*STATES(84)
      RATES(84) = ALGBRC(98)+ALGBRC(100)
      ALGBRC(25) =  CONSTS(35)*STATES(32) -  CONSTS(36)*STATES(33)
      ALGBRC(101) =  CONSTS(11)*STATES(32) -  CONSTS(12)*STATES(63)
      RATES(32) = ALGBRC(24) - (ALGBRC(25)+ALGBRC(101))
      ALGBRC(102) =  CONSTS(86)*STATES(60)*STATES(81) -  CONSTS(87)*STATES(85)
      RATES(81) = - (ALGBRC(93)+ALGBRC(95)+ALGBRC(100)+ALGBRC(102))
      ALGBRC(103) =  CONSTS(11)*STATES(33) -  CONSTS(12)*STATES(64)
      RATES(33) = (ALGBRC(25)+ALGBRC(39)) - (ALGBRC(26)+ALGBRC(43)+ALGBRC(103))
      ALGBRC(104) =  CONSTS(84)*STATES(55)*STATES(60) -  CONSTS(85)*STATES(85)
      RATES(55) = - (ALGBRC(54)+ALGBRC(61)+ALGBRC(91)+ALGBRC(104))
      RATES(60) = - (ALGBRC(58)+ALGBRC(59)+ALGBRC(60)+ALGBRC(61)+ALGBRC(98)+ALGBRC(100)+ALGBRC(102)+ALGBRC(104))
      RATES(85) = ALGBRC(102)+ALGBRC(104)
      ALGBRC(107) =  CONSTS(9)*STATES(9)*STATES(19) -  CONSTS(10)*STATES(88)
      ALGBRC(106) =  CONSTS(7)*STATES(25)*STATES(12) -  CONSTS(8)*STATES(88)
      RATES(88) = ALGBRC(106)+ALGBRC(107)
      ALGBRC(109) =  CONSTS(7)*STATES(27)*STATES(12) -  CONSTS(8)*STATES(89)
      ALGBRC(108) =  CONSTS(11)*STATES(27) -  CONSTS(12)*STATES(20)
      RATES(27) = (ALGBRC(20)+ALGBRC(21)) - (ALGBRC(108)+ALGBRC(109))
      ALGBRC(110) =  CONSTS(9)*STATES(9)*STATES(20) -  CONSTS(10)*STATES(89)
      RATES(89) = ALGBRC(109)+ALGBRC(110)
      ALGBRC(112) =  CONSTS(7)*STATES(29)*STATES(12) -  CONSTS(8)*STATES(90)
      ALGBRC(111) =  CONSTS(11)*STATES(29) -  CONSTS(12)*STATES(21)
      RATES(29) = (ALGBRC(22)+ALGBRC(23)) - (ALGBRC(111)+ALGBRC(112))
      ALGBRC(113) =  CONSTS(9)*STATES(9)*STATES(21) -  CONSTS(10)*STATES(90)
      RATES(90) = ALGBRC(112)+ALGBRC(113)
      ALGBRC(115) =  CONSTS(7)*STATES(34)*STATES(12) -  CONSTS(8)*STATES(91)
      ALGBRC(27) =  CONSTS(38)*STATES(24)*STATES(34) -  CONSTS(39)*STATES(35)
      ALGBRC(114) =  CONSTS(11)*STATES(34) -  CONSTS(12)*STATES(65)
      RATES(34) = (ALGBRC(41)+ALGBRC(26)) - (ALGBRC(27)+ALGBRC(114)+ALGBRC(115))
      ALGBRC(116) =  CONSTS(9)*STATES(9)*STATES(65) -  CONSTS(10)*STATES(91)
      RATES(91) = ALGBRC(115)+ALGBRC(116)
      ALGBRC(119) =  CONSTS(9)*STATES(9)*STATES(66) -  CONSTS(10)*STATES(92)
      ALGBRC(118) =  CONSTS(7)*STATES(35)*STATES(12) -  CONSTS(8)*STATES(92)
      RATES(92) = ALGBRC(118)+ALGBRC(119)
      ALGBRC(121) =  CONSTS(7)*STATES(36)*STATES(12) -  CONSTS(8)*STATES(93)
      ALGBRC(120) =  CONSTS(11)*STATES(36) -  CONSTS(12)*STATES(67)
      RATES(36) = (ALGBRC(29)+ALGBRC(28)) - (ALGBRC(120)+ALGBRC(121))
      ALGBRC(122) =  CONSTS(9)*STATES(9)*STATES(67) -  CONSTS(10)*STATES(93)
      RATES(93) = ALGBRC(121)+ALGBRC(122)
      ALGBRC(17) =  CONSTS(19)*STATES(9) -  CONSTS(20)*STATES(12)
      ALGBRC(124) =  CONSTS(7)*STATES(37)*STATES(12) -  CONSTS(8)*STATES(94)
      RATES(12) = ALGBRC(17) - (ALGBRC(4)+ALGBRC(106)+ALGBRC(109)+ALGBRC(112)+ALGBRC(115)+ALGBRC(118)+ALGBRC(121)+ALGBRC(124))
      ALGBRC(123) =  CONSTS(11)*STATES(37) -  CONSTS(12)*STATES(68)
      RATES(37) = (ALGBRC(32)+ALGBRC(33)) - (ALGBRC(123)+ALGBRC(124))
      ALGBRC(125) =  CONSTS(9)*STATES(68)*STATES(9) -  CONSTS(10)*STATES(94)
      RATES(9) = - (ALGBRC(5)+ALGBRC(17)+ALGBRC(107)+ALGBRC(110)+ALGBRC(113)+ALGBRC(116)+ALGBRC(119)+ALGBRC(122)+ALGBRC(125))
      RATES(94) = ALGBRC(124)+ALGBRC(125)
      ALGBRC(105) =  CONSTS(11)*STATES(25) -  CONSTS(12)*STATES(19)
      ALGBRC(126) =  CONSTS(92)*STATES(59)*STATES(25) -  CONSTS(93)*STATES(95)
      RATES(25) = (ALGBRC(19)+ALGBRC(36)) - (ALGBRC(20)+ALGBRC(21)+ALGBRC(22)+ALGBRC(23)+ALGBRC(105)+ALGBRC(106)+ALGBRC(126))
      ALGBRC(117) =  CONSTS(11)*STATES(35) -  CONSTS(12)*STATES(66)
      ALGBRC(128) =  CONSTS(92)*STATES(59)*STATES(35) -  CONSTS(93)*STATES(97)
      RATES(35) = (ALGBRC(27)+ALGBRC(34)+ALGBRC(43)) - (ALGBRC(117)+ALGBRC(118)+ALGBRC(128)+ALGBRC(28)+ALGBRC(29)+ALGBRC(32)+ALGBRC(33))
      ALGBRC(130) =  CONSTS(88)*STATES(17) -  CONSTS(89)*STATES(87)
      RATES(17) = (ALGBRC(16)+ALGBRC(99)) - (ALGBRC(63)+ALGBRC(69)+ALGBRC(79)+ALGBRC(80)+ALGBRC(81)+ALGBRC(82)+ALGBRC(130))
      ALGBRC(64) =  CONSTS(23)*STATES(24)*STATES(18) -  CONSTS(24)*STATES(19)
      ALGBRC(132) =  CONSTS(88)*STATES(18) -  CONSTS(89)*STATES(87)
      RATES(18) = (ALGBRC(10)+ALGBRC(63)) - (ALGBRC(5)+ALGBRC(64)+ALGBRC(132))
      ALGBRC(72) =  CONSTS(38)*STATES(24)*STATES(65) -  CONSTS(39)*STATES(66)
      ALGBRC(131) =  CONSTS(92)*STATES(59)*STATES(24) -  CONSTS(93)*STATES(101)
      ALGBRC(133) =  CONSTS(92)*STATES(83)*STATES(24) -  CONSTS(93)*STATES(102)
      RATES(24) = - (ALGBRC(19)+ALGBRC(27)+ALGBRC(37)+ALGBRC(42)+ALGBRC(64)+ALGBRC(72)+ALGBRC(131)+ALGBRC(133))
      ALGBRC(127) =  CONSTS(92)*STATES(83)*STATES(19) -  CONSTS(93)*STATES(96)
      ALGBRC(134) =  CONSTS(88)*STATES(19) -  CONSTS(89)*STATES(87)
      RATES(19) = (ALGBRC(80)+ALGBRC(64)+ALGBRC(105)) - (ALGBRC(65)+ALGBRC(66)+ALGBRC(67)+ALGBRC(68)+ALGBRC(107)+ALGBRC(127)+ALGBRC(134))
      ALGBRC(135) =  CONSTS(94)*STATES(59)*STATES(99) -  CONSTS(95)*STATES(95)
      RATES(95) = ALGBRC(126)+ALGBRC(135)
      ALGBRC(136) =  CONSTS(88)*STATES(20) -  CONSTS(89)*STATES(87)
      RATES(20) = (ALGBRC(65)+ALGBRC(66)+ALGBRC(108)) - (ALGBRC(110)+ALGBRC(136))
      ALGBRC(137) =  CONSTS(94)*STATES(59)*STATES(99) -  CONSTS(95)*STATES(97)
      RATES(97) = ALGBRC(128)+ALGBRC(137)
      RATES(99) = - (ALGBRC(135)+ALGBRC(137))
      ALGBRC(138) =  CONSTS(88)*STATES(21) -  CONSTS(89)*STATES(87)
      RATES(21) = (ALGBRC(67)+ALGBRC(68)+ALGBRC(111)) - (ALGBRC(113)+ALGBRC(138))
      ALGBRC(139) =  CONSTS(94)*STATES(59)*STATES(103) -  CONSTS(95)*STATES(101)
      RATES(59) = - (ALGBRC(57)+ALGBRC(58)+ALGBRC(126)+ALGBRC(128)+ALGBRC(131)+ALGBRC(135)+ALGBRC(137)+ALGBRC(139))
      RATES(101) = ALGBRC(131)+ALGBRC(139)
      ALGBRC(70) =  CONSTS(35)*STATES(63) -  CONSTS(36)*STATES(64)
      ALGBRC(140) =  CONSTS(88)*STATES(63) -  CONSTS(89)*STATES(87)
      RATES(63) = (ALGBRC(69)+ALGBRC(101)) - (ALGBRC(70)+ALGBRC(140))
      ALGBRC(141) =  CONSTS(94)*STATES(83)*STATES(100) -  CONSTS(95)*STATES(96)
      RATES(96) = ALGBRC(127)+ALGBRC(141)
      ALGBRC(142) =  CONSTS(88)*STATES(64) -  CONSTS(89)*STATES(87)
      RATES(64) = (ALGBRC(70)+ALGBRC(81)+ALGBRC(103)) - (ALGBRC(71)+ALGBRC(83)+ALGBRC(142))
      ALGBRC(129) =  CONSTS(92)*STATES(83)*STATES(66) -  CONSTS(93)*STATES(98)
      ALGBRC(143) =  CONSTS(94)*STATES(83)*STATES(100) -  CONSTS(95)*STATES(98)
      RATES(98) = ALGBRC(129)+ALGBRC(143)
      RATES(100) = - (ALGBRC(141)+ALGBRC(143))
      ALGBRC(144) =  CONSTS(88)*STATES(65) -  CONSTS(89)*STATES(87)
      RATES(65) = (ALGBRC(71)+ALGBRC(82)+ALGBRC(114)) - (ALGBRC(72)+ALGBRC(116)+ALGBRC(144))
      ALGBRC(145) =  CONSTS(94)*STATES(83)*STATES(103) -  CONSTS(95)*STATES(102)
      RATES(83) = - (ALGBRC(97)+ALGBRC(98)+ALGBRC(127)+ALGBRC(129)+ALGBRC(133)+ALGBRC(141)+ALGBRC(143)+ALGBRC(145))
      RATES(102) = ALGBRC(133)+ALGBRC(145)
      RATES(103) = - (ALGBRC(139)+ALGBRC(145))
      ALGBRC(146) =  CONSTS(88)*STATES(66) -  CONSTS(89)*STATES(87)
      RATES(66) = (ALGBRC(72)+ALGBRC(79)+ALGBRC(83)+ALGBRC(117)) - (ALGBRC(73)+ALGBRC(74)+ALGBRC(77)+ALGBRC(78)+ALGBRC(119)+ALGBRC(129)+ALGBRC(146))
      ALGBRC(147) =  CONSTS(88)*STATES(67) -  CONSTS(89)*STATES(87)
      RATES(67) = (ALGBRC(73)+ALGBRC(74)+ALGBRC(120)) - (ALGBRC(122)+ALGBRC(147))
      ALGBRC(148) =  CONSTS(88)*STATES(68) -  CONSTS(89)*STATES(87)
      RATES(68) = (ALGBRC(77)+ALGBRC(78)+ALGBRC(123)) - (ALGBRC(125)+ALGBRC(148))
      RATES(87) = ALGBRC(62)+ALGBRC(130)+ALGBRC(132)+ALGBRC(134)+ALGBRC(136)+ALGBRC(138)+ALGBRC(140)+ALGBRC(142)+ALGBRC(144)+ALGBRC(146)+ALGBRC(147)+ALGBRC(148)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(2) -  CONSTS(2)*STATES(3)
      ALGBRC(2) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      ALGBRC(3) =  CONSTS(5)*STATES(4)*1.00000 -  CONSTS(6)*STATES(5)
      ALGBRC(4) =  CONSTS(7)*STATES(23)*STATES(12) -  CONSTS(8)*STATES(7)
      ALGBRC(5) =  CONSTS(9)*STATES(18)*STATES(9) -  CONSTS(10)*STATES(7)
      ALGBRC(6) =  CONSTS(11)*STATES(2) -  CONSTS(12)*STATES(6)
      ALGBRC(7) =  CONSTS(17)*1.00000 -  CONSTS(18)*STATES(2)
      ALGBRC(8) =  CONSTS(11)*STATES(5) -  CONSTS(12)*STATES(8)
      ALGBRC(9) =  CONSTS(13)*STATES(5)*STATES(14) -  CONSTS(14)*STATES(15)
      ALGBRC(11) =  CONSTS(15)*STATES(6)*STATES(16) -  CONSTS(16)*STATES(10)
      ALGBRC(12) =  CONSTS(88)*STATES(6) -  CONSTS(89)*STATES(86)
      ALGBRC(13) =  CONSTS(3)*STATES(10)*STATES(10) -  CONSTS(4)*STATES(11)
      ALGBRC(14) =  CONSTS(90)*STATES(16) -  CONSTS(91)*STATES(13)
      ALGBRC(15) =  CONSTS(5)*STATES(11) -  CONSTS(6)*STATES(8)
      ALGBRC(16) =  CONSTS(13)*STATES(8)*STATES(14) -  CONSTS(14)*STATES(17)
      ALGBRC(18) =  CONSTS(21)*STATES(22)*STATES(15) -  CONSTS(22)*STATES(23)
      ALGBRC(10) =  CONSTS(11)*STATES(23) -  CONSTS(12)*STATES(18)
      ALGBRC(19) =  CONSTS(23)*STATES(24)*STATES(23) -  CONSTS(24)*STATES(25)
      ALGBRC(21) =  CONSTS(27)*STATES(28)*STATES(25) -  CONSTS(28)*STATES(27)
      ALGBRC(29) =  CONSTS(27)*STATES(35)*STATES(28) -  CONSTS(28)*STATES(36)
      ALGBRC(30) =  CONSTS(40)*STATES(28)*STATES(41) -  CONSTS(41)*STATES(42)
      ALGBRC(31) =  CONSTS(42)*STATES(43)*STATES(45) -  CONSTS(43)*STATES(42)
      ALGBRC(22) =  CONSTS(29)*STATES(25)*STATES(43) -  CONSTS(30)*STATES(29)
      ALGBRC(32) =  CONSTS(29)*STATES(35)*STATES(43) -  CONSTS(30)*STATES(37)
      ALGBRC(45) =  CONSTS(62)*STATES(41)*STATES(44) -  CONSTS(63)*STATES(46)
      ALGBRC(44) =  CONSTS(60)*STATES(44)*STATES(45) -  CONSTS(61)*STATES(46)
      ALGBRC(46) =  CONSTS(64)*STATES(47)*STATES(45) -  CONSTS(65)*STATES(48)
      ALGBRC(47) =  CONSTS(66)*STATES(49)*STATES(45) -  CONSTS(67)*STATES(48)
      ALGBRC(48) =  CONSTS(64)*STATES(49)*STATES(45) -  CONSTS(65)*STATES(50)
      ALGBRC(49) =  CONSTS(68)*STATES(51)*STATES(45) -  CONSTS(69)*STATES(50)
      ALGBRC(51) =  CONSTS(72)*STATES(49)*STATES(53) -  CONSTS(73)*STATES(52)
      ALGBRC(50) =  CONSTS(70)*STATES(51)*STATES(53) -  CONSTS(71)*STATES(52)
      ALGBRC(52) =  CONSTS(74)*STATES(53)*STATES(49) -  CONSTS(75)*STATES(54)
      ALGBRC(53) =  CONSTS(72)*STATES(47)*STATES(53) -  CONSTS(73)*STATES(54)
      ALGBRC(54) =  CONSTS(76)*STATES(55)*STATES(51) -  CONSTS(77)*STATES(56)
      ALGBRC(55) =  CONSTS(78)*STATES(51)*STATES(57) -  CONSTS(79)*STATES(56)
      ALGBRC(56) =  CONSTS(76)*STATES(51)*STATES(57) -  CONSTS(77)*STATES(58)
      ALGBRC(57) =  CONSTS(80)*STATES(59)*STATES(51) -  CONSTS(81)*STATES(58)
      ALGBRC(59) =  CONSTS(84)*STATES(57)*STATES(60) -  CONSTS(85)*STATES(61)
      ALGBRC(58) =  CONSTS(82)*STATES(59)*STATES(60) -  CONSTS(83)*STATES(61)
      ALGBRC(60) =  CONSTS(86)*STATES(60)*STATES(57) -  CONSTS(87)*STATES(62)
      ALGBRC(61) =  CONSTS(84)*STATES(55)*STATES(60) -  CONSTS(85)*STATES(62)
      ALGBRC(62) =  CONSTS(88)*STATES(8) -  CONSTS(89)*STATES(87)
      ALGBRC(24) =  CONSTS(33)*STATES(31)*STATES(15) -  CONSTS(34)*STATES(32)
      ALGBRC(69) =  CONSTS(33)*STATES(31)*STATES(17) -  CONSTS(34)*STATES(63)
      ALGBRC(38) =  CONSTS(52)*STATES(40) -  CONSTS(53)*STATES(31)
      ALGBRC(63) =  CONSTS(21)*STATES(17)*STATES(22) -  CONSTS(22)*STATES(18)
      ALGBRC(26) =  CONSTS(21)*STATES(22)*STATES(33) -  CONSTS(37)*STATES(34)
      ALGBRC(37) =  CONSTS(50)*STATES(24)*STATES(22) -  CONSTS(51)*STATES(30)
      ALGBRC(40) =  CONSTS(21)*STATES(22)*STATES(40) -  CONSTS(37)*STATES(39)
      ALGBRC(71) =  CONSTS(21)*STATES(22)*STATES(64) -  CONSTS(37)*STATES(65)
      ALGBRC(66) =  CONSTS(27)*STATES(69)*STATES(19) -  CONSTS(28)*STATES(20)
      ALGBRC(75) =  CONSTS(40)*STATES(69)*STATES(41) -  CONSTS(41)*STATES(70)
      ALGBRC(74) =  CONSTS(27)*STATES(66)*STATES(69) -  CONSTS(28)*STATES(67)
      ALGBRC(76) =  CONSTS(42)*STATES(71)*STATES(72) -  CONSTS(43)*STATES(70)
      ALGBRC(67) =  CONSTS(29)*STATES(71)*STATES(19) -  CONSTS(30)*STATES(21)
      ALGBRC(77) =  CONSTS(29)*STATES(71)*STATES(66) -  CONSTS(30)*STATES(68)
      ALGBRC(65) =  CONSTS(25)*STATES(26)*STATES(19) -  CONSTS(26)*STATES(20)
      ALGBRC(68) =  CONSTS(31)*STATES(19)*STATES(26) -  CONSTS(32)*STATES(21)
      ALGBRC(20) =  CONSTS(25)*STATES(26)*STATES(25) -  CONSTS(26)*STATES(27)
      ALGBRC(23) =  CONSTS(31)*STATES(25)*STATES(26) -  CONSTS(32)*STATES(29)
      ALGBRC(28) =  CONSTS(25)*STATES(26)*STATES(35) -  CONSTS(26)*STATES(36)
      ALGBRC(33) =  CONSTS(31)*STATES(35)*STATES(26) -  CONSTS(32)*STATES(37)
      ALGBRC(78) =  CONSTS(31)*STATES(66)*STATES(26) -  CONSTS(32)*STATES(68)
      ALGBRC(73) =  CONSTS(25)*STATES(26)*STATES(66) -  CONSTS(26)*STATES(67)
      ALGBRC(34) =  CONSTS(44)*STATES(38)*STATES(15) -  CONSTS(45)*STATES(35)
      ALGBRC(79) =  CONSTS(44)*STATES(17)*STATES(38) -  CONSTS(45)*STATES(66)
      ALGBRC(42) =  CONSTS(56)*STATES(24)*STATES(39) -  CONSTS(57)*STATES(38)
      ALGBRC(35) =  CONSTS(46)*STATES(40)*STATES(30) -  CONSTS(47)*STATES(38)
      ALGBRC(39) =  CONSTS(54)*STATES(15)*STATES(40) -  CONSTS(55)*STATES(33)
      ALGBRC(81) =  CONSTS(54)*STATES(17)*STATES(40) -  CONSTS(55)*STATES(64)
      ALGBRC(41) =  CONSTS(54)*STATES(15)*STATES(39) -  CONSTS(55)*STATES(34)
      ALGBRC(82) =  CONSTS(54)*STATES(17)*STATES(39) -  CONSTS(55)*STATES(65)
      ALGBRC(36) =  CONSTS(48)*STATES(15)*STATES(30) -  CONSTS(49)*STATES(25)
      ALGBRC(80) =  CONSTS(48)*STATES(17)*STATES(30) -  CONSTS(49)*STATES(19)
      ALGBRC(43) =  CONSTS(58)*STATES(30)*STATES(33) -  CONSTS(59)*STATES(35)
      ALGBRC(83) =  CONSTS(58)*STATES(30)*STATES(64) -  CONSTS(59)*STATES(66)
      ALGBRC(85) =  CONSTS(62)*STATES(41)*STATES(44) -  CONSTS(63)*STATES(73)
      ALGBRC(84) =  CONSTS(60)*STATES(44)*STATES(72) -  CONSTS(61)*STATES(73)
      ALGBRC(86) =  CONSTS(64)*STATES(47)*STATES(72) -  CONSTS(65)*STATES(74)
      ALGBRC(87) =  CONSTS(66)*STATES(75)*STATES(72) -  CONSTS(67)*STATES(74)
      ALGBRC(88) =  CONSTS(64)*STATES(72)*STATES(75) -  CONSTS(65)*STATES(76)
      ALGBRC(89) =  CONSTS(68)*STATES(72)*STATES(77) -  CONSTS(69)*STATES(76)
      ALGBRC(90) =  CONSTS(70)*STATES(77)*STATES(53) -  CONSTS(71)*STATES(78)
      ALGBRC(92) =  CONSTS(72)*STATES(75)*STATES(53) -  CONSTS(73)*STATES(78)
      ALGBRC(94) =  CONSTS(74)*STATES(53)*STATES(75) -  CONSTS(75)*STATES(79)
      ALGBRC(91) =  CONSTS(76)*STATES(55)*STATES(77) -  CONSTS(77)*STATES(80)
      ALGBRC(93) =  CONSTS(78)*STATES(81)*STATES(77) -  CONSTS(79)*STATES(80)
      ALGBRC(96) =  CONSTS(72)*STATES(47)*STATES(53) -  CONSTS(73)*STATES(79)
      ALGBRC(95) =  CONSTS(76)*STATES(77)*STATES(81) -  CONSTS(77)*STATES(82)
      ALGBRC(97) =  CONSTS(80)*STATES(83)*STATES(77) -  CONSTS(81)*STATES(82)
      ALGBRC(99) =  CONSTS(11)*STATES(15) -  CONSTS(12)*STATES(17)
      ALGBRC(98) =  CONSTS(82)*STATES(83)*STATES(60) -  CONSTS(83)*STATES(84)
      ALGBRC(100) =  CONSTS(84)*STATES(81)*STATES(60) -  CONSTS(85)*STATES(84)
      ALGBRC(25) =  CONSTS(35)*STATES(32) -  CONSTS(36)*STATES(33)
      ALGBRC(101) =  CONSTS(11)*STATES(32) -  CONSTS(12)*STATES(63)
      ALGBRC(102) =  CONSTS(86)*STATES(60)*STATES(81) -  CONSTS(87)*STATES(85)
      ALGBRC(103) =  CONSTS(11)*STATES(33) -  CONSTS(12)*STATES(64)
      ALGBRC(104) =  CONSTS(84)*STATES(55)*STATES(60) -  CONSTS(85)*STATES(85)
      ALGBRC(107) =  CONSTS(9)*STATES(9)*STATES(19) -  CONSTS(10)*STATES(88)
      ALGBRC(106) =  CONSTS(7)*STATES(25)*STATES(12) -  CONSTS(8)*STATES(88)
      ALGBRC(109) =  CONSTS(7)*STATES(27)*STATES(12) -  CONSTS(8)*STATES(89)
      ALGBRC(108) =  CONSTS(11)*STATES(27) -  CONSTS(12)*STATES(20)
      ALGBRC(110) =  CONSTS(9)*STATES(9)*STATES(20) -  CONSTS(10)*STATES(89)
      ALGBRC(112) =  CONSTS(7)*STATES(29)*STATES(12) -  CONSTS(8)*STATES(90)
      ALGBRC(111) =  CONSTS(11)*STATES(29) -  CONSTS(12)*STATES(21)
      ALGBRC(113) =  CONSTS(9)*STATES(9)*STATES(21) -  CONSTS(10)*STATES(90)
      ALGBRC(115) =  CONSTS(7)*STATES(34)*STATES(12) -  CONSTS(8)*STATES(91)
      ALGBRC(27) =  CONSTS(38)*STATES(24)*STATES(34) -  CONSTS(39)*STATES(35)
      ALGBRC(114) =  CONSTS(11)*STATES(34) -  CONSTS(12)*STATES(65)
      ALGBRC(116) =  CONSTS(9)*STATES(9)*STATES(65) -  CONSTS(10)*STATES(91)
      ALGBRC(119) =  CONSTS(9)*STATES(9)*STATES(66) -  CONSTS(10)*STATES(92)
      ALGBRC(118) =  CONSTS(7)*STATES(35)*STATES(12) -  CONSTS(8)*STATES(92)
      ALGBRC(121) =  CONSTS(7)*STATES(36)*STATES(12) -  CONSTS(8)*STATES(93)
      ALGBRC(120) =  CONSTS(11)*STATES(36) -  CONSTS(12)*STATES(67)
      ALGBRC(122) =  CONSTS(9)*STATES(9)*STATES(67) -  CONSTS(10)*STATES(93)
      ALGBRC(17) =  CONSTS(19)*STATES(9) -  CONSTS(20)*STATES(12)
      ALGBRC(124) =  CONSTS(7)*STATES(37)*STATES(12) -  CONSTS(8)*STATES(94)
      ALGBRC(123) =  CONSTS(11)*STATES(37) -  CONSTS(12)*STATES(68)
      ALGBRC(125) =  CONSTS(9)*STATES(68)*STATES(9) -  CONSTS(10)*STATES(94)
      ALGBRC(105) =  CONSTS(11)*STATES(25) -  CONSTS(12)*STATES(19)
      ALGBRC(126) =  CONSTS(92)*STATES(59)*STATES(25) -  CONSTS(93)*STATES(95)
      ALGBRC(117) =  CONSTS(11)*STATES(35) -  CONSTS(12)*STATES(66)
      ALGBRC(128) =  CONSTS(92)*STATES(59)*STATES(35) -  CONSTS(93)*STATES(97)
      ALGBRC(130) =  CONSTS(88)*STATES(17) -  CONSTS(89)*STATES(87)
      ALGBRC(64) =  CONSTS(23)*STATES(24)*STATES(18) -  CONSTS(24)*STATES(19)
      ALGBRC(132) =  CONSTS(88)*STATES(18) -  CONSTS(89)*STATES(87)
      ALGBRC(72) =  CONSTS(38)*STATES(24)*STATES(65) -  CONSTS(39)*STATES(66)
      ALGBRC(131) =  CONSTS(92)*STATES(59)*STATES(24) -  CONSTS(93)*STATES(101)
      ALGBRC(133) =  CONSTS(92)*STATES(83)*STATES(24) -  CONSTS(93)*STATES(102)
      ALGBRC(127) =  CONSTS(92)*STATES(83)*STATES(19) -  CONSTS(93)*STATES(96)
      ALGBRC(134) =  CONSTS(88)*STATES(19) -  CONSTS(89)*STATES(87)
      ALGBRC(135) =  CONSTS(94)*STATES(59)*STATES(99) -  CONSTS(95)*STATES(95)
      ALGBRC(136) =  CONSTS(88)*STATES(20) -  CONSTS(89)*STATES(87)
      ALGBRC(137) =  CONSTS(94)*STATES(59)*STATES(99) -  CONSTS(95)*STATES(97)
      ALGBRC(138) =  CONSTS(88)*STATES(21) -  CONSTS(89)*STATES(87)
      ALGBRC(139) =  CONSTS(94)*STATES(59)*STATES(103) -  CONSTS(95)*STATES(101)
      ALGBRC(70) =  CONSTS(35)*STATES(63) -  CONSTS(36)*STATES(64)
      ALGBRC(140) =  CONSTS(88)*STATES(63) -  CONSTS(89)*STATES(87)
      ALGBRC(141) =  CONSTS(94)*STATES(83)*STATES(100) -  CONSTS(95)*STATES(96)
      ALGBRC(142) =  CONSTS(88)*STATES(64) -  CONSTS(89)*STATES(87)
      ALGBRC(129) =  CONSTS(92)*STATES(83)*STATES(66) -  CONSTS(93)*STATES(98)
      ALGBRC(143) =  CONSTS(94)*STATES(83)*STATES(100) -  CONSTS(95)*STATES(98)
      ALGBRC(144) =  CONSTS(88)*STATES(65) -  CONSTS(89)*STATES(87)
      ALGBRC(145) =  CONSTS(94)*STATES(83)*STATES(103) -  CONSTS(95)*STATES(102)
      ALGBRC(146) =  CONSTS(88)*STATES(66) -  CONSTS(89)*STATES(87)
      ALGBRC(147) =  CONSTS(88)*STATES(67) -  CONSTS(89)*STATES(87)
      ALGBRC(148) =  CONSTS(88)*STATES(68) -  CONSTS(89)*STATES(87)
      RETURN
      END