function c=test(phi, n)
x=inv(toeplitz(phi.^(0:(n-1))));
sum(diag(x))