Designing a Matrix Collocation Method for Fractional Delay Integro-Differential Equations with Weakly Singular Kernels Based on Vieta–Fibonacci Polynomials

In the present work, the numerical solution of fractional delay integro-differential equations (FDIDEs) with weakly singular kernels is addressed by designing a Vieta–Fibonacci collocation method. These equations play immense roles in scientific fields, such as astrophysics, economy, control, biology, and electro-dynamics. The emerged fractional derivative is in the Caputo sense. By resultant operational matrices related to the Vieta–Fibonacci polynomials (VFPs) for the first time accompanied by the collocation method, the problem taken into consideration is converted into a system of algebraic equations, the solving of which leads to an approximate solution to the main problem. The existence and uniqueness of the solution of this category of fractional delay singular integro-differential equations (FDSIDEs) are investigated and proved using Krasnoselskii’s fixed-point theorem. A new formula for extracting the VFPs and their derivatives is given, and the orthogonality of the derivatives of VFPs is easily proved via it. An error bound of the residual function is estimated in a Vieta–Fibonacci-weighted Sobolev space, which shows that by properly choosing the number of terms of the series solution, the approximation error tends to zero. Ultimately, the designed algorithm is examined on four FDIDEs, whose results display the simple implementation and accuracy of the proposed scheme, compared to ones obtained from previous methods. Furthermore, the orthogonality of the VFPs leads to having sparse operational matrices, which makes the execution of the presented method easy.