Abstract A model of radiation by multi-atoms is discussed. System A, composed of N atoms with two energy levels, is coupled with a single-mode radiation B. A and B are in turn coupled with two thermal reservoirs D and C respectively, and we also assume that all atoms which are located within a wave length of the radiation field have the same phase. Then under the approximation of long wave, rotating wave and mean fields, by introducing Markov approximation after using the method of density matrix, we have derived the of equations motion for systems A and B respectively. These equations are concerning with the following processes: the spontaneous emission of single atom, the absorption and induced emission stimulated by thermal radiation, the induced emission of N atoms stimulated by single mode coherent state radiation as well as the coherent spontaneous emission of this multi-atoms system. The latter two processes are actually the two kinds of superradiance originated from Dicke model and had been investigated previonsly. We have analysed these five processes by numerical calculations. Besides, the influence of these processes upon the evolution of atomic system attaining thermal equilibrium and the increase of entropy is also discussed. It is pointed out that the correlations inside a multi-atoms system would not be destroyed by coherent spontaneous emission, in contrast with the incoherent one which leads to the change of entropy in atomic system. However, only when the entropy of the two reservoirs has been taken into account, will it be possible to demonstrate the mono-tonic increase in total entropy.