ToolPythonBeginner SymPy from the ground up: where the pendulum's 1.74 % comes from theta_of_phi = 2 * sp.asin(k * sp.sin(phi))cos_theta0 = sp.expand_trig(sp.cos(2 * sp.asin(k)))
ToolPythonAdvanced Von Neumann stability analysis with SymPy: the leapfrog Courant limit print(f"C = 1.1 unstable from kΔx = {sp.N(2 * sp.asin(1 / sp.Rational(11, 10)), 4)}")closed = 2 * sp.asin(C * sp.sin(theta / 2)) / (C * theta)