  We present a guide to compute the phase-boundaries of classical systems using
a dynamic Clausius-Clapeyron integration (dCCI) method within the
LAMMPS(Large-scale Atomic/Molecular Massively Parallel Simulator) code. The
advantage of the dCCI method gives coexistence curves spanning a widerange of
thermodynamic states from relatively short a single non-equilibrium
simulations. We describe the state of art of the non-equilibrium free-energy
methods as allow us compute the Gibss free-energy in wide interval of pressure
and/or temperature and the dCCI method in details, discuss their implementation
in the LAMMPS package and make available source code, scriptsas well as
auxiliary files. As an illustrate example we determine the phase diagram of
silicon in a range of pressures covering from 0 to 15 GPa and temperatures
above the melting points, to obtain the phase boundaries and triple points
between the diamond, liquid and $\beta$-Sn phases are reported.
