k1 = r k2 = K k3 = d k4 = m k5 = beta k6 = k k7 = a k8 = l1 k9 = alpha k10 = b k11 = u k12 = l2 k13 = I_plot k14 = T_plot k15 = tme