A mathematical model of dehydroisomerization of methylcyclopentane in presence of bifunctional catalysts under industrial operating conditions is presented. The model consists of six differential equations. The rate constants are expressed as polynomial functions of the catalyst blend, and coefficients of these polynomials were determined from numerous measurements, taken along the axis of a tubular reactor. The resulting polynomial coefficients led to reaction rate constants which yielded results close to the actual process.