“…A matrix spectral method was employed to approximate solutions for variable-order time-fractional weakly singular integro-partial differential equations (IPDEs) by [15] . Authors have introduced and applied successfully different operational schemes based on the Legendre polynomials to study Volterra singular partial integro-differential equations [16] , [17] , [18] . Yang et al utilized an implicit difference scheme to handle fourth-order weakly singular integro-partial differential equations (WSIPDEs) [19] .…”