Remove autotools for this branch
Remove make.inc and util.h
QR decomposition version

  - Implemented a new routine, pfitqr, which uses a QR decomposition of a
    Vandermonde matrix to solve the least-squares problem, rather than the
    normal equations. According to Higham (Accuracy and Stability of Numerical
    Algorithms) this should be more numerically stable.

  - The new QR method is the default. The normal equations method is still
    included in this version and can be used via the -n option.

  - The QR method has been tested and gives identical results to the normal
    equations version for my test data. I haven't tested any particularly
    ill-conditioned data, however.
v0.8.1: Fix uninitialized coefficient matrix. Some restructuring.

  - Remove the pfitn routine (which used the normal equations)
  - Separate pfit code into a new file, pfit.f90

  - Was testing this code and discovered that the coefficient matrix
    'coeff' was not getting initialized in the QR pfit
    subroutine. Fixed.

  Notes on testing/verification:

  - Tested d=0 case against analytic solution (agrees)
  - Tested d=1,2,3 against old normal equations version (agrees)
  - Ran valgrind (no errors)
  - Ran with -fbounds-check (no errors)