Enter a species  (e.g. Al, Fe): Pb
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: LJ_ElliottAkerson_2015_Universal__MO_959249795837_003
Enter hexagonal lattice constants 'a' and 'c' (m): 2.839879784973993e-10, 4.637596747018009e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -14.305720
         Iterations: 72
         Function evaluations: 190

Finished calculation

C11 = 1125.1734253579205 GPa
C12 = 301.1330963752054 GPa
  B = 575.8132060361104 GPa