Numerical solution of nonlinear and non-isothermal general rate model of reactive liquid chromatography