A Monte Carlo code for calculating lidar returns from clouds in regime of multiple scattering is tested, by comparing its results pertaining to second order of
scattering with those obtained by using an analytic formula developed in a completely different way. For obtaining second order of scattering all the essential parts of the
Monte Carlo code are employed. Thus the very good agreement between the results of the two procedures, which have been found in a series of different situations of scattering
media, has to be considered as a positive test of the Monte Carlo code reliability