The Haber process is the high-pressure synthesis of ammonia (1.5 H₂ + 0.5 N₂ ↔ NH₃). The Peng-Robinson equation of state and the reaction-coordinate method are used to compute the extent of reaction, ξ in moles of NH₃ per unit time, for user-set values of the reactor temperature T (K) and pressure P (bar). This simulation carries out the mass balance around the reactor and splitter in order to compute the amount of ammonia produced when chemical equilibrium is reached in the reactor and in the recycle stream. The splitter separates NH₃  perfectly from the reactants (N₂, H₂), and the fresh feed contains stoichiometric amounts of N₂ and H₂. As expected from Le Chatelier’s principle, the amount of ammonia produced is greater at low temperatures and high pressures.