function?[?Y?]?=?lfybus（?zdata?）
%UNtitleD2?此處顯示有關此函數的摘要
%???此處顯示詳細說明
nl=zdata（:1）;nr=zdata（:2）;R=zdata（:3）;X=zdata（:4）;B=zdata（:5）;A=zdata（:6）;
nbr=length（zdata（:1））;nbus=max（max（nl）max（nr））;
Z=R+j*X;
y=ones（nbr1）./Z;
Y=zeros（nbusnbus）;
for?k=1:nbr;