In this study, we investigate a class of mixed fractional partial integro-differential equations (FrPI-DE) involving symmetric singular kernels. The considered model problem involves Caputo fractional derivatives and integral operators that describe spatial interactions in a bounded domain. For the purpose of analysis, the original problem is reformulated in the form of a nonlinear Volterra–Fredholm integral equation (NV-FIE). The existence and uniqueness of the solution are established by the Banach fixed point theorem. To compute numerical solutions, a modified Toeplitz matrix method (TMM) is proposed to handle the singular kernel efficiently. The method transforms the integral equation to a system of nonlinear algebraic equations, which can be solved numerically. The convergence properties of the resulting numerical scheme are analyzed and illustrate the effectiveness of the method by providing numerical examples involving logarithmic, Cauchy-type, and weakly singular kernels. Numerical results indicate that the proposed method provides highly accurate approximations and exhibits stable convergence behavior for different parameter values. Furthermore, these results confirm the effectiveness and reliability of the proposed method for solving fractional integro-differential equations that include symmetric singular kernels.
Nasr et al. (Sat,) studied this question.