Generated Code
The following is python code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
# Size of variable arrays: sizeAlgebraic = 124 sizeStates = 19 sizeConstants = 58 from math import * from numpy import * def createLegends(): legend_states = [""] * sizeStates legend_rates = [""] * sizeStates legend_algebraic = [""] * sizeAlgebraic legend_voi = "" legend_constants = [""] * sizeConstants legend_voi = "t in component environment (second)" legend_constants[0] = "R in component environment (J_per_K_per_mol)" legend_constants[1] = "T in component environment (kelvin)" legend_constants[2] = "F in component environment (C_per_mol)" legend_constants[3] = "C_m in component environment (fF)" legend_states[0] = "q_Na_o in component environment (fmol)" legend_states[1] = "q_Na_i in component environment (fmol)" legend_states[2] = "q_S000_Na in component environment (fmol)" legend_states[3] = "q_S010_Na in component environment (fmol)" legend_states[4] = "q_S100_Na in component environment (fmol)" legend_states[5] = "q_S110_Na in component environment (fmol)" legend_states[6] = "q_S200_Na in component environment (fmol)" legend_states[7] = "q_S210_Na in component environment (fmol)" legend_states[8] = "q_S300_Na in component environment (fmol)" legend_states[9] = "q_S310_Na in component environment (fmol)" legend_states[10] = "q_S001_Na in component environment (fmol)" legend_states[11] = "q_S011_Na in component environment (fmol)" legend_states[12] = "q_S101_Na in component environment (fmol)" legend_states[13] = "q_S111_Na in component environment (fmol)" legend_states[14] = "q_S201_Na in component environment (fmol)" legend_states[15] = "q_S211_Na in component environment (fmol)" legend_states[16] = "q_S301_Na in component environment (fmol)" legend_states[17] = "q_S311_Na in component environment (fmol)" legend_algebraic[120] = "v_Na in component fast_Na (fmol_per_sec)" legend_states[18] = "q_mem in component environment (fC)" legend_algebraic[122] = "I_mem_Na in component fast_Na (fA)" legend_constants[4] = "K_000_Na in component Na_parameters (per_fmol)" legend_constants[5] = "K_001_Na in component Na_parameters (per_fmol)" legend_constants[6] = "K_010_Na in component Na_parameters (per_fmol)" legend_constants[7] = "K_011_Na in component Na_parameters (per_fmol)" legend_constants[8] = "K_100_Na in component Na_parameters (per_fmol)" legend_constants[9] = "K_101_Na in component Na_parameters (per_fmol)" legend_constants[10] = "K_110_Na in component Na_parameters (per_fmol)" legend_constants[11] = "K_111_Na in component Na_parameters (per_fmol)" legend_constants[12] = "K_200_Na in component Na_parameters (per_fmol)" legend_constants[13] = "K_201_Na in component Na_parameters (per_fmol)" legend_constants[14] = "K_210_Na in component Na_parameters (per_fmol)" legend_constants[15] = "K_211_Na in component Na_parameters (per_fmol)" legend_constants[16] = "K_300_Na in component Na_parameters (per_fmol)" legend_constants[17] = "K_301_Na in component Na_parameters (per_fmol)" legend_constants[18] = "K_310_Na in component Na_parameters (per_fmol)" legend_constants[19] = "K_311_Na in component Na_parameters (per_fmol)" legend_constants[20] = "K_Na_o in component Na_parameters (per_fmol)" legend_constants[21] = "K_Na_i in component Na_parameters (per_fmol)" legend_constants[22] = "kappa_Na in component Na_parameters (fmol_per_sec)" legend_constants[23] = "kappa_h000 in component Na_parameters (fmol_per_sec)" legend_constants[24] = "kappa_h001 in component Na_parameters (fmol_per_sec)" legend_constants[25] = "kappa_h100 in component Na_parameters (fmol_per_sec)" legend_constants[26] = "kappa_h101 in component Na_parameters (fmol_per_sec)" legend_constants[27] = "kappa_h200 in component Na_parameters (fmol_per_sec)" legend_constants[28] = "kappa_h201 in component Na_parameters (fmol_per_sec)" legend_constants[29] = "kappa_h300 in component Na_parameters (fmol_per_sec)" legend_constants[30] = "kappa_h301 in component Na_parameters (fmol_per_sec)" legend_constants[31] = "kappa_j000 in component Na_parameters (fmol_per_sec)" legend_constants[32] = "kappa_j010 in component Na_parameters (fmol_per_sec)" legend_constants[33] = "kappa_j100 in component Na_parameters (fmol_per_sec)" legend_constants[34] = "kappa_j110 in component Na_parameters (fmol_per_sec)" legend_constants[35] = "kappa_j200 in component Na_parameters (fmol_per_sec)" legend_constants[36] = "kappa_j210 in component Na_parameters (fmol_per_sec)" legend_constants[37] = "kappa_j300 in component Na_parameters (fmol_per_sec)" legend_constants[38] = "kappa_j310 in component Na_parameters (fmol_per_sec)" legend_constants[39] = "kappa_m000 in component Na_parameters (fmol_per_sec)" legend_constants[40] = "kappa_m001 in component Na_parameters (fmol_per_sec)" legend_constants[41] = "kappa_m010 in component Na_parameters (fmol_per_sec)" legend_constants[42] = "kappa_m011 in component Na_parameters (fmol_per_sec)" legend_constants[43] = "kappa_m100 in component Na_parameters (fmol_per_sec)" legend_constants[44] = "kappa_m101 in component Na_parameters (fmol_per_sec)" legend_constants[45] = "kappa_m110 in component Na_parameters (fmol_per_sec)" legend_constants[46] = "kappa_m111 in component Na_parameters (fmol_per_sec)" legend_constants[47] = "kappa_m200 in component Na_parameters (fmol_per_sec)" legend_constants[48] = "kappa_m201 in component Na_parameters (fmol_per_sec)" legend_constants[49] = "kappa_m210 in component Na_parameters (fmol_per_sec)" legend_constants[50] = "kappa_m211 in component Na_parameters (fmol_per_sec)" legend_constants[51] = "zNa in component Na_parameters (dimensionless)" legend_constants[52] = "z_fh in component Na_parameters (dimensionless)" legend_constants[53] = "z_fj in component Na_parameters (dimensionless)" legend_constants[54] = "z_fm in component Na_parameters (dimensionless)" legend_constants[55] = "z_rh in component Na_parameters (dimensionless)" legend_constants[56] = "z_rj in component Na_parameters (dimensionless)" legend_constants[57] = "z_rm in component Na_parameters (dimensionless)" legend_algebraic[16] = "mu_Na_o in component fast_Na (J_per_mol)" legend_algebraic[17] = "mu_Na_i in component fast_Na (J_per_mol)" legend_algebraic[1] = "mu_S000_Na in component fast_Na (J_per_mol)" legend_algebraic[2] = "mu_S010_Na in component fast_Na (J_per_mol)" legend_algebraic[3] = "mu_S100_Na in component fast_Na (J_per_mol)" legend_algebraic[4] = "mu_S110_Na in component fast_Na (J_per_mol)" legend_algebraic[5] = "mu_S200_Na in component fast_Na (J_per_mol)" legend_algebraic[6] = "mu_S210_Na in component fast_Na (J_per_mol)" legend_algebraic[7] = "mu_S300_Na in component fast_Na (J_per_mol)" legend_algebraic[8] = "mu_S310_Na in component fast_Na (J_per_mol)" legend_algebraic[9] = "mu_S001_Na in component fast_Na (J_per_mol)" legend_algebraic[10] = "mu_S011_Na in component fast_Na (J_per_mol)" legend_algebraic[11] = "mu_S101_Na in component fast_Na (J_per_mol)" legend_algebraic[12] = "mu_S111_Na in component fast_Na (J_per_mol)" legend_algebraic[13] = "mu_S201_Na in component fast_Na (J_per_mol)" legend_algebraic[14] = "mu_S211_Na in component fast_Na (J_per_mol)" legend_algebraic[15] = "mu_S301_Na in component fast_Na (J_per_mol)" legend_algebraic[18] = "mu_S311_Na in component fast_Na (J_per_mol)" legend_algebraic[48] = "Am_Na in component fast_Na (J_per_mol)" legend_algebraic[19] = "Af_Na in component fast_Na (J_per_mol)" legend_algebraic[49] = "Ar_Na in component fast_Na (J_per_mol)" legend_algebraic[20] = "Af_h000 in component fast_Na (J_per_mol)" legend_algebraic[50] = "Ar_h000 in component fast_Na (J_per_mol)" legend_algebraic[51] = "v_h000 in component fast_Na (fmol_per_sec)" legend_algebraic[21] = "Af_h001 in component fast_Na (J_per_mol)" legend_algebraic[52] = "Ar_h001 in component fast_Na (J_per_mol)" legend_algebraic[53] = "v_h001 in component fast_Na (fmol_per_sec)" legend_algebraic[22] = "Af_h100 in component fast_Na (J_per_mol)" legend_algebraic[54] = "Ar_h100 in component fast_Na (J_per_mol)" legend_algebraic[55] = "v_h100 in component fast_Na (fmol_per_sec)" legend_algebraic[23] = "Af_h101 in component fast_Na (J_per_mol)" legend_algebraic[56] = "Ar_h101 in component fast_Na (J_per_mol)" legend_algebraic[57] = "v_h101 in component fast_Na (fmol_per_sec)" legend_algebraic[24] = "Af_h200 in component fast_Na (J_per_mol)" legend_algebraic[58] = "Ar_h200 in component fast_Na (J_per_mol)" legend_algebraic[59] = "v_h200 in component fast_Na (fmol_per_sec)" legend_algebraic[25] = "Af_h201 in component fast_Na (J_per_mol)" legend_algebraic[60] = "Ar_h201 in component fast_Na (J_per_mol)" legend_algebraic[61] = "v_h201 in component fast_Na (fmol_per_sec)" legend_algebraic[26] = "Af_h300 in component fast_Na (J_per_mol)" legend_algebraic[62] = "Ar_h300 in component fast_Na (J_per_mol)" legend_algebraic[63] = "v_h300 in component fast_Na (fmol_per_sec)" legend_algebraic[27] = "Af_h301 in component fast_Na (J_per_mol)" legend_algebraic[64] = "Ar_h301 in component fast_Na (J_per_mol)" legend_algebraic[65] = "v_h301 in component fast_Na (fmol_per_sec)" legend_algebraic[28] = "Af_j000 in component fast_Na (J_per_mol)" legend_algebraic[66] = "Ar_j000 in component fast_Na (J_per_mol)" legend_algebraic[67] = "v_j000 in component fast_Na (fmol_per_sec)" legend_algebraic[29] = "Af_j010 in component fast_Na (J_per_mol)" legend_algebraic[68] = "Ar_j010 in component fast_Na (J_per_mol)" legend_algebraic[69] = "v_j010 in component fast_Na (fmol_per_sec)" legend_algebraic[30] = "Af_j100 in component fast_Na (J_per_mol)" legend_algebraic[70] = "Ar_j100 in component fast_Na (J_per_mol)" legend_algebraic[71] = "v_j100 in component fast_Na (fmol_per_sec)" legend_algebraic[31] = "Af_j110 in component fast_Na (J_per_mol)" legend_algebraic[72] = "Ar_j110 in component fast_Na (J_per_mol)" legend_algebraic[73] = "v_j110 in component fast_Na (fmol_per_sec)" legend_algebraic[32] = "Af_j200 in component fast_Na (J_per_mol)" legend_algebraic[74] = "Ar_j200 in component fast_Na (J_per_mol)" legend_algebraic[75] = "v_j200 in component fast_Na (fmol_per_sec)" legend_algebraic[33] = "Af_j210 in component fast_Na (J_per_mol)" legend_algebraic[76] = "Ar_j210 in component fast_Na (J_per_mol)" legend_algebraic[77] = "v_j210 in component fast_Na (fmol_per_sec)" legend_algebraic[34] = "Af_j300 in component fast_Na (J_per_mol)" legend_algebraic[78] = "Ar_j300 in component fast_Na (J_per_mol)" legend_algebraic[79] = "v_j300 in component fast_Na (fmol_per_sec)" legend_algebraic[35] = "Af_j310 in component fast_Na (J_per_mol)" legend_algebraic[80] = "Ar_j310 in component fast_Na (J_per_mol)" legend_algebraic[81] = "v_j310 in component fast_Na (fmol_per_sec)" legend_algebraic[36] = "Af_m000 in component fast_Na (J_per_mol)" legend_algebraic[82] = "Ar_m000 in component fast_Na (J_per_mol)" legend_algebraic[83] = "v_m000 in component fast_Na (fmol_per_sec)" legend_algebraic[37] = "Af_m001 in component fast_Na (J_per_mol)" legend_algebraic[84] = "Ar_m001 in component fast_Na (J_per_mol)" legend_algebraic[86] = "v_m001 in component fast_Na (fmol_per_sec)" legend_algebraic[38] = "Af_m010 in component fast_Na (J_per_mol)" legend_algebraic[87] = "Ar_m010 in component fast_Na (J_per_mol)" legend_algebraic[89] = "v_m010 in component fast_Na (fmol_per_sec)" legend_algebraic[39] = "Af_m011 in component fast_Na (J_per_mol)" legend_algebraic[90] = "Ar_m011 in component fast_Na (J_per_mol)" legend_algebraic[92] = "v_m011 in component fast_Na (fmol_per_sec)" legend_algebraic[40] = "Af_m100 in component fast_Na (J_per_mol)" legend_algebraic[93] = "Ar_m100 in component fast_Na (J_per_mol)" legend_algebraic[95] = "v_m100 in component fast_Na (fmol_per_sec)" legend_algebraic[41] = "Af_m101 in component fast_Na (J_per_mol)" legend_algebraic[96] = "Ar_m101 in component fast_Na (J_per_mol)" legend_algebraic[98] = "v_m101 in component fast_Na (fmol_per_sec)" legend_algebraic[42] = "Af_m110 in component fast_Na (J_per_mol)" legend_algebraic[99] = "Ar_m110 in component fast_Na (J_per_mol)" legend_algebraic[101] = "v_m110 in component fast_Na (fmol_per_sec)" legend_algebraic[43] = "Af_m111 in component fast_Na (J_per_mol)" legend_algebraic[102] = "Ar_m111 in component fast_Na (J_per_mol)" legend_algebraic[104] = "v_m111 in component fast_Na (fmol_per_sec)" legend_algebraic[44] = "Af_m200 in component fast_Na (J_per_mol)" legend_algebraic[105] = "Ar_m200 in component fast_Na (J_per_mol)" legend_algebraic[107] = "v_m200 in component fast_Na (fmol_per_sec)" legend_algebraic[45] = "Af_m201 in component fast_Na (J_per_mol)" legend_algebraic[108] = "Ar_m201 in component fast_Na (J_per_mol)" legend_algebraic[111] = "v_m201 in component fast_Na (fmol_per_sec)" legend_algebraic[46] = "Af_m210 in component fast_Na (J_per_mol)" legend_algebraic[112] = "Ar_m210 in component fast_Na (J_per_mol)" legend_algebraic[115] = "v_m210 in component fast_Na (fmol_per_sec)" legend_algebraic[47] = "Af_m211 in component fast_Na (J_per_mol)" legend_algebraic[116] = "Ar_m211 in component fast_Na (J_per_mol)" legend_algebraic[119] = "v_m211 in component fast_Na (fmol_per_sec)" legend_algebraic[85] = "v_S000_Na in component fast_Na (fmol_per_sec)" legend_algebraic[91] = "v_S010_Na in component fast_Na (fmol_per_sec)" legend_algebraic[97] = "v_S100_Na in component fast_Na (fmol_per_sec)" legend_algebraic[103] = "v_S110_Na in component fast_Na (fmol_per_sec)" legend_algebraic[109] = "v_S200_Na in component fast_Na (fmol_per_sec)" legend_algebraic[117] = "v_S210_Na in component fast_Na (fmol_per_sec)" legend_algebraic[110] = "v_S300_Na in component fast_Na (fmol_per_sec)" legend_algebraic[118] = "v_S310_Na in component fast_Na (fmol_per_sec)" legend_algebraic[88] = "v_S001_Na in component fast_Na (fmol_per_sec)" legend_algebraic[94] = "v_S011_Na in component fast_Na (fmol_per_sec)" legend_algebraic[100] = "v_S101_Na in component fast_Na (fmol_per_sec)" legend_algebraic[106] = "v_S111_Na in component fast_Na (fmol_per_sec)" legend_algebraic[113] = "v_S201_Na in component fast_Na (fmol_per_sec)" legend_algebraic[121] = "v_S211_Na in component fast_Na (fmol_per_sec)" legend_algebraic[114] = "v_S301_Na in component fast_Na (fmol_per_sec)" legend_algebraic[123] = "v_S311_Na in component fast_Na (fmol_per_sec)" legend_algebraic[0] = "V_mem in component fast_Na (volt)" legend_rates[0] = "d/dt q_Na_o in component environment (fmol)" legend_rates[1] = "d/dt q_Na_i in component environment (fmol)" legend_rates[18] = "d/dt q_mem in component environment (fC)" legend_rates[2] = "d/dt q_S000_Na in component environment (fmol)" legend_rates[3] = "d/dt q_S010_Na in component environment (fmol)" legend_rates[4] = "d/dt q_S100_Na in component environment (fmol)" legend_rates[5] = "d/dt q_S110_Na in component environment (fmol)" legend_rates[6] = "d/dt q_S200_Na in component environment (fmol)" legend_rates[7] = "d/dt q_S210_Na in component environment (fmol)" legend_rates[8] = "d/dt q_S300_Na in component environment (fmol)" legend_rates[9] = "d/dt q_S310_Na in component environment (fmol)" legend_rates[10] = "d/dt q_S001_Na in component environment (fmol)" legend_rates[11] = "d/dt q_S011_Na in component environment (fmol)" legend_rates[12] = "d/dt q_S101_Na in component environment (fmol)" legend_rates[13] = "d/dt q_S111_Na in component environment (fmol)" legend_rates[14] = "d/dt q_S201_Na in component environment (fmol)" legend_rates[15] = "d/dt q_S211_Na in component environment (fmol)" legend_rates[16] = "d/dt q_S301_Na in component environment (fmol)" legend_rates[17] = "d/dt q_S311_Na in component environment (fmol)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; constants[0] = 8.314 constants[1] = 310 constants[2] = 96485 constants[3] = 153400 states[0] = 725.48 states[1] = 380 states[2] = 1.72572E-05 states[3] = 1.91747E-06 states[4] = 5.75256E-06 states[5] = 6.39176E-07 states[6] = 6.39176E-07 states[7] = 7.10196E-08 states[8] = 2.36737E-08 states[9] = 2.63039E-09 states[10] = 1.91747E-06 states[11] = 2.13059E-07 states[12] = 6.39176E-07 states[13] = 7.10196E-08 states[14] = 7.10196E-08 states[15] = 7.89103E-09 states[16] = 2.63039E-09 states[17] = 2.92267E-10 states[18] = -13039 constants[4] = 0.0073376 constants[5] = 23718.0222 constants[6] = 17499.8767 constants[7] = 56566767727.4036 constants[8] = 1.5633e-05 constants[9] = 50.5317 constants[10] = 37.2839 constants[11] = 120516694.0958 constants[12] = 9.9918e-08 constants[13] = 0.32298 constants[14] = 0.2383 constants[15] = 770290.0911 constants[16] = 1.9159e-09 constants[17] = 0.006193 constants[18] = 0.0045694 constants[19] = 14770.0739 constants[20] = 0.017678 constants[21] = 0.0024108 constants[22] = 32.858 constants[23] = 0.045709 constants[24] = 1.4141e-08 constants[25] = 21.4544 constants[26] = 6.6373e-06 constants[27] = 3356.675 constants[28] = 0.0010384 constants[29] = 175057.5863 constants[30] = 0.054157 constants[31] = 0.017824 constants[32] = 7.4736e-09 constants[33] = 8.3662 constants[34] = 3.5079e-06 constants[35] = 1308.9427 constants[36] = 0.00054883 constants[37] = 68264.0876 constants[38] = 0.028623 constants[39] = 5117405.6443 constants[40] = 1.5832 constants[41] = 2.1457 constants[42] = 6.6381e-07 constants[43] = 1601300127.6001 constants[44] = 495.389 constants[45] = 671.4132 constants[46] = 0.00020771 constants[47] = 125266701374.7363 constants[48] = 38753.3515 constants[49] = 52523.3901 constants[50] = 0.016249 constants[51] = 1 constants[52] = -4.1892 constants[53] = -4.0381 constants[54] = 0.49541 constants[55] = 1.2995 constants[56] = 1.4281 constants[57] = -2.4284 return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic algebraic[0] = states[18]/constants[3] algebraic[20] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[52]*constants[2]*algebraic[0] algebraic[50] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[55]*constants[2]*algebraic[0] algebraic[51] = constants[23]*(exp(algebraic[20]/(constants[0]*constants[1]))-exp(algebraic[50]/(constants[0]*constants[1]))) algebraic[28] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[53]*constants[2]*algebraic[0] algebraic[66] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[56]*constants[2]*algebraic[0] algebraic[67] = constants[31]*(exp(algebraic[28]/(constants[0]*constants[1]))-exp(algebraic[66]/(constants[0]*constants[1]))) algebraic[36] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[54]*constants[2]*algebraic[0] algebraic[82] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[57]*constants[2]*algebraic[0] algebraic[83] = constants[39]*(exp(algebraic[36]/(constants[0]*constants[1]))-exp(algebraic[82]/(constants[0]*constants[1]))) algebraic[85] = (-algebraic[83]-algebraic[51])-algebraic[67] rates[2] = algebraic[85] algebraic[21] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[52]*constants[2]*algebraic[0] algebraic[52] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[55]*constants[2]*algebraic[0] algebraic[53] = constants[24]*(exp(algebraic[21]/(constants[0]*constants[1]))-exp(algebraic[52]/(constants[0]*constants[1]))) algebraic[37] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[54]*constants[2]*algebraic[0] algebraic[84] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[57]*constants[2]*algebraic[0] algebraic[86] = constants[40]*(exp(algebraic[37]/(constants[0]*constants[1]))-exp(algebraic[84]/(constants[0]*constants[1]))) algebraic[88] = (-algebraic[86]-algebraic[53])--algebraic[67] rates[10] = algebraic[88] algebraic[29] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[53]*constants[2]*algebraic[0] algebraic[68] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[56]*constants[2]*algebraic[0] algebraic[69] = constants[32]*(exp(algebraic[29]/(constants[0]*constants[1]))-exp(algebraic[68]/(constants[0]*constants[1]))) algebraic[38] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[54]*constants[2]*algebraic[0] algebraic[87] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[57]*constants[2]*algebraic[0] algebraic[89] = constants[41]*(exp(algebraic[38]/(constants[0]*constants[1]))-exp(algebraic[87]/(constants[0]*constants[1]))) algebraic[91] = (-algebraic[89]+algebraic[51])-algebraic[69] rates[3] = algebraic[91] algebraic[39] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[54]*constants[2]*algebraic[0] algebraic[90] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[57]*constants[2]*algebraic[0] algebraic[92] = constants[42]*(exp(algebraic[39]/(constants[0]*constants[1]))-exp(algebraic[90]/(constants[0]*constants[1]))) algebraic[94] = (-algebraic[92]+algebraic[53])--algebraic[69] rates[11] = algebraic[94] algebraic[22] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[52]*constants[2]*algebraic[0] algebraic[54] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[55]*constants[2]*algebraic[0] algebraic[55] = constants[25]*(exp(algebraic[22]/(constants[0]*constants[1]))-exp(algebraic[54]/(constants[0]*constants[1]))) algebraic[30] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[53]*constants[2]*algebraic[0] algebraic[70] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[56]*constants[2]*algebraic[0] algebraic[71] = constants[33]*(exp(algebraic[30]/(constants[0]*constants[1]))-exp(algebraic[70]/(constants[0]*constants[1]))) algebraic[40] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[54]*constants[2]*algebraic[0] algebraic[93] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[57]*constants[2]*algebraic[0] algebraic[95] = constants[43]*(exp(algebraic[40]/(constants[0]*constants[1]))-exp(algebraic[93]/(constants[0]*constants[1]))) algebraic[97] = ((algebraic[83]-algebraic[95])-algebraic[55])-algebraic[71] rates[4] = algebraic[97] algebraic[23] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[52]*constants[2]*algebraic[0] algebraic[56] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[55]*constants[2]*algebraic[0] algebraic[57] = constants[26]*(exp(algebraic[23]/(constants[0]*constants[1]))-exp(algebraic[56]/(constants[0]*constants[1]))) algebraic[41] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[54]*constants[2]*algebraic[0] algebraic[96] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[57]*constants[2]*algebraic[0] algebraic[98] = constants[44]*(exp(algebraic[41]/(constants[0]*constants[1]))-exp(algebraic[96]/(constants[0]*constants[1]))) algebraic[100] = ((algebraic[86]-algebraic[98])-algebraic[57])--algebraic[71] rates[12] = algebraic[100] algebraic[31] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[53]*constants[2]*algebraic[0] algebraic[72] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[56]*constants[2]*algebraic[0] algebraic[73] = constants[34]*(exp(algebraic[31]/(constants[0]*constants[1]))-exp(algebraic[72]/(constants[0]*constants[1]))) algebraic[42] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[54]*constants[2]*algebraic[0] algebraic[99] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[57]*constants[2]*algebraic[0] algebraic[101] = constants[45]*(exp(algebraic[42]/(constants[0]*constants[1]))-exp(algebraic[99]/(constants[0]*constants[1]))) algebraic[103] = ((algebraic[89]-algebraic[101])+algebraic[55])-algebraic[73] rates[5] = algebraic[103] algebraic[43] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[54]*constants[2]*algebraic[0] algebraic[102] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[57]*constants[2]*algebraic[0] algebraic[104] = constants[46]*(exp(algebraic[43]/(constants[0]*constants[1]))-exp(algebraic[102]/(constants[0]*constants[1]))) algebraic[106] = ((algebraic[92]-algebraic[104])+algebraic[57])--algebraic[73] rates[13] = algebraic[106] algebraic[24] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[52]*constants[2]*algebraic[0] algebraic[58] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[55]*constants[2]*algebraic[0] algebraic[59] = constants[27]*(exp(algebraic[24]/(constants[0]*constants[1]))-exp(algebraic[58]/(constants[0]*constants[1]))) algebraic[32] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[53]*constants[2]*algebraic[0] algebraic[74] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[56]*constants[2]*algebraic[0] algebraic[75] = constants[35]*(exp(algebraic[32]/(constants[0]*constants[1]))-exp(algebraic[74]/(constants[0]*constants[1]))) algebraic[44] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[54]*constants[2]*algebraic[0] algebraic[105] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[57]*constants[2]*algebraic[0] algebraic[107] = constants[47]*(exp(algebraic[44]/(constants[0]*constants[1]))-exp(algebraic[105]/(constants[0]*constants[1]))) algebraic[109] = ((algebraic[95]-algebraic[107])-algebraic[59])-algebraic[75] rates[6] = algebraic[109] algebraic[26] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[52]*constants[2]*algebraic[0] algebraic[62] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[55]*constants[2]*algebraic[0] algebraic[63] = constants[29]*(exp(algebraic[26]/(constants[0]*constants[1]))-exp(algebraic[62]/(constants[0]*constants[1]))) algebraic[34] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[53]*constants[2]*algebraic[0] algebraic[78] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[56]*constants[2]*algebraic[0] algebraic[79] = constants[37]*(exp(algebraic[34]/(constants[0]*constants[1]))-exp(algebraic[78]/(constants[0]*constants[1]))) algebraic[110] = (algebraic[107]-algebraic[63])-algebraic[79] rates[8] = algebraic[110] algebraic[25] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[52]*constants[2]*algebraic[0] algebraic[60] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[55]*constants[2]*algebraic[0] algebraic[61] = constants[28]*(exp(algebraic[25]/(constants[0]*constants[1]))-exp(algebraic[60]/(constants[0]*constants[1]))) algebraic[45] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[54]*constants[2]*algebraic[0] algebraic[108] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[57]*constants[2]*algebraic[0] algebraic[111] = constants[48]*(exp(algebraic[45]/(constants[0]*constants[1]))-exp(algebraic[108]/(constants[0]*constants[1]))) algebraic[113] = ((algebraic[98]-algebraic[111])-algebraic[61])--algebraic[75] rates[14] = algebraic[113] algebraic[27] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[52]*constants[2]*algebraic[0] algebraic[64] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[55]*constants[2]*algebraic[0] algebraic[65] = constants[30]*(exp(algebraic[27]/(constants[0]*constants[1]))-exp(algebraic[64]/(constants[0]*constants[1]))) algebraic[114] = (algebraic[111]-algebraic[65])--algebraic[79] rates[16] = algebraic[114] algebraic[33] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[53]*constants[2]*algebraic[0] algebraic[76] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[56]*constants[2]*algebraic[0] algebraic[77] = constants[36]*(exp(algebraic[33]/(constants[0]*constants[1]))-exp(algebraic[76]/(constants[0]*constants[1]))) algebraic[46] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[54]*constants[2]*algebraic[0] algebraic[112] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[57]*constants[2]*algebraic[0] algebraic[115] = constants[49]*(exp(algebraic[46]/(constants[0]*constants[1]))-exp(algebraic[112]/(constants[0]*constants[1]))) algebraic[117] = ((algebraic[101]-algebraic[115])+algebraic[59])-algebraic[77] rates[7] = algebraic[117] algebraic[35] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[53]*constants[2]*algebraic[0] algebraic[80] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[56]*constants[2]*algebraic[0] algebraic[81] = constants[38]*(exp(algebraic[35]/(constants[0]*constants[1]))-exp(algebraic[80]/(constants[0]*constants[1]))) algebraic[118] = (algebraic[115]+algebraic[63])-algebraic[81] rates[9] = algebraic[118] algebraic[48] = constants[51]*constants[2]*algebraic[0] algebraic[17] = constants[0]*constants[1]*log(constants[21]*states[1]) algebraic[18] = constants[0]*constants[1]*log(constants[19]*states[17]) algebraic[19] = algebraic[17]+constants[51]*constants[2]*algebraic[0]+algebraic[18] algebraic[16] = constants[0]*constants[1]*log(constants[20]*states[0]) algebraic[49] = algebraic[16]+algebraic[18] algebraic[120] = custom_piecewise([equal(algebraic[48] , 0.00000), constants[22]*(exp(algebraic[19]/(constants[0]*constants[1]))-exp(algebraic[49]/(constants[0]*constants[1]))) , True, (((constants[22]*algebraic[48])/(constants[0]*constants[1]))/(exp(algebraic[48]/(constants[0]*constants[1]))-1.00000))*(exp(algebraic[19]/(constants[0]*constants[1]))-exp(algebraic[49]/(constants[0]*constants[1])))]) rates[0] = algebraic[120] rates[1] = -algebraic[120] algebraic[47] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[54]*constants[2]*algebraic[0] algebraic[116] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[57]*constants[2]*algebraic[0] algebraic[119] = constants[50]*(exp(algebraic[47]/(constants[0]*constants[1]))-exp(algebraic[116]/(constants[0]*constants[1]))) algebraic[121] = ((algebraic[104]-algebraic[119])+algebraic[61])--algebraic[77] rates[15] = algebraic[121] algebraic[122] = constants[2]*((((((((((((((((((((((((((((((((((((((((((((((((((((((((-constants[51]*algebraic[120]+constants[57]*algebraic[83])-constants[54]*algebraic[83])+constants[57]*algebraic[95])-constants[54]*algebraic[95])+constants[57]*algebraic[89])-constants[54]*algebraic[89])+constants[57]*algebraic[101])-constants[54]*algebraic[101])+constants[57]*algebraic[115])-constants[54]*algebraic[115])+constants[55]*algebraic[51])-constants[52]*algebraic[51])+constants[55]*algebraic[59])-constants[52]*algebraic[59])+constants[56]*algebraic[71])-constants[53]*algebraic[71])+constants[56]*algebraic[75])-constants[53]*algebraic[75])+constants[56]*algebraic[79])-constants[53]*algebraic[79])+constants[57]*algebraic[107])-constants[54]*algebraic[107])+constants[55]*algebraic[55])-constants[52]*algebraic[55])+constants[55]*algebraic[63])-constants[52]*algebraic[63])+constants[56]*algebraic[67])-constants[53]*algebraic[67])+constants[57]*algebraic[86])-constants[54]*algebraic[86])+constants[57]*algebraic[98])-constants[54]*algebraic[98])+constants[57]*algebraic[92])-constants[54]*algebraic[92])+constants[57]*algebraic[104])-constants[54]*algebraic[104])+constants[57]*algebraic[119])-constants[54]*algebraic[119])+constants[55]*algebraic[53])-constants[52]*algebraic[53])+constants[55]*algebraic[61])-constants[52]*algebraic[61])+constants[56]*algebraic[69])-constants[53]*algebraic[69])+constants[56]*algebraic[73])-constants[53]*algebraic[73])+constants[56]*algebraic[77])-constants[53]*algebraic[77])+constants[56]*algebraic[81])-constants[53]*algebraic[81])+constants[57]*algebraic[111])-constants[54]*algebraic[111])+constants[55]*algebraic[57])-constants[52]*algebraic[57])+constants[55]*algebraic[65])-constants[52]*algebraic[65]) rates[18] = algebraic[122] algebraic[123] = (algebraic[119]+algebraic[65])-((-algebraic[81]+algebraic[120])-algebraic[120]) rates[17] = algebraic[123] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[0] = states[18]/constants[3] algebraic[20] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[52]*constants[2]*algebraic[0] algebraic[50] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[55]*constants[2]*algebraic[0] algebraic[51] = constants[23]*(exp(algebraic[20]/(constants[0]*constants[1]))-exp(algebraic[50]/(constants[0]*constants[1]))) algebraic[28] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[53]*constants[2]*algebraic[0] algebraic[66] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[56]*constants[2]*algebraic[0] algebraic[67] = constants[31]*(exp(algebraic[28]/(constants[0]*constants[1]))-exp(algebraic[66]/(constants[0]*constants[1]))) algebraic[36] = constants[0]*constants[1]*log(constants[4]*states[2])+constants[54]*constants[2]*algebraic[0] algebraic[82] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[57]*constants[2]*algebraic[0] algebraic[83] = constants[39]*(exp(algebraic[36]/(constants[0]*constants[1]))-exp(algebraic[82]/(constants[0]*constants[1]))) algebraic[85] = (-algebraic[83]-algebraic[51])-algebraic[67] algebraic[21] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[52]*constants[2]*algebraic[0] algebraic[52] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[55]*constants[2]*algebraic[0] algebraic[53] = constants[24]*(exp(algebraic[21]/(constants[0]*constants[1]))-exp(algebraic[52]/(constants[0]*constants[1]))) algebraic[37] = constants[0]*constants[1]*log(constants[5]*states[10])+constants[54]*constants[2]*algebraic[0] algebraic[84] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[57]*constants[2]*algebraic[0] algebraic[86] = constants[40]*(exp(algebraic[37]/(constants[0]*constants[1]))-exp(algebraic[84]/(constants[0]*constants[1]))) algebraic[88] = (-algebraic[86]-algebraic[53])--algebraic[67] algebraic[29] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[53]*constants[2]*algebraic[0] algebraic[68] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[56]*constants[2]*algebraic[0] algebraic[69] = constants[32]*(exp(algebraic[29]/(constants[0]*constants[1]))-exp(algebraic[68]/(constants[0]*constants[1]))) algebraic[38] = constants[0]*constants[1]*log(constants[6]*states[3])+constants[54]*constants[2]*algebraic[0] algebraic[87] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[57]*constants[2]*algebraic[0] algebraic[89] = constants[41]*(exp(algebraic[38]/(constants[0]*constants[1]))-exp(algebraic[87]/(constants[0]*constants[1]))) algebraic[91] = (-algebraic[89]+algebraic[51])-algebraic[69] algebraic[39] = constants[0]*constants[1]*log(constants[7]*states[11])+constants[54]*constants[2]*algebraic[0] algebraic[90] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[57]*constants[2]*algebraic[0] algebraic[92] = constants[42]*(exp(algebraic[39]/(constants[0]*constants[1]))-exp(algebraic[90]/(constants[0]*constants[1]))) algebraic[94] = (-algebraic[92]+algebraic[53])--algebraic[69] algebraic[22] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[52]*constants[2]*algebraic[0] algebraic[54] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[55]*constants[2]*algebraic[0] algebraic[55] = constants[25]*(exp(algebraic[22]/(constants[0]*constants[1]))-exp(algebraic[54]/(constants[0]*constants[1]))) algebraic[30] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[53]*constants[2]*algebraic[0] algebraic[70] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[56]*constants[2]*algebraic[0] algebraic[71] = constants[33]*(exp(algebraic[30]/(constants[0]*constants[1]))-exp(algebraic[70]/(constants[0]*constants[1]))) algebraic[40] = constants[0]*constants[1]*log(constants[8]*states[4])+constants[54]*constants[2]*algebraic[0] algebraic[93] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[57]*constants[2]*algebraic[0] algebraic[95] = constants[43]*(exp(algebraic[40]/(constants[0]*constants[1]))-exp(algebraic[93]/(constants[0]*constants[1]))) algebraic[97] = ((algebraic[83]-algebraic[95])-algebraic[55])-algebraic[71] algebraic[23] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[52]*constants[2]*algebraic[0] algebraic[56] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[55]*constants[2]*algebraic[0] algebraic[57] = constants[26]*(exp(algebraic[23]/(constants[0]*constants[1]))-exp(algebraic[56]/(constants[0]*constants[1]))) algebraic[41] = constants[0]*constants[1]*log(constants[9]*states[12])+constants[54]*constants[2]*algebraic[0] algebraic[96] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[57]*constants[2]*algebraic[0] algebraic[98] = constants[44]*(exp(algebraic[41]/(constants[0]*constants[1]))-exp(algebraic[96]/(constants[0]*constants[1]))) algebraic[100] = ((algebraic[86]-algebraic[98])-algebraic[57])--algebraic[71] algebraic[31] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[53]*constants[2]*algebraic[0] algebraic[72] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[56]*constants[2]*algebraic[0] algebraic[73] = constants[34]*(exp(algebraic[31]/(constants[0]*constants[1]))-exp(algebraic[72]/(constants[0]*constants[1]))) algebraic[42] = constants[0]*constants[1]*log(constants[10]*states[5])+constants[54]*constants[2]*algebraic[0] algebraic[99] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[57]*constants[2]*algebraic[0] algebraic[101] = constants[45]*(exp(algebraic[42]/(constants[0]*constants[1]))-exp(algebraic[99]/(constants[0]*constants[1]))) algebraic[103] = ((algebraic[89]-algebraic[101])+algebraic[55])-algebraic[73] algebraic[43] = constants[0]*constants[1]*log(constants[11]*states[13])+constants[54]*constants[2]*algebraic[0] algebraic[102] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[57]*constants[2]*algebraic[0] algebraic[104] = constants[46]*(exp(algebraic[43]/(constants[0]*constants[1]))-exp(algebraic[102]/(constants[0]*constants[1]))) algebraic[106] = ((algebraic[92]-algebraic[104])+algebraic[57])--algebraic[73] algebraic[24] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[52]*constants[2]*algebraic[0] algebraic[58] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[55]*constants[2]*algebraic[0] algebraic[59] = constants[27]*(exp(algebraic[24]/(constants[0]*constants[1]))-exp(algebraic[58]/(constants[0]*constants[1]))) algebraic[32] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[53]*constants[2]*algebraic[0] algebraic[74] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[56]*constants[2]*algebraic[0] algebraic[75] = constants[35]*(exp(algebraic[32]/(constants[0]*constants[1]))-exp(algebraic[74]/(constants[0]*constants[1]))) algebraic[44] = constants[0]*constants[1]*log(constants[12]*states[6])+constants[54]*constants[2]*algebraic[0] algebraic[105] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[57]*constants[2]*algebraic[0] algebraic[107] = constants[47]*(exp(algebraic[44]/(constants[0]*constants[1]))-exp(algebraic[105]/(constants[0]*constants[1]))) algebraic[109] = ((algebraic[95]-algebraic[107])-algebraic[59])-algebraic[75] algebraic[26] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[52]*constants[2]*algebraic[0] algebraic[62] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[55]*constants[2]*algebraic[0] algebraic[63] = constants[29]*(exp(algebraic[26]/(constants[0]*constants[1]))-exp(algebraic[62]/(constants[0]*constants[1]))) algebraic[34] = constants[0]*constants[1]*log(constants[16]*states[8])+constants[53]*constants[2]*algebraic[0] algebraic[78] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[56]*constants[2]*algebraic[0] algebraic[79] = constants[37]*(exp(algebraic[34]/(constants[0]*constants[1]))-exp(algebraic[78]/(constants[0]*constants[1]))) algebraic[110] = (algebraic[107]-algebraic[63])-algebraic[79] algebraic[25] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[52]*constants[2]*algebraic[0] algebraic[60] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[55]*constants[2]*algebraic[0] algebraic[61] = constants[28]*(exp(algebraic[25]/(constants[0]*constants[1]))-exp(algebraic[60]/(constants[0]*constants[1]))) algebraic[45] = constants[0]*constants[1]*log(constants[13]*states[14])+constants[54]*constants[2]*algebraic[0] algebraic[108] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[57]*constants[2]*algebraic[0] algebraic[111] = constants[48]*(exp(algebraic[45]/(constants[0]*constants[1]))-exp(algebraic[108]/(constants[0]*constants[1]))) algebraic[113] = ((algebraic[98]-algebraic[111])-algebraic[61])--algebraic[75] algebraic[27] = constants[0]*constants[1]*log(constants[17]*states[16])+constants[52]*constants[2]*algebraic[0] algebraic[64] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[55]*constants[2]*algebraic[0] algebraic[65] = constants[30]*(exp(algebraic[27]/(constants[0]*constants[1]))-exp(algebraic[64]/(constants[0]*constants[1]))) algebraic[114] = (algebraic[111]-algebraic[65])--algebraic[79] algebraic[33] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[53]*constants[2]*algebraic[0] algebraic[76] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[56]*constants[2]*algebraic[0] algebraic[77] = constants[36]*(exp(algebraic[33]/(constants[0]*constants[1]))-exp(algebraic[76]/(constants[0]*constants[1]))) algebraic[46] = constants[0]*constants[1]*log(constants[14]*states[7])+constants[54]*constants[2]*algebraic[0] algebraic[112] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[57]*constants[2]*algebraic[0] algebraic[115] = constants[49]*(exp(algebraic[46]/(constants[0]*constants[1]))-exp(algebraic[112]/(constants[0]*constants[1]))) algebraic[117] = ((algebraic[101]-algebraic[115])+algebraic[59])-algebraic[77] algebraic[35] = constants[0]*constants[1]*log(constants[18]*states[9])+constants[53]*constants[2]*algebraic[0] algebraic[80] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[56]*constants[2]*algebraic[0] algebraic[81] = constants[38]*(exp(algebraic[35]/(constants[0]*constants[1]))-exp(algebraic[80]/(constants[0]*constants[1]))) algebraic[118] = (algebraic[115]+algebraic[63])-algebraic[81] algebraic[48] = constants[51]*constants[2]*algebraic[0] algebraic[17] = constants[0]*constants[1]*log(constants[21]*states[1]) algebraic[18] = constants[0]*constants[1]*log(constants[19]*states[17]) algebraic[19] = algebraic[17]+constants[51]*constants[2]*algebraic[0]+algebraic[18] algebraic[16] = constants[0]*constants[1]*log(constants[20]*states[0]) algebraic[49] = algebraic[16]+algebraic[18] algebraic[120] = custom_piecewise([equal(algebraic[48] , 0.00000), constants[22]*(exp(algebraic[19]/(constants[0]*constants[1]))-exp(algebraic[49]/(constants[0]*constants[1]))) , True, (((constants[22]*algebraic[48])/(constants[0]*constants[1]))/(exp(algebraic[48]/(constants[0]*constants[1]))-1.00000))*(exp(algebraic[19]/(constants[0]*constants[1]))-exp(algebraic[49]/(constants[0]*constants[1])))]) algebraic[47] = constants[0]*constants[1]*log(constants[15]*states[15])+constants[54]*constants[2]*algebraic[0] algebraic[116] = constants[0]*constants[1]*log(constants[19]*states[17])+constants[57]*constants[2]*algebraic[0] algebraic[119] = constants[50]*(exp(algebraic[47]/(constants[0]*constants[1]))-exp(algebraic[116]/(constants[0]*constants[1]))) algebraic[121] = ((algebraic[104]-algebraic[119])+algebraic[61])--algebraic[77] algebraic[122] = constants[2]*((((((((((((((((((((((((((((((((((((((((((((((((((((((((-constants[51]*algebraic[120]+constants[57]*algebraic[83])-constants[54]*algebraic[83])+constants[57]*algebraic[95])-constants[54]*algebraic[95])+constants[57]*algebraic[89])-constants[54]*algebraic[89])+constants[57]*algebraic[101])-constants[54]*algebraic[101])+constants[57]*algebraic[115])-constants[54]*algebraic[115])+constants[55]*algebraic[51])-constants[52]*algebraic[51])+constants[55]*algebraic[59])-constants[52]*algebraic[59])+constants[56]*algebraic[71])-constants[53]*algebraic[71])+constants[56]*algebraic[75])-constants[53]*algebraic[75])+constants[56]*algebraic[79])-constants[53]*algebraic[79])+constants[57]*algebraic[107])-constants[54]*algebraic[107])+constants[55]*algebraic[55])-constants[52]*algebraic[55])+constants[55]*algebraic[63])-constants[52]*algebraic[63])+constants[56]*algebraic[67])-constants[53]*algebraic[67])+constants[57]*algebraic[86])-constants[54]*algebraic[86])+constants[57]*algebraic[98])-constants[54]*algebraic[98])+constants[57]*algebraic[92])-constants[54]*algebraic[92])+constants[57]*algebraic[104])-constants[54]*algebraic[104])+constants[57]*algebraic[119])-constants[54]*algebraic[119])+constants[55]*algebraic[53])-constants[52]*algebraic[53])+constants[55]*algebraic[61])-constants[52]*algebraic[61])+constants[56]*algebraic[69])-constants[53]*algebraic[69])+constants[56]*algebraic[73])-constants[53]*algebraic[73])+constants[56]*algebraic[77])-constants[53]*algebraic[77])+constants[56]*algebraic[81])-constants[53]*algebraic[81])+constants[57]*algebraic[111])-constants[54]*algebraic[111])+constants[55]*algebraic[57])-constants[52]*algebraic[57])+constants[55]*algebraic[65])-constants[52]*algebraic[65]) algebraic[123] = (algebraic[119]+algebraic[65])-((-algebraic[81]+algebraic[120])-algebraic[120]) algebraic[1] = constants[0]*constants[1]*log(constants[4]*states[2]) algebraic[2] = constants[0]*constants[1]*log(constants[6]*states[3]) algebraic[3] = constants[0]*constants[1]*log(constants[8]*states[4]) algebraic[4] = constants[0]*constants[1]*log(constants[10]*states[5]) algebraic[5] = constants[0]*constants[1]*log(constants[12]*states[6]) algebraic[6] = constants[0]*constants[1]*log(constants[14]*states[7]) algebraic[7] = constants[0]*constants[1]*log(constants[16]*states[8]) algebraic[8] = constants[0]*constants[1]*log(constants[18]*states[9]) algebraic[9] = constants[0]*constants[1]*log(constants[5]*states[10]) algebraic[10] = constants[0]*constants[1]*log(constants[7]*states[11]) algebraic[11] = constants[0]*constants[1]*log(constants[9]*states[12]) algebraic[12] = constants[0]*constants[1]*log(constants[11]*states[13]) algebraic[13] = constants[0]*constants[1]*log(constants[13]*states[14]) algebraic[14] = constants[0]*constants[1]*log(constants[15]*states[15]) algebraic[15] = constants[0]*constants[1]*log(constants[17]*states[16]) return algebraic def custom_piecewise(cases): """Compute result of a piecewise function""" return select(cases[0::2],cases[1::2]) 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)