ConceptPythonBeginner Computer algebra: formulas as trees, exact until you ask for digits M = sp.Matrix(n, n, lambda i, j: sp.Symbol(f"a{i}{j}"))H = sp.Matrix(10, 10, lambda i, j: sp.Rational(1, i + j + 1))
RecipePythonIntermediate Derive a Jacobian with SymPy and pass it to solve_ivp with lambdify rhs = sp.Matrix([-k1 * A + k3 * B * C,