n the Molecular Dynamics program as the below code, replace the velocity Verlet algorithm by the so-called Euler algorithm in which the positions and velocities are updated according to the following rule:

Note that the difference between these algorithms seems quite small. Then, run two simulations, one using the original program and the other using the Euler algorithm program. Compare the dependence of the systemss total energy on time, as calculated by these two programs. Which algorithm would you rather use in your reserach?

For each simulation you should plot the dependence of the kinetic, potential and total energy (all per particle) on time. You should also turn in the relevant portion of the code, i.e. the main loop that updates the positions, velocities and forces.

(6t)2 (6t)2