ToolPythonIntermediate Solve recurrences with SymPy's rsolve: fixation in the Moran model char = sp.powsimp(sp.expand(eq.replace(phi, sp.Lambda(k, x**k)) / x**(k - 1)))
ToolPythonAdvanced Von Neumann stability analysis with SymPy: the leapfrog Courant limit P = sp.powsimp(sp.expand(mode / (G**n * sp.exp(sp.I * k * j * dx))))expr = sp.powsimp(sp.expand(mode / (G**n * sp.exp(sp.I * k * j * dx))))