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