n=2000
p=0,004
np=8
р(брак=0)=1/1 * exp(-8)
p(брак=1)=8 * exp(-8)
p(брак=2)=32 * exp(-8)
p(брак=3)=256/3 * exp(-8)
р(А)=p(брак<=3)=(1+8+32+256\3)*exp(-8)=0.042
=============================================================
допустим, партию послали в магазин.
покупатель тянет детали.

2000 --->>> 50:
с(2000 50)= 2*10^100
вероятность появления среди них бракованных:
р(БРАК=0)=1/с(2000 50)=5,01*10^-101
p(БРАК=1)=с(2000 1)/с(2000 50)=10^-97
p(БРАК=2)=с(2000 2)/с(2000 50)=10^-94
p(БРАК=3)=с(2000 3)/с(2000 50)=6.67*10^-92

р(В)=р(среди 50 0 или 1 БРАК)=р(БРАК=0)+р(БРАК=1)=10^-97
р(БРАК0,1,2,3)=р(БРАК=0)*р(брак=0)+...+р(БРАК=3)*р(брак=3)=1,96*10^-92

но поскольку больше 3х браков быть не может, пересчитываем через условную вероятность...:

р(В|А)=р(ВА)/р(А) = р(БРАК0,1,2,3)*р(В)/р(A)= р(БРАК0,1,2,3) =

похоже на правду?