import math import numpy as np def kinetic_parameters(M, include_type2_reactions, dims, V): fast_kinetic_constant = 1e9 num_cols = dims['num_cols'] num_rows = dims['num_rows'] pK1 = 6.03 Km_B1 = pow(10,-pK1) # unit mM kf_B1 = fast_kinetic_constant kr_B1 = fast_kinetic_constant * Km_B1 W = [1]+ [V['V_myo']]*num_rows k_kinetic = [kf_B1, kr_B1] N_cT = [] K_C = [] return (k_kinetic, N_cT, K_C, W)