Y BUS MATRIX

 n=input('enter number of buses=');

Zbus=[0 0.15 + 0.6i 0.1 + 0.4i 0.15 + 0.6i; 0.15 + 0.6i 0 0.05 + 0.2i 0;

0.1 + 0.4i 0.05 + 0.2i 0 0.05 + 0.2i; 0.15 + 0.6i 0 0.05 + 0.2i 0]

for m=1:n

 for k=1:n

 if(m==k)

 Ybus(m,k)=0;

 elseif(Zbus(m,k)==0)

 Ybus(m,k)=0;

 else

 Ybus(m,k)=-1/Zbus(m,k)

 end

 end

end

Y_sum=sum(Ybus,2)

for m=1:n

 for k=1:n

 if(m==k)

 Ybus(m,k)= -Y_sum(m)

 end

 end

end