%initialing?the?parameters
clear;
%input（‘Please?input?the?parameters?of?the?Line?and?the?Transformers‘）;
Nline=input（‘Please?input?the?number?of?the?bus\nNline=‘）;
Ntrans=input（‘Please?input?the?number?of?the?transformer\nNtrans=‘）;
Line=zeros（Nline5）;
Transformer=zeros（Ntrans6）;
for?m=1:Nline
????Line（m:）=input（‘Please?input?the?parameters?of?the?bus\nLine=?‘）;
end
for?m=1:Ntrans
????Transformer（m:）=input（‘Please?input?the?parameters?of?the?transformer\nTransformer=?‘）;
end
clc;
PQ0=input（‘Please?input?the?value?of?the?PQ0\nPQ0=‘）;
PQ0=PQ0‘;
EF=input（‘Please?input?the?value?of?the?EF\nEF=‘）;
EF=EF‘;
%initial?some?variables?that?will?be?used?later
Nnode=Nline+Ntrans;
k=0;
error=1;
Sslack=0;
%form?the?Matrix?of?Y
Y=zeros（Nnode）;
if?Nline>=1???????????????%compute?by?adding?the?nodes?of?the?lines
????for?k=1:Nline
????????t1=Line（k2）;
????????t2=Line（k3）;
????????Y1=1/（Line（k4）+j*Line（k5））;
????????Y（t1t1）=Y（t1t1）+Y1;
????????Y（t1t2）=Y（t1t2）-Y1