Optical properties of biological tissue can change as a result of thermal denaturation due to temperature rise; a familiar example is whitening observed in cooking egg-white. Changes in optical properties with temperature have been reported in the literature. Temperature rise due to laser irradiation is a function of the optical properties of tissue which themselves are a function of temperature of the tissue. This creates a coupling between light and temperature fields for biological tissue under laser irradiation. The effects of this coupling on the temperature response and light distribution may play an important role in dosimetry consideration for therapeutic as well as diagnostic application of lasers in medicine. In a previous study this problem was addressed in one dimension, for short irradiation exposures, using certain simplifying assumptions. The purpose of this research was to develop a mathematical model for dynamic optical changes with thermal denaturation and a computer program for simulation of these effects for a multi-dimensional geometry.