ToolPythonIntermediate Solve recurrences with SymPy's rsolve: fixation in the Moran model print("roots at r = 1:", sp.roots(char.subs(r, 1), x))
ToolPythonAdvanced Von Neumann stability analysis with SymPy: the leapfrog Courant limit print(sp.solve(P_edge, G), sp.roots(P_edge, G))