We present a density functional scheme for calculating the frequency‐dependent linear response of superconductors. The central result is a set of integral equations determining the linear response of the normal and anomalous densities to external perturbations. Analytic solutions of these integral equations are obtained for homogeneous systems with separable effective interactions. For inhomogeneous superconductors, the formalism leads to a scheme for calculating the critical temperature without explicitly solving the gap equation. © 1995 John Wiley & Sons, Inc.