Concentrations on chemical species from mole balances are solved for 1, 2, and 4 species in Python with the Scipy.Integrate package ODEINT.