四维理想气体速度分布定律
1，将分子考虑成质量 m = 1，半径 r = 1 的小球，初始时粒子具有某
种速度分布 (均匀分布)
2，蒙特卡洛方法并不需要模拟真实粒子轨迹，假设任意两个粒子
在某处，以某种方式发生了碰撞即可，碰撞足够多的次数后，分子
的速率分布会达到稳定状态
分子数不小于 10的4次方
完全弹性碰撞，系统能量动量要守恒，温度可以由粒子平均平动动
能得出
给出不同温度的粒子速度分布 (作图表示)，并对比理论值
p(v) ∝ v3
exp(−
mv2／
2kT）
尝试并行化程序加快计算速度，如何更好的表示分布函数