In this paper, we compare the performances of two Butcher-based block hybrid methods for the numerical integration of initial value problems. We compare the condition numbers of the linear system of equations arising from both methods and the absolute errors of the solution obtained. The results of the numerical experiments illustrate that the better conditioned method outperformed its less conditioned counterpart based on the absolute errors. In addition, after applying our method on some examples, it was discovered that the absolute errors in this work were better than those of a recent study in the literature. Hence, we recommend this method for the numerical solution of stiff and non-stiff initial value problems.
