nant [29,30,31]. (1) dSdt=μN−S−βAA+βIIN+βWWS (2) dEdt=βAA+βIIN+βWWS−ε+μE (3) dAdt=εE−ω+μA (4) dIdt=1−pωA−v+μII (5) dRdt=pωA+vI−μR (6) dWdt=σAA+σIIN