%%?IEEE-33節點系統潮流計算
%輸入的支路信息是有名值
%輸入Bus矩陣，第一列節點號，第二列首節點/有功，第三列尾節點/無功
%輸入Branch矩陣第一列存支路號，第2列存支路首節點，
%%%%%%%%%%第3列存支路尾節點，第4存支路電阻，第5列存支路電抗）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%format?long;
clc
clear?all
tic
h=waitbar（0‘please?wait‘）;

for?i=1:1000
????%computation?here%
????str=[‘運行中...‘num2str（i/1000*100）‘%‘];
????waitbar（i/1000hstr）
end
delete（h）;
Bus=[100;210060;39040;412080;56030;66020;7200100;8200100;96020;106020;114530;126035;136035;1412080;156010;166020;176020;189040;199040;209040;219040;229040;239040;24420200;25420200;266025;276025;286020;2912070;30200600;3115070;32210100;336040;];
Branch=[1120.09220.0407;2230.49300.2511;3340.36600.1864;4450.38110.1941;5560.81900.7070;6670.18720.6188;7780.71440.2351;8891.03000.7400;99101.04400.7400;1010110.19660.065;1111120.37440.1238;1212131.46801.1550;1313140.54160.7129;1414150.59100.5260;1515160.74630.5450;1616171.28901.7210;1717180.73200.5740;182190.16400.1565;1919201.50421.3554;2020210.40950.4784;2121220.70890.9373;223230.45120.3083;2323240.89800.7091;2424250.89600.7011;256260.20300.1034;2626270.28420.1447;2727281.05900.9337;2828290.80420.7006;2929300.50750.2585;3030310.97440.9630;3131320.31050.3619;3232330.34100.5302;];
[busnumrow]=size（Bus）;
[branchnumrow]=size（Branch）;
soubus=Branch（:2）;
mobus=Branch（:3）;
Vbus=ones（busnum1）;
Vbus（:1）=12.66;
Vbus1=Vbus;
Ploss=zeros（busnum1）;
Qloss=zeros（busnum1）;
e=1;
k=0;
Branch1=Branch;
n=1;
%%%%%%%%%%%%%%%支路重新排序，各個分支線同時進行計算
while?~isempty（Branch1）%%%%T