The present work discusses the mathematical modeling and the control design of the steam reforming of natural gas. The developed model comprises a set of differential and algebraic equations, based on energy and mass balances for reactions performed in a fixed catalyst bed reactor, where natural gas and water are transformed mainly into a mixture of hydrogen and carbon oxides. Normally, after removal of hydrogen and purification of the output stream, the residual gas can be also directed to the furnace to provide heat to the reactor. This is a common practice in industrial sites in order to minimize losses. As the global reactions are exothermic, the reactor temperature may reach prohibitive high values, leading to coke formation and catalyst deactivation. For this reason, a control scheme is proposed to account for regulation of the reactor outlet temperature, using residual and fuel gas streams as manipulated variables, allowing the analysis of effect of several process variables in reactor performance. The obtained results indicate that the proposed mathematical model can accurately represent the steam reforming process and that the proposed control scheme can allow for efficient operation of the reactor, even when the residual gas stream is not sufficient to reach the desired operation temperature.