kronrod_test, a FORTRAN90 code which calls kronrod(), which computes both a Gauss quadrature rule of order N, and the Gauss-Kronrod rule of order 2*N+1.
The computer code and data files made available on this web page are distributed under the MIT license
kronrod, a FORTRAN90 code which computes a Gauss and Gauss-Kronrod pair of quadrature rules of arbitrary order, by Robert Piessens, Maria Branders.