Teoría de la dinámica molecular (1) -Mecánica newtoniana, dinámica hamiltoniana y espacio de fase

Resumen de la teoría de la dinámica molecular

A continuación, dibujamos la trayectoria del resorte armónico en el espacio de fase: El
código MALAB es el siguiente:

clc;clear;
%--------------------------------%
% initial setting 
step = 0.1;
X_initial  = -4;
V_initial = 0;
lambda    = 1;
mass      = 1;
X_current = X_initial;
V_current = V_initial;
t         = 0;run = 0;
%--------------------------------%
% Run
while (t<=100)
  t   = t+step;
  half_t    = 0.5*step;
  a = -X_current;
  V_current = V_current+a*step;
  X_current = X_current+V_current*step;
  momentum  = V_current;
 %---------------------------------------%
 % plot 
  hold on;  
  s = scatter(X_current,momentum);
  s.LineWidth = 0.6;
  s.MarkerEdgeColor = 'b';
  s.MarkerFaceColor = [0 0.5 0.5];
  run = run +1;  
 %---------------------------------------%
 % output 
  if(mod(run,100)==0)
      fprintf("The step %d is finished!\n",run);      
  end 
end
disp("--------------------------");
disp("--------ALL DONE!!!-------");
disp("--------------------------");