lammps案例：液体平衡态rdf计算模拟练习

分享一个2维液态平衡模拟案例代码。
随机生成400*0.8个原子并进行温度初始化，原子在LJ力场作用下运动到一定位置后达到平衡，并输出rdf数据。

代码已经注释，以供参考。
当原子数量超过400*0.9时，体系接近固态，可自行更改参数并结合rdf图像进行验证。
in文件代码如下：

#设置原子个数
variable afrac  equal 0.80
variable npart  equal ${
   afrac}*400
<