Trailing-Edge
-
PDP-10 Archives
-
decus_20tap2_198111
-
decus/20-0026/deli1.doc
There are 2 other files named deli1.doc in the archive. Click here to see a list.
SUBROUTINE DELI1
PURPOSE
COMPUTES THE ELLIPTIC INTEGRAL OF FIRST KIND
USAGE
CALL DELI1(RES,X,CK)
DESCRIPTION OF PARAMETERS
RES - RESULT VALUE IN DOUBLE PRECISION
X - UPPER INTEGRATION BOUND (ARGUMENT OF ELLIPTIC
INTEGRAL OF FIRST KIND) IN DOUBLE PRECISION
CK - COMPLEMENTARY MODULUS IN DOUBLE PRECISION
REMARKS
DOUBLE PRECISION MODULUS K = DSQRT(1.D0-CK*CK).
SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
NONE
METHOD
DEFINITION
RES=INTEGRAL(1/SQRT((1+T*T)*(1+(CK*T)**2)), SUMMED
OVER T FROM 0 TO X).
EQUIVALENT ARE THE DEFINITIONS
RES=INTEGRAL(1/(COS(T)*SQRT(1+(CK*TAN(T))**2)), SUMMED
OVER T FROM 0 TO ATAN(X)),
RES=INTEGRAL(1/SQRT(1-(K*SIN(T))**2), SUMMED OVER
T FROM 0 TO ATAN(X)).
EVALUATION
LANDENS TRANSFORMATION IS USED FOR CALCULATION.
REFERENCE
R. BULIRSCH, NUMERICAL CALCULATION OF ELLIPTIC INTEGRALS AND
ELLIPTIC FUNCTIONS.
HANDBOOK SERIES OF SPECIAL FUNCTIONS
NUMERISCHE MATHEMATIK VOL. 7, 1965, PP. 78-90.