s21ca evaluates the Jacobian elliptic functions sn, cn and dn.

Syntax

C#
public static void s21ca( double u, double m, out double sn, out double cn, out double dn, out int ifail )

Visual Basic
Public Shared Sub s21ca ( _ u As Double, _ m As Double, _ <OutAttribute> ByRef sn As Double, _ <OutAttribute> ByRef cn As Double, _ <OutAttribute> ByRef dn As Double, _ <OutAttribute> ByRef ifail As Integer _ )

Visual Basic

Public Shared Sub s21ca ( _
	u As Double, _
	m As Double, _
	<OutAttribute> ByRef sn As Double, _
	<OutAttribute> ByRef cn As Double, _
	<OutAttribute> ByRef dn As Double, _
	<OutAttribute> ByRef ifail As Integer _
)

Visual C++
public: static void s21ca( double u, double m, [OutAttribute] double% sn, [OutAttribute] double% cn, [OutAttribute] double% dn, [OutAttribute] int% ifail )

F#
static member s21ca : u : float * m : float * sn : float byref * cn : float byref * dn : float byref * ifail : int byref -> unit

Parameters

u

Type: System..::..Double

On entry: the argument

u

and the parameter

m

of the functions, respectively.

Constraints:

$abs (u) \leq \sqrt{λ}$ , where $λ = 1 / x02am$ ;
if $abs (u) < 1 / \sqrt{λ}$ , $abs (m) \leq \sqrt{λ}$ .

m

Type: System..::..Double

On entry: the argument

u

and the parameter

m

of the functions, respectively.

Constraints:

$abs (u) \leq \sqrt{λ}$ , where $λ = 1 / x02am$ ;
if $abs (u) < 1 / \sqrt{λ}$ , $abs (m) \leq \sqrt{λ}$ .

sn: Type: System..::..Double%
On exit: the values of the functions $sn u$ , $cn u$ and $dn u$ , respectively.

cn: Type: System..::..Double%
On exit: the values of the functions $sn u$ , $cn u$ and $dn u$ , respectively.

dn: Type: System..::..Double%
On exit: the values of the functions $sn u$ , $cn u$ and $dn u$ , respectively.

ifail: Type: System..::..Int32%
On exit: $ifail = 0$ unless the method detects an error or a warning has been flagged (see [Error Indicators and Warnings]).

Description

s21ca evaluates the Jacobian elliptic functions of argument

u

and parameter

m

\begin{array}{lcl} sn (u ∣ m) & = & \sin ϕ, \\ cn (u ∣ m) & = & \cos ϕ, \\ dn (u ∣ m) & = & \sqrt{1 - m \sin^{2} ϕ}, \end{array}

where

ϕ

, called the amplitude of

u

, is defined by the integral

u = \int_{0}^{ϕ} \frac{d θ}{\sqrt{1 - m \sin^{2} θ}} .

The elliptic functions are sometimes written simply as

sn u

cn u

and

dn u

, avoiding explicit reference to the parameter

m

Another nine elliptic functions may be computed via the formulae

\begin{array}{lcl} cd u & = & cn u / dn u \\ sd u & = & sn u / dn u \\ nd u & = & 1 / dn u \\ dc u & = & dn u / cn u \\ nc u & = & 1 / cn u \\ sc u & = & sn u / cn u \\ ns u & = & 1 / sn u \\ ds u & = & dn u / sn u \\ cs u & = & cn u / sn u \end{array}

(see Abramowitz and Stegun (1972)).

s21ca is based on a procedure given by Bulirsch (1960), and uses the process of the arithmetic-geometric mean (16.9 in Abramowitz and Stegun (1972)). Constraints are placed on the values of

u

and

m

in order to avoid the possibility of machine overflow.

References

Abramowitz M and Stegun I A (1972) Handbook of Mathematical Functions (3rd Edition) Dover Publications

Bulirsch R (1960) Numerical calculation of elliptic integrals and elliptic functions Numer. Math. 7 76–90

Error Indicators and Warnings

Errors or warnings detected by the method:

$ifail = 1$

On entry,

abs (u) > \sqrt{λ}

, where

λ = 1 / x02am ()

$ifail = 2$

On entry,

abs (m) > \sqrt{λ}

and

abs (u) < 1 / \sqrt{λ}

$ifail = -9000$: An error occured, see message report.

Accuracy

In principle the method is capable of achieving full relative precision in the computed values. However, the accuracy obtainable in practice depends on the accuracy of the standard elementary functions such as SIN and COS.