ClearAll["Global`*"]
(*R = reds, G = greys*)
(*S = susceptible, I = infected, R = recovered*)
tseries = {t, 0, 3};
vars = {HG[t], HR[t], SG[t], IG[t], RG[t], SR[t], IR[t], aG[t], qG[t]};
b = 0.4; (*natural mortality rate, both species*)
\[Beta] = 0.7; (*rate of virus transmission*)
aR =...