We present an integral density method for calculating the multifractal dimension spectrum for nucleon distribution in atomic nuclei. This method is then applied to analyze the non-uniformity of density distribution in several typical types of nuclear matter distributions, including the Woods–Saxon distribution, halo structure, and tetrahedral α clustering. The subsequent discussion provides a comprehensive and detailed exploration of the results obtained. The multifractal dimension spectrum shows a remarkable sensitivity to the density distribution, establishing it as a simple and novel tool for studying the distribution of nucleons in nuclear multibody systems.