Numerical Construction of Taylor Series Approximations for a Set of Simultaneous First Order Differential Equations