Derive a Jacobian with SymPy and pass it to solve_ivp with lambdify

Afterwards you can derive a Jacobian with SymPy, simplify it, turn it into a NumPy function with lambdify, and pass it to solve_ivp as jac.