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