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

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -3.835971
         Iterations: 72
         Function evaluations: 186

Finished calculation

C11 = 84.74044287034224 GPa
C12 = 21.50763328116767 GPa
  B = 42.58523647755919 GPa