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

%Zostavenie rovnic MSP v maticovom tvare:
A=[R(1)+R(2),           -R(2),                0;
  -R(2)+R(4),        R(2)+R(3)+R(4),          0;  
     0,              -R(4),               R(5)+R(4)];
 
B=[U(1);        0;               -U(2)];

%Vypocet sluckovych prudov:    A*Is=B
Is=A\B;

%Vypocet vetvovych prudov:
I(1)=Is(1);
I(2)=Is(1)-Is(2);
I(3)=Is(2);
I(4)=Is(2)+Is(3);
I(5)=Is(3);



return

