Improving the Efficiency of Runge-Kutta Methods for the Solution of BVPs for Higher-Order ODEs