Model equations:dSdt=−β(t)⋅S(I1+I2)dE1dt=β(t)⋅S(I1+I2)−2σE1dE2dt=2σE1−2σE2dI1dt=2σE2−2γI1dI2dt=2γI1−2γI2dRdt=2γI2