A first-principles method is presented to calculate elastic constants up the fourth order of crystals with cubic and hexagonal symmetries. The relies on numerical differentiation second Piola-Kirchhoff stress tensor a density functional theory approach Cauchy for list deformed configurations reference state. number strained required independent second, third, 24 37 symmetries, respectively. Her...