A two-component model of reactive liquid chromatography is presented considering M!N type reaction. The model incorporates surface and pore diffusions in the adsorbates, axial dispersion, interfacial mass transfer, first order chemical reactions in the liquid and particle phases, and two sets of boundary conditions. The model contains a system of four coupled partial differential equations(PDEs) describing the dynamics of reactants and products in both phases. The Laplace transformation and eigen-decomposition technique are jointly applied to solve the model equations analytically. An efficient and accurate numerical Laplace inversion technique is utilized to retrieve back solutions in the original time domain. The developed semi-analytical results are verified against the numerical results of a high resolution finite volume scheme (HR-FVS). A good agreement between the solutions not only confirms the accuracy of semi-analytical results but also validates the accuracy of proposed numerical scheme. This study extends and generalizes our previous analysis on heterogeneous reactions in the liquid chromatography. In order to analyze the behavior of a chromatographic reactor, different case studies are presented showing the effects of model parameters on the process.