In order to perform the large-scale molecular dynamics simulation of the Yukawa system, a mathematical expression for molecular dynamics using the fast multipole method is described. The model simulations are also performed to test
the performance of our implementation of the FMM
