In this contribution, a Monte-Carlo computer simulation of phonon-assisted relaxation processes of excitons in CdSe/ZnSe quantum islands is presented. With the same set of parameters, it was possible to reproduce both the temperature dependence of the energy position and of the full width at half maximum (FWHM) of the luminescence emitted from a ZnCdSe quantum film containing 5 to 10 nm wide Cd-rich islands. With help of the model, the energy dependence of the distribution of the localized states is shown to be weaker than Gaussian. Furthermore, the obtained fitting parameters give evidence for an intra-island excitonic relaxation. Finally, our results, based on a realistic description of the exciton kinetics, are compared with those of a rate equation model. 1 Introduction In localized systems, and among them in quantum-dot systems, thermally activated excitons can transfer to adjacent states with the help of phonon-assisted processes. As a result, the energy position of the photoluminescence (PL) maximum exhibits a "red-blue-red"-shift with increasing temperature [1]. This anormalous shift is well known for various quantum-well systems, where it was observed experimentally [2-4] and described theoretically [5]. For quantum dot systems a rate equation approach has been used [6]. In these systems, the shift has been interpreted in terms of an inter-island relaxation [7]. Interpretation of optical spectra has, however, indicated the fractal nature of the island shape [8,9]. In the resulting potential, intra-island hopping relaxation of excitons would be possible.In this paper, the scenario of an intra-island spatial relaxation is tested. Experimental results for the anomalous PL shift are compared with a Monte-Carlo simulation of phonon-assisted hopping [10]. After a description of the sample and of the experimental data, the operating principle of the simulation and its application to our experimental data are presented. Then the obtained parameter sets and their meaning for the distribution of localization centers for hopping processes will be discussed. Finally, the model is compared with the results of a calculation based on rate equations.