SUBROUTINE QHFE

PURPOSE
   TO COMPUTE THE VECTOR OF INTEGRAL VALUES FOR A GIVEN
   EQUIDISTANT TABLE OF FUNCTION AND DERIVATIVE VALUES.

USAGE
   CALL QHFE (H,Y,DERY,Z,NDIM)

DESCRIPTION OF PARAMETERS
   H	  - THE INCREMENT OF ARGUMENT VALUES.
   Y	  - THE INPUT VECTOR OF FUNCTION VALUES.
   DERY   - THE INPUT VECTOR OF DERIVATIVE VALUES.
   Z	  - THE RESULTING VECTOR OF INTEGRAL VALUES. Z MAY BE
	    IDENTICAL WITH Y OR DERY.
   NDIM   - THE DIMENSION OF VECTORS Y,DERY,Z.

REMARKS
   NO ACTION IN CASE NDIM LESS THAN 1.

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   NONE

METHOD
   BEGINNING WITH Z(1)=0, EVALUATION OF VECTOR Z IS DONE BY
   MEANS OF HERMITEAN FOURTH ORDER INTEGRATION FORMULA.
   FOR REFERENCE, SEE
   (1) F.B.HILDEBRAND, INTRODUCTION TO NUMERICAL ANALYSIS,
       MCGRAW-HILL, NEW YORK/TORONTO/LONDON, 1956, PP.314-319.
   (2) R.ZURMUEHL, PRAKTISCHE MATHEMATIK FUER INGENIEURE UND
       PHYSIKER, SPRINGER, BERLIN/GOETTINGEN/HEIDELBERG, 1963,
       PP.227-230.