function flag=identitymat(A)
sz=size(A);
n=sz(1);

flag=22;
for i=1:1:n
    for j=1:1:n
        if (A(i,j)~=1) & (i==j)
            flag=66;
            break;
        end
       if (A(i,j)~=0) & (i~=j)
           flag=66;
           break;
       end
   end
end