Falkner hybrid block methods for second-order IVPs: A novel approach to enhancing accuracy and stability properties

Second-order initial value problems (IVPs) in ordinary differential equations (ODEs) are ubiquitous in various fields, including physics, engineering, and economics. However, their numerical integration poses significant challenges, particularly when dealing with oscillatory or stiff problems. This article introduces a novel Falkner hybrid block method for the numerical integration of second-order IVPs in ODEs. The newly developed method is of order six with a large interval of absolute stability and is implemented using a fixed step size technique. The numerical experiments show the accuracy of our methods when compared with Falkner linear multistep methods, block methods, and other hybrid codes proposed in the scientific literature. This innovative approach demonstrates improved accuracy and stability in solving second-order IVPs, making it a valuable tool for researchers and practitioners.