Abstract
This paper proposes a new pair of block techniques (NPBT) for the direct solution of third-order singular initial-value problems (IVPs). The proposed method uses a polynomial and two intermediate points to approximate the theoretical solution of third-order singular IVPs, resulting in a reasonable approximation within the integration interval. The method's essential features, including stability and convergence order, are analyzed. The proposed NPBT method is improved by using an embedding-like strategy that allows it to be executed in a variable step size mode in order to gain better efficiency. The effectiveness of the proposed method is assessed using various model problems. The approximate solution provided by the proposed NPBT method is more accurate than that of the existing methods utilized for comparison. This efficient solution positions NPBT as a good numerical method for integrating third-order singular IVP models in the fields of applied sciences and engineering.