Abstract
Volterra integro-partial differential equations with weakly singular kernels (VIPDEWSK) are utilized to model diverse physical phenomena. A matrix collocation method is proposed for determining the approximate solution of this functional equation category. The method employs shifted Chebyshev polynomials of the fifth kind (SCPFK) to construct two-dimensional pseudo-operational matrices of integration, avoiding the need for explicit integration and thereby speeding up computations. Error bounds are examined in a Chebyshev-weighted space, providing insights into approximation accuracy. The approach is applied to several experimental examples, and the results are compared with those obtained using the Bernoulli wavelets and Legendre wavelets methods.