. In this paper, we develop a quadratic spline collocation method for integrating the nonlinear partial differential equations (PDEs) of a plug flow reactor model. The method is proposed in order to be used for the operation of control design and/or numerical simulations. We first present the Crank-Nicolson method to temporally discretize the state variable. Then, we develop and analyze the proposed spline collocation method for the spatial discretization. The design of the collocation method is interpreted as one order error convergent. This scheme is applied on some test examples, the numerical results illustrate the efficiency of the method and confirm the theoretical behavior of the rates of convergence.