lab

num = 0;
W = [0,7,num,3,4,num,7,num;
    5,0,num,8,num,5,num,8;
    num,3,0,6,num,4,num,7;
    5,8,6,0,9,7,num,num;
    5,8,7,6,0,5,num,3;
    3,num,num,8,num,0,8,4;
    num,7,num,num,num,6,0,9;
    7,num,6,num,num,num,4,0];
A = [0,1,0,1,1,0,1,0;
    1,0,0,1,0,1,0,1;
    0,1,0,1,0,1,0,1;
    1,1,1,0,1,1,0,0;
    1,1,1,1,0,1,0,1;
    1,0,0,1,0,0,1,1;
    0,1,0,0,0,1,0,1;
    1,0,1,0,0,0,1,0];
names = {'A','B','C','D','E','F','G','H',};
G = digraph(W, names);
p = plot(G,'EdgeLabel', G.Edges.Weight);
num = 10000;
W = [0,7,num,3,4,num,7,num;
    5,0,num,8,num,5,num,8;
    num,3,0,6,num,4,num,7;
    5,8,6,0,9,7,num,num;
    5,8,7,6,0,5,num,3;
    3,num,num,8,num,0,8,4;
    num,7,num,num,num,6,0,9;
    7,num,6,num,num,num,4,0];
W1 = W;
n = size(W,1);
for k = 1:n
    for i = 1:n
        for j = 1:n
            W1(i,j) = min(W1(i,j), W(i,k) + W(k,j));
        end
    end
end


sim_A1 = triu(A) + transpose(triu(A,1));
sim_A2 = triu(transpose(A)) + transpose(triu(transpose(A),1));
A3 = sim_A2+sim_A1;
sim_A = A3;
for i = 1:n
    for j = 1:n
        if A3(i,j) ~= 0
            sim_A(i,j) = 1;
        end
    end
end

D = A3;
sum_line = 0;

for i = 1:n
    for j = 1:n
        if i == j
            for k = 1:n
                sum_line = sum_line + sim_A(i,k);
            end
            D(i,j) = sum_line;
            sum_line = 0;
        else
            D(i,j) = 0;
        end
    end
end



% D = [6,0,0,0,0,0,0,0;
%     0,7,0,0,0,0,0,0;
%     0,0,5,0,0,0,0,0;
%     0,0,0,5,0,0,0,0;
%     0,0,0,0,6,0,0,0;
%     0,0,0,0,0,7,0,0;
%     0,0,0,0,0,0,4,0;
%     0,0,0,0,0,0,0,6];
K = D - sim_A;
N = ((-1)^16)*det(K(1:7,1:7));