The purpose of this project was to compare the different physical models behind four algorithms in computational chemistry: Molecular dynamics with a thermostat (specifically simple velocity rescaling, Berendsen, and Nosé-Hoover), Langevin dynamics, Brownian dynamics, and Monte Carlo. These algorithms were programmed in C and the impact of specific parameters, such as the coupling parameter and time step, were studied. Their results were compared based on their radial distribution functions and, when the thermostats were in use, fluctuations in temperature.