A real-time algorithm for rendering of an inhomogeneous scattering medium such as fog with a surface object immersed therein is described. An input media animation is represented as a sequence of density fields. The algorithm computes surface reflectance of the surface object in the inhomogeneous scattering medium. The algorithm may also compute airlight of the inhomogeneous scattering medium. Several approximations are taken which lead to analytical solutions of quantities such as optical depth integrations and single scattering integrations, and a reduced number of integrations that need to be calculated. The resultant algorithm is able to render inhomogeneous media including their shadowing and scattering effects in the real time. The algorithm may be adopted for a variety of light sources including point lights and environmental lights.