This paper describes an algorithm for the simulation of diffuse first reflections using a hybrid of empirical and physical models. The system uses modified components of the Phong shading algorithm developed for three-dimensional graphics models in conjunction with physical models of components of the Schroeder diffusers. The development and analysis of the algorithm are discussed.