[lammps-users] elastic constants?

Hello,
does anyone, by a chance, has a code to calculate elastic constants
and bulk modulus in LAMMPS?
If not, what would be the best way to add it? As a "fix" command?

Marcin