Authors’ summary: We propose a method for the numerical solution of ordinary differential systems. The system is considered geometrically as a submanifold in a jet space. The solutions are then certain integral manifolds that can be computed numerically when the system has been transformed into involutive form.