On the Numerical Solution of Involutive Ordinary Differential Systems: Higher Order Methods