  For a family of near banded Toeplitz matrices, generalized characteristic
polynomials are shown to be orthogonal polynomials of two variables, which
include the Chebyshev polynomials of the second kind on the deltoid as a
special case. These orthogonal polynomials possess maximal number of real
common zeros, which generate a family of Gaussian cubature rules in two
variables.
