Abstract
Recently, direct methods that involve higher derivatives to numerically approximate higher order initial value problems (IVPs) have been explored, which aim to construct numerical methods with higher order and very high precision of the solutions. This article aims to construct a fourth and fifth derivative, three-point implicit block method to tackle general third-order ordinary differential equations directly. As a consequence of the increase in order acquired via the implicit block method of higher derivatives, a significant improvement in efficiency has been observed. The new method is derived in a block mode to simultaneously evaluate the approximations at three points. The derivation of the new method can be easily implemented. We established the proposed method's characteristics, including order, zero-stability, and convergence. Numerical experiments are used to confirm the superiority of the method. Applications to problems in physics and engineering are given to assess the significance of the method.