RecipePythonIntermediate Derive a Jacobian with SymPy and pass it to solve_ivp with lambdify y = A, B, C = sp.symbols("A B C")k1, k2, k3 = sp.symbols("k1 k2 k3", positive=True)
ToolPythonBeginner SymPy from the ground up: where the pendulum's 1.74 % comes from theta, theta0, theta_dot, L, g = sp.symbols("theta theta0 theta_dot L g", positive=True)t = sp.symbols("t")k, phi = sp.symbols("k phi", positive=True)x = sp.symbols("x")