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