“…where f ðxÞ is the unknown solution, T D ðα r ,β r Þ x is tempered fractional derivatives, q r ∈ ℝ, r = 1, ⋯, l are constants, and α r , β r ≥ 0 are real derivative orders which β denotes the tempered coefficient, while hðxÞ is the unhomogeneous terms. To date, various numerical or analytical methods were derived to find the solution for different fractional calculus problems, such as [11][12][13]. On top of that, the operational matrix method via different types of the polynomial is one of the common numerical schemes which had been widely used in solving various types of fractional calculus problems, such as the poly-Bernoulli operational matrix for solving fractional delay differential equation [14], poly-Genocchi operational matrix for solving fractional differential equation [15], Jacobi wavelet operational matrix of fractional integration for solving fractional integro-differential equation [16], and Fibonacci wavelet operational matrix of integration for solving of nonlinear Stratonovich Volterra integral equations [17].…”