I have a non-linear ordinary differential equation system of dimension 4. I have found one steady state of interest for this system and I would like to test numerically the stability of this steady state. Analytically the Jacobian would be very large, thereby I am looking for a numerical way to study the stability of steady states for such system of ODE. Do you know if there exists such a numerical method ?
Thank you.