We present and analyse numerical methods with operator splitting procedures, applied to an epidemic model which takes into account the space-dependence of infection. derive conditions on time step, under preserve non-negativity monotonicity properties exact solution. Our results are illustrated by experiments.