Wed Oct  8 08:25:36 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     4933.550417        0.000000               0
         1     4933.550417        0.083610     1.69471e-05
         2     4932.876095        0.752486     1.58433e-05
         3     4931.496801        2.103634     1.01383e-05
         4     4929.343102        4.188770    -3.75896e-06
         5     4926.312332        7.102099    -2.75635e-05
         6     4922.330738       10.964650     -5.1693e-05
         7     4917.410877       15.870429    -5.45471e-05
         8     4911.758516       21.799253     1.49008e-06
         9     4905.955477       28.500721     0.000183596
        10     4901.057162       35.349015     0.000578845

    Elapsed wall clock time =  0.8705315589904785  seconds.

md_test():
  Normal end of execution.
Wed Oct  8 08:25:37 2025