Solve recurrences with SymPy's rsolve: fixation in the Moran model

Afterwards you can solve a linear recurrence with boundary values using SymPy's rsolve, handle its degenerate case with limit, and check it by simulation.