Само применение приближённой интегральной формулы Муавра-Лапласа здесь некорректно, как бы точно потом не вычислялось всё остальное. Для малых p точнее будет работать приближённое распределение Пуассона.