Introduction - If you have any usage issues, please Google them yourself
D01AKF IS AN ADAPTIVE INTEGRATOR, ESPECIALLY SUITED TO
C STRONGLY OSCILLATORY NON-SINGULAR INTEGRANDS, WHICH CALCULATES
C AN APPROXIMATION TO THE INTEGRAL OF A FUNCTION OVER A FINITE
C INTERVAL (A,B)
C
C D01AKF ITSELF IS ESSENTIALLY A DUMMY ROUTINE WHOSE FUNCTION IS TO
C PARTITION THE WORK ARRAYS WORK AND IWORK FOR USE BY D01AKV.
C WORK IS PARTITIONED INTO 4 ARRAYS EACH OF SIZE LIMIT, WHERE
C LIMIT = MIN(LWORK/4, LIWORK).
C IWORK IS A SINGLE ARRAY IN D01AKV OF SIZE LIMIT.