r=zeros(9,7);p=0.2628;
for l=1:9
    x=0;
    for k=1:7
        x=x+1;
        suc=0;
        for i=1:1000
            count=0;
            count1=0;
            for j=1:8
                v=1;
                for h=1:5
                    if v==1
                       a=rand();
                       if a<p^k
                          v=0;
                       end
                    end
                end
                if v==1
                   count=count+1;
                end
            end
            for j=8:8+l-1
                v=1;
                for h=1:5
                    if v==1
                       a=rand();
                       if a<p^k
                          v=0;
                       end
                    end
                end
                if v==1
                   count1=count1+1;
                end
            end
            if count1+count>=8
               suc=suc+8;
            else
                suc=suc+count;
            end
        end
        r(l,x)=suc/8000;
    end
end
i=[0:1:6];
plot(i,r(1,:),'-o',i,r(2,:),'-*',i,r(3,:),'-+',i,r(4,:),'-d',i,r(5,:),'-p',i,r(6,:),'-o',i,r(7,:),'-o',i,r(8,:),'-o',i,r(9,:),'-o')
axis([0,6,0,1.2]);
legend 0 1 2 3 4 5 6 7 8
grid;