In this paper, a numerical method for solving is the non-holonomic system is given. The non-holonomic constrains is satisfied within the time step and the solutions of the system on the discrete time are transformed by using the impulses (Lagrange multiplies). The numerical method is of second order precision and it can satisfy precisely the non-holonomic constraints and gives good performance for long time integration.