where f = β(δ A A + I U + δ T I T ), g = δ r β(δ A A r + I r,U + I r,T ), and the prime ' ′ ' denotes the derivative of the numbers in compartments with respect to the time; β is the baseline transmission rate of the sensitive strain; p is the probability of developing clinical symptoms;