In this article, we used a novel semi-analytical approach, named the optimal auxiliary function method (OAFM), to solve integro-differential equations (IDEs). The OAFM includes an auxiliary function and convergence control parameters, which expedite the convergence of the method. To apply the proposed method, some assumptions regarding small or large parameters in the problem are necessary. We present numerical outcomes acquired via the OAFM alongside those obtained from other numerical techniques in tables. Furthermore, we demonstrate the efficacy and ease of implementing the proposed method for various IDEs using 2D graphs.
