ToolPythonBeginner SymPy from the ground up: where the pendulum's 1.74 % comes from theta_dot_phi = theta_dot_phi.subs(sp.Abs(sp.cos(phi)), sp.cos(phi))
ToolPythonAdvanced Von Neumann stability analysis with SymPy: the leapfrog Courant limit print(sp.Abs(sp.exp(sp.I * k_plain * dx_plain)), " vs ", sp.Abs(sp.exp(sp.I * k * dx)))