A discussion on the approximate solutions of first order systems of non-linear ordinary equations

We develop a one step matrix method in order to obtain approximate solutions of first order systems and non-linear ordinary differential equations, reducible to first order systems. We find a sequence of such solutions that converge to the exact solution. We study the precision, in terms of the local error, of the method by applying it to different well known examples. The advantage of the method over others widely used lies on the simplicity of its implementation.