This repository offers basic implementations in Python of Runge-Kutta methods (explicit and implicit) for educational purposes. A damped Newton method is also available with efficient Jacobian reuse. A few test cases are given (heat equation, spring-mass system, Curtiss-Hirschfelder).

Dynamic time step adaptation is implemented for explicit embedded methods.

A script is also provided to study the stability domains, order stars and relative precision contours of any Runge-Kutta scheme. Here is the result for the well-known Radau5 scheme: