Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/qhse.doc
There are 2 other files named qhse.doc in the archive. Click here to see a list.
SUBROUTINE QHSE

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

USAGE
   CALL QHSE (H,Y,FDY,SDY,Z,NDIM)

DESCRIPTION OF PARAMETERS
   H	  - THE INCREMENT OF ARGUMENT VALUES.
   Y	  - THE INPUT VECTOR OF FUNCTION VALUES.
   FDY	  - THE INPUT VECTOR OF FIRST DERIVATIVE.
   SDY	  - THE INPUT VECTOR OF SECOND DERIVATIVE.
   Z	  - THE RESULTING VECTOR OF INTEGRAL VALUES. Z MAY BE
	    IDENTICAL WITH Y,FDY OR SDY.
   NDIM   - THE DIMENSION OF VECTORS Y,FDY,SDY,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 SIXTH ORDER INTEGRATION FORMULA.
   FOR REFERENCE, SEE
   R.ZURMUEHL, PRAKTISCHE MATHEMATIK FUER INGENIEURE UND
   PHYSIKER, SPRINGER, BERLIN/GOETTINGEN/HEIDELBERG, 1963,
   PP.227-230.