%N lsqfill %F $MIR/src/subs/polyfit.for
%D Add sums-of-squares to least squares matrix
%: polynomials, calibration, least-squares
%P pjt
%B
        SUBROUTINE lsqfill( count, x, y, pmax )

        INTEGER count, pmax
        REAL    x(*), y(*)
  Inputs:
       count -- the number of data points
       x     -- independent variable
       y     -- dependent variable
       pmax  -- polynomial order of fit
