The use of multilinear interpolators to approximate the responses obtained by a rigorous dynamic simulation can be a viable option for the integration of dynamic and steady state models within the framework of an equation-oriented simulator such as EMSO. This work presents a methodology to generate look-up tables that may work as surrogate models of the cellulose enzymatic hydrolysis reactor operating in batch and fed-batch modes. The reaction kinetics was described by Michaelis-Menten model with product inhibition. The input variables of the interpolator for the batch processes were: enzyme concentration, substrate concentration and time. For the fed batch process, an operation time of 120 h was used. From the results, it was observed that for larger batch times the interpolation was more accurate. For the fed batch reactor, it was observed that for low initial concentrations of substrate and enzyme, the results were less accurate. In any case, the proposed methodology was able to approximate the response of dynamic models, allowing the results to be incorporated to steady state simulations of the overall process.