We applied a previously established technique to calculate R c [1] , and considered infectious classes in the order X = (A, I U , I T , A r , I r,U , I r,T , I T ,r ) to simplify computations. The system (2)-(4) can be written as