  A quasi-Toeplitz (QT) matrix is a semi-infinite matrix of the form $A=T(a)+E$
where $T(a)$ is the Toeplitz matrix with entries $(T(a))_{i,j}=a_{j-i}$, for
$a_{j-i}\in\mathbb C$, $i,j\ge 1$, while $E$ is a matrix representing a compact
operator in $\ell^2$. The matrix $A$ is finitely representable if $a_k=0$ for
$k<-m$ and for $k>n$, given $m,n>0$, and if $E$ has a finite number of nonzero
entries. The problem of numerically computing eigenpairs of a finitely
representable QT matrix is investigated, i.e., pairs $(\lambda,{\bf v})$ such
that $A{\bf v}=\lambda {\bf v}$, with $\lambda\in\mathbb C$, ${\bf
v}=(v_j)_{j\in\mathbb Z^+}$, ${\bf v}\ne 0$, and $\sum_{j=1}^\infty
|v_j|^2<\infty$. It is shown that the problem is reduced to a finite nonlinear
eigenvalue problem of the kind $ WU(\lambda){\pmb \beta}=0$, where $W$ is a
constant matrix and $U$ depends on $\lambda$ and can be given in terms of
either a Vandermonde matrix or a companion matrix. Algorithms relying on
Newton's method applied to the equation $\det WU(\lambda)=0$ are analyzed.
Numerical experiments show the effectiveness of this approach. The algorithms
have been included in the CQT-Toolbox [Numer. Algorithms 81 (2019), no. 2,
741--769].
