% A = crsdiag(M)
% Provide a cross diagnal matrix. The diagnal
% elements of the matrix will be 1,-1,1,-1 ...

function A = crsdiag(M)
for I = 1:M
  a(I) = (-1)^(I-1);
end
A = diag(a);