close all
clear all
clc

Zimp= [2.5 1.5 0.5 3 2 1 2 5.5];   %Impedancia
Zv= diag(Zimp)      %diag. matica vetových impedancií


C= [0 1 0 0 0 0 0 0;   %rozsírena druhá incidencna mat.
    0 0 0 0 0 0 0 1;
   0 0 0 0 0 0 1 0;
    0 -1 0 -1 0 0 -1 -1;
   -1 0 -1 0 -1 -1 1 1;
   0 0 0 0 1 1 -1 -1;
   -1 0 0 0 0 0 0 0;
   0 0 0 0 -1 0 0 1;
   ]
       
Zs= transpose(C)*Zv*C %matica sluckovych impedancii

Z11= Zs(1:2,1:2)
Z12= Zs(1:2,3:6)
Z13= Zs(1:2,7:8)
Z21= Zs(3:6,1:2)
Z22= Zs(3:6,3:6)
Z23= Zs(3:6,7:8)
Z31= Zs(7:8,1:2)
Z32= Zs(7:8,3:6)
Z33= Zs(7:8,7:8)

Izat=[-5; %matica zatazovych prudov
      -55;
      -40;
      -39;]

Inap= [10;
        5;
] %matica uzlových napajacích napätí

Ins= inv(Z33)*(-Z31*Inap-Z32*Izat) %matica prúdov v lin. nezávislých sluckach

Ics= [Inap; %celkova matica sluckovych prúdov
      Izat;
      Ins;]
  
Iv= C*Ics; %prúdy vo vetvách

deltaU= Zv*Iv;  %úbytky napätí vo vetvách