Google
 

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.