[n m]=size(a);
z=n*m;
for i=1:z
p{i}=i;
end
for i=';;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;'
for k=1:z
e=p{k}(end);
h = a(e)+1;
g = e-1;
w = e+n;
x = p{k};
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>0 & a(g)==h
p{k} = [x; g];
end
d(k) = nnz(p{k});
end
end
[v b]=max(d);
b = p{b}