A model for the calculation of the scalar, position-dependent static dielectric function of an ensemble consisting of polarizable sites and of free relaxing point dipoles at fixed positions is ,proposed. It is based on classical electrostatics and leads to an iteration equation or nonlinear partial differential equation for the local dielectric constant. The expressions contain the equations of Debye, Onsager, and Neumann as special cases and thus might be considered as an extension to inhomogeneous matter. The model may have applications in the case of biological macromolecules in particular proteins where the polar side chains can be identified with the model's dipoles. The algorithm leads to a position dependent dielectric constant in th...