ToolPythonIntermediate Linear ODE systems with SymPy's dsolve: the Bateman function of an oral dose print(sp.limit(C, k_a, k_e))sp.limit(sp.exp(-k_plain * n_plain * tau_plain), n_plain, sp.oo)C_equal = sp.limit(C, k_a, k_e)
ToolPythonIntermediate Solve recurrences with SymPy's rsolve: fixation in the Moran model print(sp.limit(sol, r, 1))sp.limit(phi1, N, sp.oo)print("beneficial, r = 1 + s: ", sp.limit(phi1.subs(r, 1 + s), N, sp.oo))print("harmful, r = 1/(1+s):", sp.limit(phi1.subs(r, 1 / (1 + s)), N, sp.oo))