To implement the algorithm for a real-life reaction optimization problem, an esterification reaction was selected as a model reaction, which is one of the most essential reactions in chemical and pharmaceutical industry. The selected reaction scheme is shown in Figure , where 1,4-benzenedimethanol (1) reacts with acetic anhydride (2) to form the desired 4-hydroxymethylbenzyl acetate (3) and undesired 1,4-bis(acetoxymethyl)benzene (4).