15-May-2025 13:53:32

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

md():
  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.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    48619.057192    73491.288689    0.000000e+00       22.629454      383.383069        0.000000
        50   123698.269424    24398.987412    2.128150e-01      109.223731      220.902637       99.810845
       100   124573.668578    23486.260980    2.125093e-01      213.899172      216.731451       31.453355
       150   124676.708604    23383.116925    2.125084e-01      320.637833      216.255020       18.329064
       200   124701.951227    23358.451669    2.125132e-01      427.969093      216.140934       14.332262
       250   124713.495149    23347.653628    2.125193e-01      535.552111      216.090970       11.923611
       300   124719.989140    23341.685937    2.125236e-01      643.263646      216.063352       10.715223
       350   124722.310590    23339.843310    2.125275e-01      751.049014      216.054823       10.101923
       400   124722.533035    23340.095191    2.125314e-01      858.881251      216.055989       10.036518
       450   124721.284139    23341.859990    2.125356e-01      966.745287      216.064157       10.537136
       500   124721.705073    23342.072527    2.125408e-01     1074.632195      216.065141       11.677495

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

md():
  Normal end of execution.

md_test():
  Normal end of execution.

15-May-2025 15:40:51