function?f=myfile（nod8）?
clc
clear?all;
format?long;
Z=[0??1??0.54+j*0.65??0.6+j*0.45
???1??2??0.62+j*0.5??0.4+j*0.3
???2??3??0.6+j*0.35??0.4+j*0.28
???1??4??0.72+j*0.75??0.6+j*0.4
???4??5??1+j*0.55??0.4+j*0.3
???4??6??0.65+j*0.35??0.5+j*0.35
???1??7??0.9+j*0.5????0.5+j*0.4];%33nodes?datas?第一列為首節點編號?第二列為末節點編號?第三列為支路參數?第四列為節點負荷

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%開始計算
V=ones（71）;%賦節點電壓值
for?i=1:7
????I=（Z（:4）.‘）*（diag（1./V））;
????%V1=V;
????I=I‘;
????I1=I;
???for?k=7:-1:1???%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%后推支路電流及功率分布?
??????A=（find（Z（:1）==k））;%判斷矩陣Z第1列是否存在值為k的元素，并把行值賦給矩陣A，亦即找出網絡的末節點
????????if?size（A）==0??%為末節點時
???????????I（k1）;??????
????????else??%不為末節點
???????????I（k1）=sum（I（A1））+I（k1）;?%利用KCL求支路電流
????????end
???end
???V（11）=10.5-Z（13）*I（11）;%%%%%%%%%%%