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