P=[]; % prazen array
for R2=0:1:50 % 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 R2

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
axis([0 100 0 50]); axis off
hold on
arrow3([0 0],[100 0])
arrow3([0 0],[0 50])
xlabel('R2 / Ohm')
ylabel('Moc na R2 / W')