clear all
close all
clc

Z= [2 3 2 1 0.5 5.5 1.5 2.5]; % impedancie vo vetvách vetiev

Y= inv(diag(Z)); % diagonálna matica vetvových admitancií

Kr= [ 1 0 -1 0 0 0;  % redukovaná prvá incidenčná matica
 1 0 0 -1 0 0; 
0 0 -1 0 0 0; 
0 0 -1 0 1 0;
0 0 0 -1 1 0;
0 0 0 0 -1 1;
0 0 0 1 0 -1;
0 1 0 -1 0 0;
;];

Inap=[ 5; 
 55;];

Izat= [ -5; % stĺpcová matica uzlových zaťažových prúdov
 -55;
 -39;
 -40;];

Yadm= transpose(Kr)*Y*Kr % uzlová admitančná matica
Y11= Yadm(1:2,1:2);
Y12= Yadm(1:2,3:6);
Y21= Yadm(3:6,1:2);
Y22= Yadm(3:6,3:6);

Icelkove= [Inap; % matica uzlových prúdov
Izat; ]
U= inv(Yadm)*Icelkove

deltaU= Kr*U % matica úbytkov napätí vo vetvách
I= Y*deltaU % matia prúdov vo vetvách