gorithm.
        </li>
      </ul>
    </p>

    <p>
      The routines, as originally developed, did not include a routine to
      solve linear systems A*x=b once the QR factors had been computed.
      A routine <b>DGBBQRS</b> has been written to carry out this computation
      in cases where the QR factors were computed using <b>DGBBQR2</b>,
      that is, the unblocked code.
    </p>

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

    <p>
      This software library includes routines which are 