The linear equations method is proposed to calculate the complete modal content of the partially coherent laser beam using only the intensity information. This method could give not only the incoherent expansion coefficients of the modal decomposition but also the cross-correlation expansion coefficients using the intensity profiles in several planes of finite distance along the propagation direction. A simulation is also presented to verify the validity of this theory. In our algorithm, the minimum and maximum mode orders should be known a priori, so we provide an estimation method for the two parameters.