07-Jan-2022 23:02:34

md_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test md().

MD
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 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    47267.968454    76180.769741    0.000000e+00       22.254108      390.335163        0.000000
        50   123788.201652    25433.726464    2.087765e-01      113.935309      225.538141      103.508241
       100   124565.583041    24626.496436    2.085347e-01      221.577815      221.930153       33.649327
       150   124682.647413    24508.832273    2.085298e-01      331.015538      221.399333       17.198647
       200   124708.402932    24483.472602    2.085330e-01      440.983347      221.284760       12.500238
       250   124719.985887    24472.586823    2.085387e-01      551.179360      221.235561       10.891803
       300   124727.605253    24465.296717    2.085413e-01      661.491816      221.202607        9.944385
       350   124725.960607    24467.358326    2.085447e-01      771.870584      221.211927       10.144949
       400   124724.418788    24469.455013    2.085492e-01      882.293329      221.221405        9.822649
       450   124724.782289    24469.589358    2.085532e-01      992.747505      221.222012        9.181400
       500   124727.387300    24467.451167    2.085570e-01     1103.224497      221.212347        9.335848

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

MD
  Normal end of execution.

md_test():
  Normal end of execution.

07-Jan-2022 23:02:48