import math import numpy as np def kinetic_parameters(M, include_type2_reactions, dims, V): fast_kinetic_constant = 1e6 num_cols = dims['num_cols'] num_rows = dims['num_rows'] pK2 = 7.57 Km_B2 = pow(10,-pK2) # unit mM kf_B2 = fast_kinetic_constant kr_B2 = fast_kinetic_constant * Km_B2 W = [1]+ [V['V_myo']]*num_rows k_kinetic = [kf_B2, kr_B2] N_cT = [] K_C = [] return (k_kinetic, N_cT, K_C, W)