The area under the concentration-time curve (AUC) is a simple but very important measure in pharmacokinetics. The estimate of AUC may be obtained as an immediate result of nonlinear regression analysis, in which a set of observed concentrations are fitted to a certain pharmacokinetic function by estimating its parameters. However, the role of AUC is more significant in cases where the estimate is obtained by direct integration of data with the spirit of noncompartmental pharmacokinetic analysis.1) Despite the importance of AUC in such a situation, there has only been a single class of methods for estimating AUC: piecewise interpolatory methods. This class of methods includes trapezoidal/log-trapezoidal rules, spline fitting, and the so-called Lagrange method, as well as hybrids of each of them. There have already been several reports regarding comparative performances of these methods.2-6) Among the class of methods, Purves reported that "parabolas-through-the-origin then logtrapezoidal rule" (PTTO) performed best, 4) whereas Jawień pointed out the possible drawback of Purves' method in theoretical situations.6) It remains controversial as to which method is optimal in actual settings.In the piecewise interpolatory methods, an interpolating polynomial is constructed for each sub-interval specified by two adjacent data points. Then, the approximate value of AUC is given as a sum of the partial areas bounded by the polynomials. Accuracy of approximation depends on the locations of the sampling points, and several criteria were proposed as part of an optimal sampling strategy. 7,8) In general, as is the case with parameter estimation, 9-12) these criteria were driven so as to minimize the variance of AUC estimates. However, these optimizing strategies seem to be useful in restricted settings, because these strategies require prior knowledge concerning both pharmacokinetic and variance models.Another difficulty involved in the piecewise interpolatory methods is the accuracy of the extrapolated portions of the AUC from the final sampling time to time infinity. To accurately compute the extrapolated portions as completely as possible, the portions are usually calculated from several terminal data points by means of log-/non-linear regressions. However, this operation assumes a regression model and departs to some extent from the framework of noncompartmental analysis. Furthermore, the two operations (quadrature up to the final point and extrapolation) appear to lack theoretical consistency between them. For example, optimizing the locations of sampling times may increase errors in estimated parameters for the extrapolated region, thus increasing the error of the total AUC. In this article, the author proposes a different numerical integration scheme for estimating the AUC over the infinite time interval [0, ∞). The method is based on the Gaussian quadrature.13) Unlike the piecewise interpolatory quadratures, the Gaussian quadrature gains its accuracy by specifying abscissas (sampling times), which in turn can be viewe...