clear
n=33;%n=input（‘請輸入節(jié)點數(shù)‘）;
b=32;%b=input（‘請輸入支路數(shù)‘）;
%disp（‘請輸入支路阻抗和節(jié)點功率矩陣，第一列存支路號，第二列存首節(jié)點號，第三列存尾節(jié)點號，第四列存支路阻抗，第五列存節(jié)點給定功率’）;
Sb=10;
Ub=12.66;
Zb=Ub^2/Sb;
Z=[?1?0?10.0922/Zb+0.0470i/Zb0.100/Sb+0.060i/Sb
????2?1?20.4930/Zb+0.2511i/Zb0.090/Sb+0.040i/Sb
????3?2?30.3660/Zb+0.1864i/Zb0.120/Sb+0.080i/Sb
????4?3?40.3811/Zb+0.1941i/Zb0.060/Sb+0.030i/Sb
????5?4?50.8190/Zb+0.7070i/Zb0.060/Sb+0.020i/Sb
????6?5?60.1872/Zb+0.6188i/Zb0.200/Sb+0.100i/Sb
????7?6?70.7114/Zb+0.2351i/Zb0.200/Sb+0.100i/Sb
????8?7?81.0300/Zb+0.7400i/Zb0.060/Sb+0.020i/Sb
????9?8?91.0440/Zb+0.7400i/Zb0.060/Sb+0.020i/Sb
????10?9100.1966/Zb+0.0650i/Zb0.045/Sb+0.030i/Sb
????1110110.3744/Zb+0.1238i/Zb0.060/Sb+0.035i/Sb
????1211121.4680/Zb+1.1550i/Zb0.060/Sb+0.035i/Sb
????1312130.5416/Zb+0.7129i/Zb0.120/Sb+0.080i/Sb
????1413140.5910/Zb+0.5260i/Zb0.060/Sb+0.010i/Sb
????1514150.7463/Zb+0.5450i/Zb0.060/Sb+0.020i/Sb
????1615161.2890/Zb+1.7210i/Zb0.