Hybrid Monte-Carlo simulation algorithm of fluorescence anisotropy decay in the presence of forward and back energy transfer is presented for rigid disordered systems. The results of hybrid Monte-Carlo simulations are compared with the results of SCDM theoretical model and with Monte-Carlo step by step simulations. Very good agreement between theoretical and simulation data as well as Monte-Carlo simulations has been obtained. The hybrid Monte-Carlo simulation of emission anisotropy is characterized by much better accuracy and little time consumption, especially for long times and also for low donor and at the same time high acceptor concentrations. (C) 2008 Elsevier B.V. All rights reserved.