Exponentially Fitted Multi-Order Block Methods for Solving Initial Value Problems

This study presents a novel exponentially fitted multi-order block numerical method designed to efficiently solve first- and second-order initial value problems (IVPs) arising in ordinary differential equations. Leveraging interpolation and collocation techniques with Chebyshev polynomials as basis functions, the proposed method incorporates exponential terms to enhance precision and performance.