The focus of this investigation centers on the formulation of a novel spectral method deployed to numerically solve partial integro-differential equations that possess memory features. The approach leverages the hexic shifted Chebyshev polynomials, and addresses the nonlinearity of the computational output using iterative methods. A comprehensive explanation of the methodology is presented, and a thorough convergence analysis is conducted. This technique has the capacity to be effortlessly modified and used to solve a plethora of linear and nonlinear issues, while minimizing computational time. Ultimately, the effectiveness of this novel strategy is demonstrated through the successful resolution of two exemplary problems. The findings of this study suggest that this spectral approach holds significant promise for solving partial integro-differential equations.
