C348: Elliptic Integral for Complex Argument

Author(s): K.S. Kölbig Library: MATHLIB
Submitter: Submitted: 07.06.1992
Language: Fortran Revised:

Function subprograms CELINT and WELINT calculate, for complex argument z=x+iy and real complementary modulus k' a general elliptic integral of the second kind:

displaymath121

which contains the elliptic integrals of the first and second kind as special cases:

displaymath123

The double-precision version WELINT is available only on computers which support a COMPLEX*16 Fortran data type.

Structure:

FUNCTION subprograms
User Entry Names: CELINT, WELINT
Files Referenced: Unit 6
External References: MTLMTR, ABEND

Usage:

In any arithmetic expression, with tex2html_wrap_inline125 ,

CELINT(Z,AKP,A,B) or WELINT(Z,AKP,A,B) has the value tex2html_wrap_inline127 ,

where CELINT is of type COMPLEX, WELINT is of type COMPLEX*16, Z is of the same type as the function name, and AKP, A, B are of type REAL for CELINT and of type DOUBLE PRECISION for WELINT.

Method:

The evaluation of tex2html_wrap_inline129 is based on the Gauss transformation. For details, in particular for the conformal mapping provided by tex2html_wrap_inline131 , see Ref. 1.

Accuracy:

CELINT (except on CDC and Cray computers) has full single-precision accuracy. For most values of the arguments, WELINT (and CELINT on CDC and Cray computers) has an accuracy of approximately one significant digit less than the machine precision.

Error handling:

Error C348.1: Re tex2html_wrap_inline133 . The function value is set equal to zero, and a message is written on Unit 6, unless subroutine MTLSET (N002) has been called.

Notes:

For other forms of the elliptic integrals see the write-up for RELI1 (C346).

The subprogram is based on the Algol60 procedure elco2 given in Ref. 1.

References:

  1. R. Bulirsch, Numerical calculation of elliptic integrals and elliptic functions, Numer. Math. 7 (1965) 78-90.
tex2html_wrap_inline135

Michel Goossens Tue Jun 4 23:15:57 METDST 1996