Trailing-Edge
-
PDP-10 Archives
-
decuslib10-02
-
43,50145/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.