We propose a new scheme to calculate the elastic constants of solids which is based on linear-response theory. Elastic constants are given directly by a single self-consistent calculation (i.e., they are not obtained by numerical differentiation of total energies or stresses). As an illustration, we apply our procedure to the determination of the equilibrium lattice constant, bulk modulus, and pressure derivative of the bulk modulus of silicon.
