function sunni(maxp)
f1=0.10;
f2=0.05;
f3=0.01;
for p=1:maxp
    s1(p)=s(f1, p);
    s2(p)=s(f2, p);
    s3(p)=s(f3, p);
end
ip=1:maxp;
plot(ip,s1,ip,s2,ip,s3);
legend('f=10%','f=5%','f=1%');
xlabel('Number of processors used');
ylabel('Speedup');

function res = s(f, p)
res = (f + (1-f) * g(p)) / (f + (1-f) * g(p)/p);

function res = g(p)
res = p ^ (2/3);