This paper focuses on presenting an accurate, efficient, and fast pseudospectral method to solve tempered fractional differential equations (TFDEs) in both spatial and temporal dimensions. We use a tempered integration matrix that can be employed to evaluate ๐-fold tempered integrals of a real function ๐ for ๐ โ โ + . Also, it may be used to compute any non-integer order ๐ < 0 tempered derivatives of ๐. We employ the Chebyshev interpolating polynomial for ๐ at Gauss-Lobatto (GL) points in the range [โ1, 1] and any identically shifted range. The proposed method carries with it a recast of the TFDE into integration formulas to actually take advantage of the adaptation of the integral operators, hence avoiding the ill-conditioning and reduction in the convergence rate of integer differential operators. Via various tempered fractional differential examples, the present technique shows many advantages, for instance, high accuracy, a much higher rate of running, fewer computational hurdles and programming, calculating the tempered-derivative/integral of fractional order and its exceptional accuracy in comparison with other competitive numerical schemes. The study includes the elapsed times taken to construct the collocation matrices and obtain the numerical solutions. Also, numerical examination of the produced condition number ๐
(๐ด) of the resulted linear systems. The accuracy and efficiency of the proposed method are studied from the standpoint of the ๐ฟ 2 and ๐ฟ โ -norms error and fast rate of spectral convergence.