При вычислении стоит использовать асимптотическую формулу Пуассона, (т.к. npq<9) - даёт хорошую точность в 0.001, Муавра-Лапласом считать не вздумайте! Погрешность большая, 5 точно не поставят, сразу спалят, что теорию не доучили...

...Можно и Бернулли, но нужен Mathcad, а лучше maple там в 999 и 1000 степень надо будет возводить, благо не будет факториалов, которые замечательно сократятся...