Propriedades dos sistemas bidimensionais Lennard-Jones são estudadas com o objetivo de apresentarmos o método de dinâmica molecular. Em particular, damos atenção ao cálculo da função radial de distribuição de pares. A boa concordância de nossos resultados com a teoria e resultados de outra simulação mostra que essé e um método não só de fácil implementação, mas também confiável. Além disso, mostramos que o método de dinâmica molecular pode nos ajudar na interpretação dos resultados e aumentar o nosso entendimento dos mesmos. Palavras-chave: sistema de muitos corpos, Lennard-Jones, distribuição de pares, Argônio, dinâmica molecular.Properties of two-dimensional Lennard-Jones systems are studied with the aim of introducing the molecular dynamics method. In particular, we give attention to the calculation of the radial distribution function. The good agreement of our results with theory and results of another simulation shows that this is a method not only easy to implement, but also reliable. In addition, we have shown that the molecular dynamics method can help us with the interpretation of results and increase our understanding of them.