y=zeros(1,401);
for k=0:400
    a=0;
    b=10;
    for i=1:100
        c=(a+b)/2;
        if max(max(poisson5(0,4,0,4,15,15,c,1+(k/100))))<80
            a=c;
        else
            b=c;
        end
    end
    y(1,k+1)=c;
end
x=1:.01:5;
plot(x,y)