The model is solved numerically using finite difference schemes, in particular, central difference schemes are applied to diffusion and chemotaxis terms and a forward Euler method is used for time stepping.