Computing Moore-Penrose inverses of Toeplitz matrices by Newtonʹs iteration

We modify the algorithm of [1], based on Newtonʹs iteration and on the concept of ε-displacement rank, to the computation of the Moore-Penrose inverse of a rank-deficient Toeplitz matrix. Numerical results are presented to illustrate the effectiveness of the method.