Но, в серии n испытаний будет такое:
c = 1 - a * n
v = p + ( 1 - c ) * ( 1 - c ) * ( 1 - c ) * k
событие НЕ произойдет, если в rnd(), возвращающем значения от 0 до 1 выкинуть больше v.
Так какая, блин, вероятность события, если проводить серию в 1000 испытаний ?
Что-то от данной формулировки задачи случился subj. Я просто не понимаю, что хотел этим сказать автор.
P.S. если гонять авторскую реализацию, то при p = 0.001, k = 0.05, a = 0.0001 событие стабильно получаем при n > 100.
Upd: цитирую того же Злобина: "Обычно солнце восходит каждое утро. Но стоит мне заявить, что реализация моих планов зависит от того, взойдет ли завтра утром солнце, и я узнаю, что с вероятностью 0.5 солнце завтра НЕ ВЗОЙДЕТ."
Так вот, мне надо показать пользователю, как автор исходного кода представляет себе вероятность восхода солнца завтра.