Abstract
By using the molecular dynamics (MD) simulation, we study the second order phase transition of a model of crystalline membranes. We employ three kinds of techniques to solve the equation of motion : first two are based on the global coordinate system and the third one is based on local coordinate systems. The results are compatible with those obtained by Monte Carlo (MC) simulations.