In highly scattering media, light energy fluence rate distributions can be described by diffusion theory. Boundary conditions, appropriate to the diffusion approximation, are derived for surfaces where reflection of diffuse light occurs. Both outer surfaces and interfaces separating media with different indices of refraction can be treated. The diffusion equation together with its boundary conditions is solved using the finite element method. This numerical method allows much freedom of geometry.