To overcome these obstacles, we provide a newly developed routine solely relying on opensource software. The code based on Python 3.7, NumPy, Matplotlib, pandas, and SciPy accepts a plain-text file of the chemical reaction set as input and automatically generates the corresponding matrix of coupled differential equations. Its outcome has been validated against the MATLAB implementation of the reaction set of Schneider et al., revealing excellent agreement in both temporal (Figure illustrates that this is achieved by using natural language processing via regular expressions (regex). Next, the reaction set is translated into a matrix to provide suitable input for a numerical ODE-solver (Figure ).