...; y =complex(r,-x); epsilon =0.001; r =epsilon+1; epsilon; ybus =zeros(nbus,nbus); for k=1:nbr if n1(k)>0&a...