[n m]=size(a);
z=n*m;
for i=1:z
p{i}=i;
end
for i=1:99
for k=1:z
x = p{k};
e=x(end);
h = a(e)+1;
g = e-1;
w = e+n;
f = mod(g,n);
if w<=z & a(w)==h
p{k} = [x w];
elseif e>n & a(e-n)==h
p{k} = [x e-n];
elseif f<n-1 & a(e+1)==h
p{k} = [x e+1];
elseif f & a(g)==h
p{k} = [x g];
end
d(k) = nnz(p{k});
end
end
[v b]=max(d);
b = p{b}'