function chksd(A);
lam = abs(eig(A));
A   = A - diag(diag(A));
sm  = sum(A')';
if any(lam<sm); 
  disp('not strictly diagonally dominant')
else
  disp('matrix is strictly diagonally dominant')
end
