b=0 [m n]=size(a); for k=1:9 for i=1:m-k for j=1:n-k c=a(i:i+k,j:j+k); b=max(b,sum(all(c(1)==c(:))*c(:))); for i=1:m-k for j=1:n-k c=a(i:i+k,j:j+k); b=max(b,sum(all(c==c(1))*c(:))); end end end end end