A computer program is described for examination of areas (AUC) and moments of serum concentration vs time data using the Lagrange technique alone or in conjunction with linear or log-trapezoidal methods. The suitability of the Lagrange polynomial approximations to the experimental data during AUC computation is possible through output of interpolated values between consecutive data points. Ill-fitting AUC's between any data points can be replaced with values generated by either trapezoidal method. Along with the partial and total AUC and moment, the program generates model-independent pharmacokinetic parameters such as plasma clearance, terminal slope, half-life, volume of distribution at steady-state, mean residence time, and variance of the residence time.