A = [0 1 0 0 1; 0 0 1 1 0; 1 0 0 0 0; 1 0 0 0 1; 0 1 0 1 0];

M=A+A^2 = [0 2 1 2 1; 2 0 1 1 1; 1 1 0 0 1; 1 2 0 1 2; 1 1 1 2 1]

go = true;
while (go == true)
    for i = 1:5
        for j = 1:5
            if( A(i,j) == 0)
                A(i,j) = 1;
                M = A+A^2;
                if ((numel(M)-nnz(M)) == 0)
                    disp(M)
                    go = false;
                end
            end
        end
    end
end

while (1)
    for i = 1:5
        for j = 1:5
            if( A(i,j) == 0)
                A(i,j) = 1;
                M = A+A^2;
                if ((numel(M)-nnz(M)) == 0)
                    disp(M)
                    go = false;
                end
            end
        end
    end
end