Model Mathematics

\frac{d}{d time} (T) = s + lambda T (1 - (\frac{T + T1_star + T2_star}{Tmax})) - (mu T) - (k1 (1 - r1) T V1_I) - (k2 (1 - r2) T V2_I)

\frac{d}{d time} (M) = lambda_M - (mu_M M) - (k1_M (1 - r1) M V1_I) - (k2_M (1 - r2) M V2_I)

\frac{d}{d time} (T1_star) = m11 k1 (1 - r1) T V1_I + m21 k2 (1 - r2) T V2_I - (delta1 T1_star) \frac{d}{d time} (T2_star) = m22 k2 (1 - r2) T V2_I + m12 k1 (1 - r1) T V1_I - (delta2 T2_star)

\frac{d}{d time} (M1_star) = m11 k1_M (1 - r1) M V1_I + m21 k2_M (1 - r2) M V2_I - (delta1_M M1_star) \frac{d}{d time} (M2_star) = m22 k2_M (1 - r2) M V2_I + m12 k1_M (1 - r1) M V1_I - (delta2_M M2_star)

\frac{d}{d time} (V1_I) = (1 - p1) N1 delta1 T1_star + (1 - p1) N1_M delta1_M M1_star - (c1 V1_I) \frac{d}{d time} (V1) = N1 delta1 T1_star + N1_M delta1_M M1_star - (c1 V1)

\frac{d}{d time} (V2_I) = (1 - p2) N2 delta2 T2_star + (1 - p2) N2_M delta2_M M2_star - (c2 V2_I) \frac{d}{d time} (V2) = N2 delta2 T2_star + N2_M delta2_M M2_star - (c2 V2)