Tue May 20 22:02:57 2025

md_test():
  python version: 3.10.12
  numpy version:  1.26.4
  md() is a molecular dynamics program.

md():
  A molecular dynamics program.

  D_NUM, the spatial dimension, is 3
  P_NUM, the number of particles in the simulation is 100.
  STEP_NUM, the number of time steps, is 10.
  DT, the time step size, is 0.1 seconds.

  At each step, we report the potential and kinetic energies.
  The sum of these energies should be a constant.
  As an accuracy check, we also print the relative error
  in the total energy.

      Step      Potential       Kinetic        (P+K-E0)/E0
                Energy P        Energy K       Relative Energy Error

         0     4937.003520        0.000000               0
         1     4937.003520        0.067196     1.36107e-05
         2     4936.462365        0.604766     1.28845e-05
         3     4935.361417        1.688723     9.44289e-06
         4     4933.665169        3.350359     2.43229e-06
         5     4931.342982        5.634867    -5.19981e-06
         6     4928.416578        8.574675    -2.48467e-06
         7     4925.017893       12.136570     3.05738e-05
         8     4921.441948       16.161401     0.000121497
         9     4918.118685       20.324686     0.000291645
        10     4915.453569       24.200834     0.000536942

    Elapsed wall clock time =  0.8224389553070068  seconds.

md_test():
  Normal end of execution.
Tue May 20 22:02:58 2025