function diag = domdiag(mat) [f,c]=size(mat); diag=1; for i=1:f sum=0; for j=1:c sum= sum+ abs(mat(i,j))-abs(mat(i,j)); if(abs(mat(i,j))<=sum) diag=0; return end end end if (diag==1) return end