“…In practice, the coefficients r * ij of the orthonormal polynomial (17) can be derived directly from the elements r ij of the matrix R, avoiding the inversion of the matrix R. The entries r ij are used to compute the coefficients a j and b j of the three-term recurrence relation (18) for 1 ≤ j ≤ p, where ψ −1 (ξ) = 0 and ψ 0 (ξ) = 1. More specifically, the coefficients a j and b j are written in terms of the elements r ij as a j = r j,j+1 r j,j − r j−1,j r j−1,j−1 and b j = r j+1,j+1 r j,j , for 1 ≤ j ≤ p, with r 0,0 = 1 and r 0,1 = 0.…”