Fixation probability by simulation: how often a beneficial mutation survives

Afterwards you can simulate selection and drift with numpy.random, estimate a fixation probability with an error bar, and check it against Kimura's formula.