Я запустил Монте Карло еще пару раз. Среднее значение вероятности получалось проводя 10 экспериментов, в которых создавалось N последовательностей ДНК, длиной в 100 нуклеотидов:
| N | Время работы | Вероятность | Истинное значение (теор.) | Отклонение |
|---|---|---|---|---|
| 100 | 0.074 с | 0.991 | 0.9950793852734019 | 0.004079385273401903 |
| 1 000 | 0.633 с | 0.9949 | 0.9950793852734019 | 0.0001793852734018886 |
| 10 000 | 5.924 с | 0.99506 | 0.9950793852734019 | 1.938527340195062e-05 |
| 100 000 | 59.021 с | 0.9950130000000001 | 0.9950793852734019 | 6.638527340174782e-05 |
| 1 000 000 | 5849.2 с | 0.9950823200000001 | 0.9950793852734019 | -2.9347265981805037e-06 |
В общем, я думаю мы нашли правильное решение))
P.S. 5849.2 секунды это 97.5 минут