function V=powerflow(A0)
A0T = A0';                              % 形成 A0 的轉置
S = [0;-NodeData(:,2) - i*NodeData(:,3)];       % 形成 S
ZL = [0;BranchData(:,3) + i*BranchData(:,4)];       % 形成 ZL

V = ones(NN,1);
V(1) = 1.05;                            % 各個節點電壓賦初值

IL(NN,1) = -conj(S(NN) / V(NN));          % 最末支路電流賦初值

Delta = 1;                              % 收斂判據賦初值
TempV = V;                              % 賦初值，用于記憶上次迭代結果
while Delta > 1e-8
    IN = conj(S ./ V);                    % 節點注入電流
    
    for n = 1:NN-1
        IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n);
    end                                 % 電流回代過程
    
    for n = 2:NN
        V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n);
    end                                 % 電壓前推過程
    
    Delta = max(abs(V-TempV));          % 更新收斂判據
    TempV = V;                          % 記憶迭代結果
end