import math import numpy as np def kinetic_parameters(M, include_type2_reactions, dims, V): fast_kinetic_constant = 1e6 Km_TRPN = 0.5e-3 # unit mM kf_TRPN = fast_kinetic_constant kr_TRPN = fast_kinetic_constant * Km_TRPN W = [1, V['V_myo'], V['V_myo'], V['V_myo']] k_kinetic = [kf_TRPN, kr_TRPN] N_cT = [] K_C = [] return (k_kinetic, N_cT, K_C, W)