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

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -4.034333
         Iterations: 74
         Function evaluations: 195

Finished calculation

C11 = 86.52147145614002 GPa
C12 = 21.779441553649846 GPa
  B = 43.36011818781323 GPa