vs = 5;
r1 = 1000;
r2 = 0:1:10000;
is = vs./(r1+r2);
v2 = is.*r2;
p2 = is.*v2;
plot (r2, p2)
title ('power dissipated by r2')
xlabel ('r2 in ohm')
ylabel ('power in watt')
grid on
mp = max (p2);
mp_index = find (p2 == mp);
mp_r2 = r2 (mp_index);
fprintf ('Maximum power = %g, occurs at r2 = %g \n', mp, mp_r2);