配电网潮流计算程序matlab

clear
n=40;%n=input（‘请输入节点数‘）;
b=39;%b=input（‘请输入支路数‘）;
%disp（‘请输入支路阻抗和节点功率矩阵，第一列存支路号，第二列存首节点号，第三列存尾节点号，第四列存支路自阻抗第五列存尾节点给定功率‘）;
Sb=15;
Ub=23;
Zb=Ub^2/Sb;
Z=[1010.169/Zb+i*0.655/Zb0
2120.864/Zb+i*0.751/Zb0
3130.196/Zb+i*0.655/Zb0
4241.374/Zb+i*0.774/Zb0.522/Sb+i*0.174/Sb
5250.864/Zb+i*0.751/Zb0
6360.444/Zb+i*0.439/Zb0
7370.196/Zb+i*0.655/Zb0
8580.864/Zb+i*0.751/Zb0
9590.864/Zb+i*0.751/Zb0
106101.374/Zb+i*0.774/Zb1.08/Sb+i*0.36/Sb
116110.864/Zb+i*0.751/Zb0
127120.279/Zb+i*0.051/Zb0
137130.279/Zb+i*0.051/Zb0
148141.374/Zb+i*0.774/Zb0.36/Sb+i*0.12/Sb
158151.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
169161.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
179170.864/Zb+i*0.751/Zb0
1811181.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
1911191.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2012201.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
2112210.444/Zb+i*0.439/Zb0
2213220.444/Zb+i*0.439/Zb1.8/Sb+i*0.6/Sb
2313230.444/Zb+i*0.439/Zb0
2417241.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2517251.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
2621261.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
2721270.444/Zb+i*0.439/Zb0
2823281.374/Zb+i*0.774/Zb0.855/Sb+i*0.285/Sb
2923290.864/Zb+i*0.751/Zb0
3027301.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
3127310.864/Zb+i*0.751/Zb0
3229321.374/Zb+i*