“…These include the variational iteration method [6], automatic differentiation method [7], finite difference schemes [8], modified cubic B-splines collocation and cubic hyperbolic Bspline method [9,10], least-squares quadratic B-spline finite element method [11], finite element collocation method [12], mixed finite volume element methods [13], sinc-Galerkin method [14], and Laplace transform decomposition method [15]. Orthogonal polynomials play a vital role in solving a wide range of mathematical problems, including differential and integro-differential equations [16][17][18][19][20][21][22][23][24][25]. Chelyshkov polynomials, a class of orthogonal polynomials, possess properties similar to classical orthogonal polynomials and are associated with hypergeometric functions, orthogonal exponential polynomials, and Jacobi polynomials [26,27].…”