Abstract: Zone method has been commonly used in radiative heat transfer calculation in industry, but it has difficulty in calculating multi integration for solution of direct radiative exchange areas, and also in treating with non uniformal gas absorbing and scattering. this paper presents Monte-Carlo method for solution of direct radiative exchange areas, compiling the corresponding program such areas between the zones of enclosed square system filled with uniform absorbing medium are calculated, and comparing the numerical results with analytical ones.