The subject codes use minimization of Gibbs energy to determine equilibrium composition given reactants and pressure. The literature (decades old NASA docs: SP-273, RP-1311, etc.) describes the process for setting up the matrices to solve for the correction factors needed for each iteration...