a=5;
b=100;
for i=1:100
    c=(a+b)/2;
    if max(max(poisson6(0,4,0,4,20,20,c)))<80
        a=c;
    else
        b=c;
    end
end
c