Lennard-Jones potential Matlab code

It is to program a realistic Molecular Dynamics simulation of an Argon gas in 2 spatial dimension. Details: - Program a 2D molecular Dynamics simulation using the Lennard-Jones potential for the atom pairs interaction - Use specific Lennard-Jones parameters for Argon - Implement the simulation by using the Verlet algorithm. There are two versions, the Verlet and the Velocity Verlet. Either one is fine to use, but it may be more useful if you make use of the Velocity Verlet algorithm. - Use a cutoff limit for the interaction distance to control run times - Implement periodic boundary conditions - The initial positions and velocities of the atoms can be random. - your simulation is capable of extracting physical quantities (i.e. pressure, speed distribution, etc.).