24-Dec-2022 11:45:42

md_test():
  MATLAB/Octave version 4.2.2
  Test md().

md():
  MATLAB/Octave version 4.2.2
  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.01 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    49592.454842    70683.995895    0.000000e+00       22.676983      375.989351        0.000000
        50   123854.387615    23743.398091    2.271545e-01      108.976338      217.914653       89.219546
       100   124603.618385    22960.895980    2.268779e-01      212.718509      214.293705       26.707299
       150   124684.281399    22879.582778    2.268725e-01      318.371357      213.913921       17.590421
       200   124707.645569    22856.938057    2.268784e-01      424.581257      213.808036       13.515011
       250   124714.882319    22850.429389    2.268845e-01      531.025420      213.777592       12.146371
       300   124718.107992    22847.746243    2.268890e-01      637.591105      213.765040       10.480620
       350   124721.283798    22845.145666    2.268938e-01      744.229592      213.752874        9.451243
       400   124726.671560    22840.164903    2.268972e-01      850.914203      213.729572        8.851539
       450   124730.023446    22837.117002    2.268997e-01      957.628147      213.715311        8.916556
       500   124728.085071    22839.383037    2.269024e-01     1064.361794      213.725913        8.700796

  Wall clock time = 8929.058893 seconds.
  Graphics saved as "md.png"

md():
  Normal end of execution.

md_test():
  Normal end of execution.

24-Dec-2022 14:14:32