An analytical benchmark and a Mathematica program for MD codes: testing LAMMPS on the 2nd generation Brenner potential