improvement,
      but, if a polynomial interpolant is desired, a better strategy is to change
      the basis, which is what is done with the Lagrange interpolation method,
      in which case, essentially, the linear system to be solved becomes the
      identity matrix.
    </p>

    <p>
      <b>VANDERMONDE_INTERP_1D</b> needs access to the QR_SOLVE and R8LIB libraries.
      The test code also needs access to the TEST_INTERP library. 
    </p>

    <h3 align = "center">
      Licensing:
    </h3>

    