毕业需求 跪求大师 给指点下程序，要不然不毕业我就自杀

楚香饭

楼主的代码有何问题呢？

在我这里可以运行，输入 TIMESTEP 为 300，得到一组变化的图形（我不知道是不是你希望的晶体）

我截了两个时间的图形如下：



另外得到 LUKE.dat 文件，部分如下：

           1   1.000000   
           2   1.732051   
           3   2.449490   
           4   3.464102   
           5   4.472136   
           6   5.099020   
           7   5.744563   
           8   6.480741   
           9   6.855655   
          10   7.681146   
          11   8.366600   
          12   9.273619   
          13   10.53565   
          14   11.57584   
          15   12.40967   
          16   13.30413   
          17   13.92839   
          18   14.93318   
          19   15.81139   
          20   16.91154   
          21   18.05547   
          22   19.13113   
          23   20.04994   
          24   21.11871   
          25   22.00000   
          26   23.04344   
          27   23.93742   
          28   24.95997   
          29   26.13427   
          30   27.11088   
          31   28.12472   

如果楼主对计算结果不满意，可以进行 Debug 调试，查看是哪里的问题。

如果楼主是希望在这里找到对应专业的人来修改计算结果不正确，那难度就很高了。真的。