(Оффтоп)
А не обобщить ли проблему на кидающую монетку тётю? А ещё есть грантопригодные темы с изучением монетометания тётями, которые на деле дяди, и дядями, которые тёти...
Вообще же мне видится источник ошибки ТС в том, что, рассчитывая вероятность получить в (n+1)-м броске "орла" при условии , что ранее их выпало k, и совершенно верно приравнивая её к вероятности получить в (n+1) бросках (k+1) "орлов" при условии, что в первых n бросках их выпало k, вычисляет её по вероятности для другой задачи - получить в (n+1) бросках (k+1) "орла", не оговаривая условие, что в первых n бросках их выпало k. Эта величина больше потому, что в неё входит и возможность получить (k+1) "орла" уже в первых n бросках. Что и приводит к завышению вычисленной вероятности по сравнению с правильной p.