Numerical methods for the solution of large kinetic systems