这组MATLAB例程模拟了具有周期性边界条件的三维硬球堆积。它从一个均匀晶格开始,然后利用Metropolis算法对粒子位置进行多次退火迭代。随后计算成对径向分布函数$g_2(r)$,在集合中多个中心上取平均值。对于较高的初始粒子密度,得到的$g_2$与Percus-Yevick近似解析结果非常匹配,能够以任意精度描述足够高的粒子密度。