The moment equations of the concentration distribution for the multi-region model are derived using the method of moment. The method originally devised by Aris is to obtain the concentration moments satisfying a given PDE (Partial Differential Equation. The method of moment is used to obtain the first five moments (0th to 4to) that satisfy the model PDE. Each moment of the concentration distribution for the model equation is plotted for the dimensionless time and gave similar results except the skewness and the kurtosis. The results of the analysis show the physical meaning of each moment. The comparisons with the number of regions or the global interaction coefficient give a possibility to determine the parameters of the multi-region model with the analytical concepts.