The charge distribution is a basic aspect of electrical transport. In this work we investigate the self-consistent charge response of normal-superconducting heterostructures. Of interest is the variation of the charge density due to voltage changes at contacts and due to changes in the potential. We present response functions in terms of functional derivatives of the scattering matrix. We discuss corrections to the Lindhard function due to the proximity of the superconductor. We use these results to find the dynamic conductance matrix to lowest order in frequency. We illustrate similarities and differences between normal systems and heterostructures for specific examples like a ballistic wire, a resonator, and a quantum point contact.