Numerical simulation of a kinetic model for chemotaxis