  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.
