Формула Бернулли:
Р=Р(n,0)+P(n,1)
n=1000, p=0.003, q=1-p
Возможно, при вычислении придется воспользоваться приближенными формулами Пуассона или Лапласа.