Abstract
In this paper we construct an algorithm to approximate the solution of the initial value problem y (t) = f (t, y) with y(t0) = y0. The method is implicit and combines the classical Simpson’s rule with the Simpson’s
3/8 rule to yield an unconditionally A-stable method of order 4.