NUMERICAL INTEGRATION OF LARGE KINETIC SYSTEMS