Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/dlaps.doc
There are 2 other files named dlaps.doc in the archive. Click here to see a list.
SUBROUTINE DLAPS

PURPOSE
   COMPUTES THE VALUE OF AN N-TERM EXPANSION IN LAGUERRE
   POLYNOMIALS WITH COEFFICIENT VECTOR C FOR ARGUMENT VALUE X.

USAGE
   CALL DLAPS(Y,X,C,N)

DESCRIPTION OF PARAMETERS
   Y	 - RESULT VALUE
	   DOUBLE PRECISION VARIABLE
   X	 - ARGUMENT VALUE
	   DOUBLE PRECISION VARIABLE
   C	 - COEFFICIENT VECTOR OF GIVEN EXPANSION
	   COEFFICIENTS ARE ORDERED FROM LOW TO HIGH
	   DOUBLE PRECISION VECTOR
   N	 - DIMENSION OF COEFFICIENT VECTOR C

REMARKS
   OPERATION IS BYPASSED IN CASE N LESS THAN 1

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   NONE

METHOD
   DEFINITION
   Y=SUM(C(I)*L(I-1,X), SUMMED OVER I FROM 1 TO N).
   EVALUATION IS DONE BY MEANS OF UPWARD RECURSION
   USING THE RECURRENCE EQUATION FOR LAGUERRE POLYNOMIALS
   L(N+1,X)=2*L(N,X)-L(N-1,X)-((1+X)*L(N,X)-L(N-1,X))/(N+1).