In conclusion, the algorithm by Wall allows for compact and fast MATLAB and Python implementations to study systems involving multiple equilibrating components. The performances are superior with respect to more general solvers readily available to the scientific community, enabling among others the study of complex kinetic processes where the system quickly equilibrates its components as the reaction progresses. We suggest that equpy could be used to handle problems where a slow solver and/or a software that relies on poorly formatted input would be a bottleneck.