% Parameterek [V, KOhm, mA]
R = 2.73;
n = 1.89;
I0 = 6;
Rt = 2*R;
% Egyutthato matrix
A = [3*R 2*R 0 0 0 1 0 0;2*R 3*R+Rt -R 0 0 0 0 0;...
  0 -R 3*R 0 1 0 0 0;0 0 0 R -1 1 0 0;0 0 0 0 1 -n 0 0;...
  0 0 0 0 0 0 n 1;0 0 -1 1 0 0 1 0;1 0 0 1 0 0 0 -1];
% Konstansok vektora (gerjesztes)
B = [I0*R;0;0;0;0;0;0;0];
% Megoldas 
x = A \ B;
% Eredmeny extrahalas
j2 = x(2); 
ut = -j2*Rt
Pt = j2^2*Rt

%% Rovidzaras
% Egyutthato matrix
A = [3*R 2*R 0 0 0 1 0 0;2*R 3*R -R 0 0 0 0 0;...
  0 -R 3*R 0 1 0 0 0;0 0 0 R -1 1 0 0;0 0 0 0 1 -n 0 0;...
  0 0 0 0 0 0 n 1;0 0 -1 1 0 0 1 0;1 0 0 1 0 0 0 -1];
% Konstansok vektora (gerjesztes)
B = [I0*R;0;0;0;0;0;0;0];
% Megoldas 
xrz = A \ B;
IN = -xrz(2)

%% Szakadas
% Egyutthato matrix
A = [3*R 2*R 0 0 0 1 0 0;0 1 0 0 0 0 0 0;...
  0 -R 3*R 0 1 0 0 0;0 0 0 R -1 1 0 0;0 0 0 0 1 -n 0 0;...
  0 0 0 0 0 0 n 1;0 0 -1 1 0 0 1 0;1 0 0 1 0 0 0 -1];
% Konstansok vektora (gerjesztes)
B = [I0*R;0;0;0;0;0;0;0];
% Megoldas 
xuj = A \ B;
UT = (-xuj(3))*R+2*R*xuj(1)

%% 
RB = UT / IN
Pmax = IN^2*RB/4
UT^2/(4*RB)