n=5;%節點數為5
L=5;%支路數為5
B1=[125.931.50;1531101.05;230.8230.954545;34651000;45651000];
%每行是一條支路
%第一列是支路的一個端點I
%第二列是支路的另一個端點J
%第三列是支路的電阻R
%第四列是支路的電抗X
%第五列是支路的對地電納或變壓器變比k
%規定k*始終在I邊Y=zeros（5）;
G=zeros（5）;
B=zeros（5）;
for?m1=1:L
????I=B1（m11）;
????J=B1（m12）;
????R=B1（m13）;
????X=B1（m14）;
????k=B1（m15）;
????if?J==0????????????%變壓器勵磁支路
????????????G（II）=?G（II）+R/（R^2+X^2）;
????????????B（II）=??B（II）-X/（R^2+X^2）;
????????end
????if?I>0????????????%線路支路
????????????G（IJ）=-R/（R^2+X^2）;
????????????B（IJ）=X/（R^2+X^2）;
????????????G（JI）=G（IJ）;
????????????B（