Google
 

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

PURPOSE
   COMPUTES THE THREE JACOBIAN ELLIPTIC FUNCTIONS SN, CN, DN.

USAGE
   CALL DJELF(SN,CN,DN,X,SCK)

DESCRIPTION OF PARAMETERS
   SN	 - RESULT VALUE SN(X) IN DOUBLE PRECISION
   CN	 - RESULT VALUE CN(X) IN DOUBLE PRECISION
   DN	 - RESULT VALUE DN(X) IN DOUBLE PRECISION
   X	 - DOUBLE PRECISION ARGUMENT OF JACOBIAN ELLIPTIC
	   FUNCTIONS
   SCK	 - SQUARE OF COMPLEMENTARY MODULUS IN DOUBLE PRECISION

REMARKS
   NONE

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   NONE

METHOD
   DEFINITION
   X=INTEGRAL(1/SQRT((1-T*T)*(1-(K*T)**2)), SUMMED OVER
   T FROM 0 TO SN), WHERE K=SQRT(1-SCK).
   SN*SN + CN*CN = 1
   (K*SN)**2 + DN**2 = 1.
   EVALUATION
   CALCULATION IS DONE USING THE PROCESS OF THE ARITHMETIC
   GEOMETRIC MEAN TOGETHER WITH GAUSS DESCENDING TRANSFORMATION
   BEFORE INVERSION OF THE INTEGRAL TAKES PLACE.
   REFERENCE
   R. BULIRSCH, NUMERICAL CALCULATION OF ELLIPTIC INTEGRALS AND
	  ELLIPTIC FUNCTIOMS.
	  HANDBOOK SERIES OF SPECIAL FUNCTIONS
	  NUMERISCHE MATHEMATIK VOL. 7, 1965, PP. 78-90.