Save Jacobian While Solving with NDSolve

How does ones save the Jacobian matrix while solving a system of ODEs with NDSolve?

Hi,

You can define a tangent system and solve it with NDSolve, that will give you the Jacobian.
Also, ParametricNDSolve may be of help to you.
I.M.

You can look up NDSolve`EvaluateJacobianWithSolutionData in the tutorial Components and Data Structures.