In the paper, an algorithm for reconstruction of the intensity distribution over the laser beam cross section from the temperature field of a thin target is described. The reasons for the solution instability and methods for their elimination are analyzed. The results of numerical reconstruction experiment are presented for various boundary conditions.