ybuswithmutual.m

clc;
clear all;
disp( ' formation of bus admittance matrix ')
nb = input ( 'enter the no of buses including the slack  bus ' );
nl = input ( 'enter the no of lines' );
mb=input ('enter the no of mutal links ');

for j= 1:nl
    sb(j)=input ('enter the starting bus');
    eb(j)=input ('enter the ending bus  ');
    z(j) =input ('enter the input values');
end
y=zeros(nl,nb);
for i=1:nl
    k=sb(i);
    m=eb(i);
    y(i,k)=1;
    y(i,m)=-1;
end
y(:, 1) = [];
u=zeros(nl,nl);
for i=1:nl
    u(i,i)=1/z(i);
end
for j=1:mb
    e=input ('enter the first line');
    f=input ('enter the second line');
    g=input ('enter the value'); 
       u(e,e)= z(f)/(z(e)*z(f)-g^2);
       u(f,f)=z(e)/(z(f)*z(e)-g^2); 
       u(f,e)=-g/(z(f)*z(e)-g^2);
       u(e,f)=u(f,e);
end
t=y';
x=t*u;
x;