Lateral adsorbate-adsorbate interactions result in variation of the desorption rate constants with coverage. This effect can be studied in great detail from the shape of a multi-isotherm. To produce the multi-isotherm, the temperature is increased in a (semi)stepwise fashion to some temperature, followed by maintaining this temperature for a prolonged time. Then, the temperature is stepped to a higher value and held constant at this new temperature. This cycle is continued until all of the adsorbates have desorbed. Using a detailed kinetic Monte Carlo model and an optimization algorithm based on Evolutionary Strategy, we are able to reproduce the shape of the experimentally measured multi-isotherm of nitrogen on Rh(111) and obtain the lateral interactions between the nitrogen atoms.