“…The intimate relation between Chebyshev polynomials and Fibonacci, Lucas, and more generally Horadam's numbers (i. e., solutions to a general second-order linear recursion with constant coefficients) has been intensively explored, among other areas, in enumerative combinatorics for computing generating functions (see, for instance, [20,27,30]), number theory for representation of integer sequences (see, for instance, an early work [6] and recent [33,35]), and applied linear algebra for computing determinants and inverse matrices (see, for instance, [13,14,17,40]). The formula (3) also provides a natural link between the theory of Chebyshev polynomials and a discrete Poisson equation and its Sturm-Liouville theory [7,[13][14][15]22], and subsequently to the harmonic functions and first passage time of random walks (see, for instance, the classical [24,29] and more recent [10,12,18]). Altogether, the link fundamentally bridges between analytic and probabilistic methods exploited in this paper.…”