The Toeplitz matrix
Tn+1 is not positive definite to working accuracy. If, on exit,
XN is close to unity, then the principal minor was probably close to being singular, and the sequence
τ0,τ1,…,τN may be a valid sequence nevertheless.
X returns the solution of the equations
and
V returns
vn, but it may not be positive.