Existence, Stability, and Numerical Methods for Multi-Fractional Integro-Differential Equations with Singular Kernel

This work investigates the solutions of fractional integro-differential equations (FIDEs) using a unique kernel operator within the Caputo framework. The problem is addressed using both analytical and numerical techniques. First, the two-step Adomian decomposition method (TSADM) is applied to obtain an exact solution (if it exists). In the second part, numerical methods are used to generate approximate solutions, complementing the analytical approach based on the Adomian decomposition method (ADM), which is further extended using the Sumudu and Shehu transform techniques in cases where TSADM fails to yield an exact solution. Additionally, we establish the existence and uniqueness of the solution via fixed-point theorems. Furthermore, the Ulam–Hyers stability of the solution is analyzed. A detailed error analysis is performed to assess the precision and performance of the developed approaches. The results are demonstrated through validated examples, supported by comparative graphs and detailed error norm tables (L∞, L2, and L1). The graphical and tabular comparisons indicate that the Sumudu-Adomian decomposition method (Sumudu-ADM) and the Shehu-Adomian decomposition method (Shehu-ADM) approaches provide highly accurate approximations, with Shehu-ADM often delivering enhanced performance due to its weighted formulation. The suggested approach is simple and effective, often producing accurate estimates in a few iterations. Compared to conventional numerical and analytical techniques, the presented methods are computationally less intensive and more adaptable to a broad class of fractional-order differential equations encountered in scientific applications. The adopted methods offer high accuracy, low computational cost, and strong adaptability, with potential for extension to variable-order fractional models. They are suitable for a wide range of complex systems exhibiting evolving memory behavior.