On the numerical solution of involutive ordinary differential systems: enhanced linear algebra