Calculate the integral (density * dl), where l is the beam length when the beam passes through the volume. The algorithm comes down to which integration method, how you interpolate densities, interpolate between rays, and if so, how and which optimal set of rays to use if you use interpolation. And of course, how do you imagine geometry.