This code refers to the work A necessary condition for non oscillatory and positivity preserving time-integration schemes by Thomas Izgin, Philipp Öffner, Davide Torlo available in arXiv:2211.08905.
In this repository, there are the codes used for the work for Modified Patankar schemes on the connection between oscillations studied in Torlo et al. APNUM 182 (2022) and Lyapunov stability function studied in Izgin et al. ESAIM: M2AN 56 (2022).
The numerical study with the related results can be found in the Julia Jupyter Notebook Oscillation Numerical Bound.
The stability functions and their stability regions are in the Maple notebook