Traditional methods for evaluating the low-albedo volume rendering integral do not include bounds on the magnitude of approximation error. In this paper, we examine three techniques for solving this integral with error bounds: trapezoid rule, Simpson's rule, and a power series method. In each case, the expression for the error bound provides a mechanism for computing the integral to any specified precision. The formulations presented are appropriate for polynomial reconstruction from point samples; however, the approach is considerably more general. The three techniques we present differ in relative efficiency for computing results to a given precision. The trapezoid rule and Simpson's rule are the most efficient for low- to medium-precision solutions. The power series method converges rapidly to a machine precision solution, providing both an efficient means for high-accuracy volume rendering, and a reference standard by which other approximations may be measured.