from scipy.integrate import solve_ivp
import numpy as np
import matplotlib.pyplot as plt
def rhs(t,S):
# S = [[Ss]0, [Sg]1, [D]2, [Gl]3, [Mal]4, [Mlt]5, [Ea]6, [Ea]7]
# [Ss] = [g/kg]
# [D], [Gl], [Mal], [Mlt] = [g/kg]
# [E] = [U/kg]
# Constants
R = 8.31 # [J/mol...