The numerical method for solving the mathematical model of chromatography process coupled with implicit isotherm has been proposed. The exemplary predictions of elution band profiles were performed for competitive adsorption data of 2-phenylethanol and 3-phenylpropanol on ODS-silica with methanol-water as the mobile phase. The simulations of chromatography process with various isotherm models taking into account lateral interactions in adsorbed phase and surface heterogeneity were discussed.