% Program v Matlabu za izris moči na uporu R2 
P=[]; % prazen array
for R2=0:1:100 % povečujem upornost od 0 do 100 
A=[1,0,0,1,0;-1,1,1,0,0;0,-1,0,0,1;10,0,15j,-5j,0;0,R2,-15j,0,-20j]; % matrika
b=[-j;0;j;100;0]; 
I=inv(A)*b; % resitev tokov, I(2) je tok skozi upor RP
PR2=0.5*I(2).*conj(I(2))*R2 % izracun moci
P=[P PR2] % shranjevanje vrednosti moci v vektor P
end
plot(0:1:100,P) % izris
xlabel('R2 / Ohm')
ylabel('Moc na R2 / W')

%%