%PoissonBisect
x=1; y=10;
for i=1:75
    z=(x+y)/2;
    m = max(poisson3(50,50,z));
    if max(m) <80
        x=z;
    else
        y=z;
    end
end
z
