Qualitative properties of numerical solutions of PDE models of disease propagation