function amdahl(maxp)
f1=0.10;
f2=0.05;
f3=0.01;
for p=1:maxp
    s1(p)=1/(f1+(1-f1)/p);
    s2(p)=1/(f2+(1-f2)/p);
    s3(p)=1/(f3+(1-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');