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