tra51_1

[n m]=size(a);
z=n*m;
p=num2cell(1:z)
for i=1:z
    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}'