The increasing demand for optimization in process control of reactive systems necessitates the development of fast and reliable software for the numerical computation of optimal controls taking into account the special structures of the chemical systems under investigation. A recently developed algorithm based on partially reduced SQP methods is used for the computation of optimal controls following the direct approach. In this work, this algorithm is coupled with a package for the simulation of homogeneous reactive systems, which provides the specific routines for the description of the problem at hand. The combined code takes advantage from the fact that the number of influence variables to be determined is rather low compared with the overall number of variables. For example, the computation tool developed is applied to study the homogenous oxidative coupling of methane.