function ShowCS()
% Illustrates the orthogonality that exists among the columns of
% P = CSmatrix(n)

clc
format short
disp('If P = CSmatrix(n), then P''*P is diagonal...')
P = CSmatrix(10)
disp('P''P = ')
disp(P'*P)
disp('D = P''*P = (n/2)I except that D(1,1) = D(m+1,m+1) = n...')
n = 10;
P = CSmatrix(n);
d = diag(P'*P);
fprintf('\n')
for k=1:n
fprintf('D(%2d,%2d) = %2d\n',k,k,round(d(k)))
end