% STAVKI
% THEVENIN
Rb=0:0.1:50 	% tvorimo niz vrednosti Rb od 0 do 50 s korakom 0,1
Uth=30.91
Rth=14.32
P=Rb*Uth^2./(Rth+Rb).^2    % Izracun moci 
plot(Rb,P)               % izris
xlabel('Rb     (Ohm)')
ylabel('Pb     (W)')


maxP=max(P)
i=find(P==max(P))
Rb(i)
% break



b=[-2;0;2;10;0] % vektor znanih vrednosti / desna stran enacbe
II=[]  % prazen vektor, potreben za shranjevanje izračunanih vrednosti moči
for R3=0:0.1:50   % zanka povecuje upornosti
    A=[1,0,0,1,0;-1,1,1,0,0;0,-1,0,0,1;20,0,R3,-1,0;0,5,-R3,0,40]
    
    I=A\b     % izracun tokov za dolocen R3
    
    II=[II,I(3)]  % shranjevanje vrednosti toka I3 v vektor, ki se zaporedno polni
end

R3=0:0.1:50   % vektor upornosti
P=II.^2.*R3     % izracun moci
plot(R3,P,'+g')   % izris