We report the melting point of normal alkanes, from molecular dynamics calculations using the NPT ensemble. From the radial distribution function reflecting the arrangement of carbons, the solid–liquid phase transformation was identified in calculating a melting point in good agreement with the measured value. The effectiveness of periodic boundary conditions in the molecular dynamics calculations was also confirmed.