%N solve %F $MIR/src/subs/polyfit.for
%D Solve a matrix for lsqfill
%: polynomials, calibration, least-squares
%P pjt
%B
        SUBROUTINE solve( nlsq )

        INTEGER nlsq

  Inputs:
       nlsq -- order of final least squares fit
