S. cerevisiae glucose fermentation with Monod kinetics

S. cerevisiae glucose fermentation with Monod kinetics. Includes product inhibition. Primary and secondary fermentations.

This code allows for a very basic simulation of the fermentation of glucose in batch mode using parameters found in the literature.
The simulation assumes a linear product inhibition (1-P/Pmax).
It also assumes an extent of fermentation of 81%. That means when 19% of the initial substrate is reached, the fermentation and growth of biomass will stop. This can be easily undone by removing an if condition.
The code generates 4 plots:
1. Primary fermentation.
2. Secondary fermentation (lower umax).
3. Both fermentations to scale.
4. Both fermentations to scale, superimposed.

I am a Biotech graduate new to coding, so please leave a comment if you see any mistakes / possible improvements!

References:
Shopska, V., Denkova, R., Lyubenova, V., & Kostov, G. (2019). Kinetic Characteristics of Alcohol Fermentation in Brewing: State of Art and Control of the Fermentation Process. Fermented Beverages, 529–575. https://doi.org/10.1016/b978-0-12-815271-3.00013-0