function I=MUN(U,R)
%Vypocet prudov vo vetvach obvodu podla prilozenej schemy na zaklade
%metody uzlovych napati.
%MUN => Metoda uzlovych napeti
%I=MUN(U,R)
%I...vektor vetvovych prudov
%U...vektor napati zdrojov
%R...vektor odporov

%Zostavenie rovnic MUN v maticovom tvare:
A=[1/R(1)+1/R(2)+1/R(3),                      -1/R(3);
   
   -1/R(3),                          1/R(3)+1/R(4)+1/R(5)];    

B=[U(1)/R(1);                      U(2)/R(5)];

%Vypocet uzlovych napati:    A*Uzn=B
Uu=A\B;

%Vypocet vetvovych prudov:
I(1)=(U(1)-Uu(1))/R(1);
I(2)=(Uu(1))/R(2);
I(3)=(Uu(1)-Uu(2))/R(3);
I(4)=(-Uu(2))/R(4);
I(5)=(Uu(2)-U(2))/R(5);

return