In this work, the boundary control of a distributed parameter system modelled by linear parabolic partial differential equations (PDEs) with spatially varying coefficients is studied. An infinite-dimensional state space setting is considered and an exact transformation of the boundary actuation is realised to obtain an evolutionary model. The evolutionary model which incorporates the spatially varying coefficients of the underlying set of the PDEs is used for subsequent linear quadratic regulator synthesis. The formulated linear quadratic-state feedback controller is applied to a nonlinear model of the reactor and its performance is studied.