We introduce a numerical scheme to approximate a quasilinear hyperbolic system which models the movement of cells under the influence of chemotaxis. Since we expect to find solutions which contain vacuum parts, we propose an upwinding scheme which properly handles the presence of vacuum and which gives a good approximation of the time asymptotic states of the system. For this scheme we prove some basic analytical properties and study its stability near some of the steady states of the system.