Enter a species  (e.g. Al, Fe): Dy
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): 3.734641352234526e-10, 6.098765495196522e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -11.492281
         Iterations: 70
         Function evaluations: 184

Finished calculation

C11 = 397.4378400266734 GPa
C12 = 106.3673169528909 GPa
  B = 203.39082464415173 GPa