08-Oct-2025 10:06:26

md_fast_test():
  MATLAB/Octave version 6.4.0
  Test md_fast().
08-Oct-2025 10:06:26

MD_FAST
  MATLAB/Octave version 6.4.0
  A molecular dynamics program.

  ND, the spatial dimension, is 3
  NP, the number of particles in the simulation is 500.
  STEP_NUM, the number of time steps, is 500.
  DT, the time step size, is 0.100000 seconds.

  Initializing positions, velocities, and accelerations.

  Computing initial forces and energies.

  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   124374.076322      771.804782    0.000000e+00
        50   124675.518655     1161.922266    5.526029e-03
       100   124739.201442     1159.465692    6.015268e-03
       150   124749.536034     1158.701033    6.091738e-03
       200   124749.104000     1160.489450    6.102577e-03
       250   124750.000000     1161.041187    6.114145e-03
       300   124749.851534     1161.422578    6.116006e-03
       350   124749.938333     1161.592940    6.118061e-03
       400   124750.000000     1161.781644    6.120062e-03
       450   124750.000000     1161.781644    6.120062e-03
       500   124750.000000     1161.781644    6.120062e-03

  Main computation:
    Wall clock time = 36.898195 seconds.

MD_FAST
  Normal end of execution.

08-Oct-2025 10:07:03

md_fast_test():
  Normal end of execution.

08-Oct-2025 10:07:03