v2 = [50.5 124.1 438 808 1808 2086 2126 2134.5 2134 2120 2176]/1000;
r2 = [3.2 10.1 20 32.9 50.9 98.3 195.8 293.1 501.8 983 4628 9952];
rs = 62.5;
rp = 62.16;
isc = 500;
vd = 2000;
pv = zeros(1,11);
pi = zeros(1,11);

for j = 1:11
    pv(j) = (isc - v2(j)/(rs+rp))^2*r2(j);
end

for j = 1:11
    pi(j) = (vd - rs * (isc - v2(j)/(rs+rp))) * (isc - v2(j)/(rs+rp));
end

pf = zeros(1,11);

for j = 1:11
    pf(j) = min(pi(j), pv(j));
end
pf