Jacobi elliptic functions

Short description: Mathematical function

In mathematics, the Jacobi elliptic functions are a set of basic elliptic functions. They are found in the description of the motion of a pendulum (see also pendulum (mathematics)), as well as in the design of electronic elliptic filters. While trigonometric functions are defined with reference to a circle, the Jacobi elliptic functions are a generalization which refer to other conic sections, the ellipse in particular. The relation to trigonometric functions is contained in the notation, for example, by the matching notation sn for sin. The Jacobi elliptic functions are used more often in practical problems than the Weierstrass elliptic functions as they do not require notions of complex analysis to be defined and/or understood. They were introduced by Carl Gustav Jakob Jacobi (1829). Carl Friedrich Gauss had already studied special Jacobi elliptic functions in 1797, the lemniscate elliptic functions in particular,[1] but his work was published much later.

## Contents

* 1 Overview
* 2 Notation
* 3 Definition in terms of inverses of elliptic integrals
* 4 Definition as trigonometry: the Jacobi ellipse
* 5 Definition in terms of Jacobi theta functions
* 5.1 Jacobi theta function description
* 5.2 Jacobi zn function
* 6 Comparison between sums and products
* 7 The elliptic nome and its series
* 7.1 Elliptic integral and elliptic nome
* 7.2 Kneser integer sequence
* 7.3 Schellbach Schwarz integer sequence
* 8 Definition in terms of Neville theta functions
* 9 Jacobi transformations
* 9.1 The Jacobi imaginary transformations
* 9.2 The Jacobi real transformations
* 9.3 Other Jacobi transformations
* 10 The Jacobi hyperbola
* 11 Minor functions
* 12 Periodicity, poles, and residues
* 13 Special values
* 13.1 Lemniscatic values
* 13.2 Equianharmonic values
* 13.3 Elementary functions
* 14 Identities
* 14.1 Addition theorems
* 14.2 Half Angle formula
* 14.3 K formulas
* 14.4 Relations between squares of the functions
* 15 Representations of the function values via Theta functions
* 16 Jacobi elliptic functions as solutions of nonlinear ordinary differential equations
* 17 Lambert series expansion in terms of the nome
* 18 Fast computation
* 19 Approximation in terms of hyperbolic functions
* 20 Continued fractions
* 21 Inverse functions
* 22 Map projection
* 23 See also
* 24 Notes
* 25 References
* 26 External links

## Overview

The fundamental rectangle in the complex plane of u

There are twelve Jacobi elliptic functions denoted by pq⁡(u,m), where p and q are any of the letters c, s, n, and d. (Functions of the form pp⁡(u,m) are trivially set to unity for notational completeness.) u is the argument, and m is the parameter, both of which may be complex. In fact, the Jacobi elliptic functions are meromorphic in both u and m.[2] The distribution of the zeros and poles in the u-plane is well-known. However, questions of the distribution of the zeros and poles in the m-plane remain to be investigated.[2]

In the complex plane of the argument u, the twelve functions form a repeating lattice of simple poles and zeroes.[3] Depending on the function, one repeating parallelogram, or unit cell, will have sides of length 2K or 4K on the real axis, and 2K′ or 4K′ on the imaginary axis, where K=K(m) and K′=K(1−m) are known as the quarter periods with K(⋅) being the elliptic integral of the first kind. The nature of the unit cell can be determined by inspecting the "auxiliary rectangle" (generally a parallelogram), which is a rectangle formed by the origin (0,0) at one corner, and (K,K′) as the diagonally opposite corner. As in the diagram, the four corners of the auxiliary rectangle are named s, c, d, and n, going counter-clockwise from the origin. The function pq⁡(u,m) will have a zero at the p corner and a pole at the q corner. The twelve functions correspond to the twelve ways of arranging these poles and zeroes in the corners of the rectangle.

When the argument u and parameter m are real, with 0<m<1, K and K′ will be real and the auxiliary parallelogram will in fact be a rectangle, and the Jacobi elliptic functions will all be real valued on the real line.

Since the Jacobian elliptic functions are doubly periodic in u, they factor through a torus – in effect, their domain can be taken to be a torus, just as cosine and sine are in effect defined on a circle. Instead of having only one circle, we now have the product of two circles, one real and the other imaginary. The complex plane can be replaced by a complex torus. The circumference of the first circle is 4K and the second 4K′, where K and K′ are the quarter periods. Each function has two zeroes and two poles at opposite positions on the torus. Among the points 0, K, K+iK′, iK′ there is one zero and one pole.

The Jacobian elliptic functions are then doubly periodic, meromorphic functions satisfying the following properties:

* There is a simple zero at the corner p, and a simple pole at the corner q.
* The complex number p−q is equal to half the period of the function pq⁡u; that is, the function pq⁡u is periodic in the direction pq, with the period being 2(p−q). The function pq⁡u is also periodic in the other two directions pp′ and pq′, with periods such that p−p′ and p−q′ are quarter periods.

Jacobi elliptic function sn

Jacobi elliptic function cn

Jacobi elliptic function dn

Jacobi elliptic function sc

Plots of four Jacobi Elliptic Functions in the complex plane of u, illustrating their double periodic behavior. Images generated using a version of the domain coloring method.[4] All have values of k=m equal to 0.8.

## Notation

The elliptic functions can be given in a variety of notations, which can make the subject unnecessarily confusing. Elliptic functions are functions of two variables. The first variable might be given in terms of the amplitude φ, or more commonly, in terms of u given below. The second variable might be given in terms of the parameter m, or as the elliptic modulus k, where k2=m, or in terms of the modular angle α, where m=sin2α. The complements of k and m are defined as m′=1−m and k′=m′. These four terms are used below without comment to simplify various expressions.

The twelve Jacobi elliptic functions are generally written as pq⁡(u,m) where p and q are any of the letters c, s, n, and d. Functions of the form pp⁡(u,m) are trivially set to unity for notational completeness. The “major” functions are generally taken to be cn⁡(u,m), sn⁡(u,m) and dn⁡(u,m) from which all other functions can be derived and expressions are often written solely in terms of these three functions, however, various symmetries and generalizations are often most conveniently expressed using the full set. (This notation is due to Gudermann and Glaisher and is not Jacobi's original notation.)

Throughout this article, pq⁡(u,t2)=pq⁡(u;t).

The functions are notationally related to each other by the multiplication rule: (arguments suppressed)

pq⁡⋅p'q'⁡=pq'⁡⋅p'q

from which other commonly used relationships can be derived:

prqr=pq
pr⁡⋅rq⁡=pq
1qp=pq

The multiplication rule follows immediately from the identification of the elliptic functions with the Neville theta functions[5]

pq⁡(u,m)=θp(u,m)θq(u,m)

Also note that:

K(m)=K(k2)=∫01dt(1−t2)(1−mt2)=∫01dt(1−t2)(1−k2t2).

## Definition in terms of inverses of elliptic integrals

Model of the Jacobi amplitude (measured along vertical axis) as a function of independent variables u and the modulus k

There is a definition, relating the elliptic functions to the inverse of the incomplete elliptic integral of the first kind. These functions take the parameters u and m as inputs. The φ that satisfies

u=∫0φdθ1−msin2θ

is called the Jacobi amplitude:

am⁡(u,m)=φ.

In this framework, the elliptic sine sn u (Latin: sinus amplitudinis) is given by

sn⁡(u,m)=sin⁡am⁡(u,m)

and the elliptic cosine cn u (Latin: cosinus amplitudinis) is given by

cn⁡(u,m)=cos⁡am⁡(u,m)

and the delta amplitude dn u (Latin: delta amplitudinis)[note 1]

dn⁡(u,m)=dduam⁡(u,m).

In the above, the value m is a free parameter, usually taken to be real such that 0≤m≤1, and so the elliptic functions can be thought of as being given by two variables, u and the parameter m. The remaining nine elliptic functions are easily built from the above three (sn, cn, dn), and are given in a section below.

In the most general setting, am⁡(u,m) is a multivalued function (in u) with infinitely many logarithmic branch points (the branches differ by integer multiples of 2π), namely the points 2sK(m)+(4t+1)K(1−m)i and 2sK(m)+(4t+3)K(1−m)i where s,t∈ℤ.[6] This multivalued function can be made single-valued by cutting the complex plane along the line segments joining these branch points (the cutting can be done in non-equivalent ways, giving non-equivalent single-valued functions), thus making am⁡(u,m) analytic everywhere except on the branch cuts. In contrast, sin⁡am⁡(u,m) and other elliptic functions have no branch points, give consistent values for every branch of am, and are meromorphic in the whole complex plane. Since every elliptic function is meromorphic in the whole complex plane (by definition), am⁡(u,m) (when considered as a single-valued function) is not an elliptic function.

However, the integral inversion above defines a unique single-valued real-analytic function in a real neighborhood of u=0 if m is real. There is a unique analytic continuation of this function from that neighborhood to u∈ℝ. The analytic continuation of this function is periodic in u if and only if m>1 (with the minimal period 4K(1/m)/m), and it is denoted by am⁡(u,m) in the rest of this article.

Jacobi also introduced the coamplitude function:

coam⁡(u,m)=am⁡(K(m)−u,m).

The Jacobi epsilon function can be defined as[7]

ℰ(u,m)=∫0udn2(t,m)dt

and relates the incomplete elliptic integral of the first kind to the incomplete elliptic integral of the second kind (with parameter m):

E(φ,m)=ℰ(F(φ,m),m).

The Jacobi epsilon function is not an elliptic function. However, unlike the Jacobi amplitude and coamplitude, the Jacobi epsilon function is meromorphic in the whole complex plane (in both u and m).

The Jacobi zn function is defined by

zn⁡(u,m)=∫0u(dn⁡(t,m)2−E(m)K(m))dt.

It is a singly periodic function which is meromorphic in u. Its minimal period is 2K(m). It is related to the Jacobi zeta function by Z(φ,m)=zn⁡(F(φ,m),m).

Note that when φ=π/2, that u then equals the quarter period K.

## Definition as trigonometry: the Jacobi ellipse

Plot of the Jacobi ellipse (x2 + y2/b2 = 1, b real) and the twelve Jacobi elliptic functions pq(u,m) for particular values of angle φ and parameter b. The solid curve is the ellipse, with m = 1 − 1/b2 and u = F(φ,m) where F(*,*) is the elliptic integral of the first kind (with parameter m=k2). The dotted curve is the unit circle. Tangent lines from the circle and ellipse at x = cd crossing the x-axis at dc are shown in light grey.

cos⁡φ,sin⁡φ are defined on the unit circle, with radius r = 1 and angle φ= arc length of the unit circle measured from the positive x-axis. Similarly, Jacobi elliptic functions are defined on the unit ellipse,[citation needed] with a = 1\. Let

x2+y2b2=1,b>1,m=1−1b2,0<m<1,x=rcos⁡φ,y=rsin⁡φ

then:

r(φ,m)=11−msin2φ.

For each angle φ the parameter

u=u(φ,m)=∫0φr(θ,m)dθ

(the incomplete elliptic integral of the first kind) is computed. On the unit circle (a=b=1), u would be an arc length. The quantity u[φ,k]=u(φ,k2) is related to the incomplete elliptic integral of the second kind (with modulus k) by[8]

u[φ,k]=11−k2(1+1−k22E⁡(φ+arctan⁡(1−k2tan⁡φ),1−1−k21+1−k2)−E⁡(φ,k)+k2sin⁡φcos⁡φ21−k2sin2φ),

and therefore is related to the arc length of an ellipse. Let P=(x,y)=(rcos⁡φ,rsin⁡φ) be a point on the ellipse, and let P′=(x′,y′)=(cos⁡φ,sin⁡φ) be the point where the unit circle intersects the line between P and the origin O. Then the familiar relations from the unit circle:

x′=cos⁡φ,y′=sin⁡φ

read for the ellipse:

x′=cn⁡(u,m),y′=sn⁡(u,m).

So the projections of the intersection point P′ of the line OP with the unit circle on the x\- and y-axes are simply cn⁡(u,m) and sn⁡(u,m). These projections may be interpreted as 'definition as trigonometry'. In short:

cn⁡(u,m)=xr(φ,m),sn⁡(u,m)=yr(φ,m),dn⁡(u,m)=1r(φ,m).

For the x and y value of the point P with u and parameter m we get, after inserting the relation:

r(φ,m)=1dn⁡(u,m)

into: x=r(φ,m)cos⁡(φ),y=r(φ,m)sin⁡(φ) that:

x=cn⁡(u,m)dn⁡(u,m),y=sn⁡(u,m)dn⁡(u,m).

The latter relations for the x\- and y-coordinates of points on the unit ellipse may be considered as generalization of the relations x=cos⁡φ,y=sin⁡φ for the coordinates of points on the unit circle.

The following table summarizes the expressions for all Jacobi elliptic functions pq(u,m) in the variables (x,y,r) and (φ,dn) with r=x2+y2

Jacobi elliptic functions pq[u,m] as functions of {x,y,r} and {φ,dn}  | q

|

c  | s  | n  | d

c  | 1 | x/y=cot⁡φ | x/r=cos⁡φ | x=cos⁡φ/dn
s  | y/x=tan⁡φ | 1 | y/r=sin⁡φ | y=sin⁡φ/dn
n  | r/x=sec⁡φ | r/y=csc⁡φ | 1 | r=1/dn
d  | 1/x=sec⁡φdn | 1/y=csc⁡φdn | 1/r=dn | 1

## Definition in terms of Jacobi theta functions

### Jacobi theta function description

Equivalently, Jacobi's elliptic functions can be defined in terms of his theta functions. If we abbreviate ϑ00(0;q) as ϑ00(q), and ϑ01(0;q),ϑ10(0;q),ϑ11(0;q) respectively as ϑ01(q),ϑ10(q),ϑ11(q) (the theta constants) then the theta function elliptic modulus k is k={ϑ10[q(k)]ϑ00[q(k)]}2. We define the nome as q=exp⁡(πiτ) in relation to the period ratio. We have

sn⁡(u;k)=−ϑ00[q(k)]ϑ11[u÷K¯(k);q(k)]ϑ10[q(k)]ϑ01[u÷K¯(k);q(k)]cn⁡(u;k)=ϑ01[q(k)]ϑ10[u÷K¯(k);q(k)]ϑ10[q(k)]ϑ01[u÷K¯(k);q(k)]dn⁡(u;k)=ϑ01[q(k)]ϑ00[u÷K¯(k);q(k)]ϑ00[q(k)]ϑ01[u÷K¯(k);q(k)]

where K¯(k)=2πK(k) and q=exp⁡[−πK′(k)/K(k)].

Edmund Whittaker and George Watson defined the Jacobi theta functions this way in their textbook A Course of Modern Analysis:[9]

ϑ00(v;w)=∏n=1∞(1−w2n)[1+2cos⁡(2v)w2n−1+w4n−2]
ϑ01(v;w)=∏n=1∞(1−w2n)[1−2cos⁡(2v)w2n−1+w4n−2]
ϑ10(v;w)=2w1/4cos⁡v∏n=1∞(1−w2n)[1+2cos⁡(2v)w2n+w4n]
ϑ11(v;w)=−2w1/4sin⁡v∏n=1∞(1−w2n)[1−2cos⁡(2v)w2n+w4n]

### Jacobi zn function

The Jacobi zn function can be expressed by theta functions as well:

zn⁡(u;k)=π2Kϑ0′1[u÷K¯(k);q(k)]ϑ01[u÷K¯(k);q(k)]=π2Kϑ0′0[u÷K¯(k);q(k)]ϑ00[u÷K¯(k);q(k)]+k2sn⁡(u;k)cn⁡(u;k)dn⁡(u;k)=π2Kϑ1′0[u÷K¯(k);q(k)]ϑ10[u÷K¯(k);q(k)]+dn⁡(u;k)sn⁡(u;k)cn⁡(u;k)=π2Kϑ1′1[u÷K¯(k);q(k)]ϑ11[u÷K¯(k);q(k)]−cn⁡(u;k)dn⁡(u;k)sn⁡(u;k)

where  denotes the partial derivative with respect to the left bracket entry:

ϑ'00(v;w)=∂∂vϑ00(v;w)ϑ'01(v;w)=∂∂vϑ01(v;w)

and so on.

The following definition of the Jacobi zn function is identical to the now mentioned formulas:

zn(u;k)=∑n=1∞2πK(k)−1sin⁡[πK(k)−1u]q(k)2n−11−2cos⁡[πK(k)−1u]q(k)2n−1+q(k)4n−2

In a successive way the amplitude sine sn can be generated as follows:

sn⁡(u;k)=2{zn(12u;k)+zn[K(k)−12u;k]}k2+{zn(12u;k)+zn[K(k)−12u;k]}2

## Comparison between sums and products

The reduced elliptic integral of first kind shall be defined as follows again:

K¯(ε)=2πK(ε)

And the reduced elliptic nome shall be defined after this pattern:

q¯(ε)=ε−2q(ε)4

The brothers Peter and Jonathan Borwein also gave these two following formulas for the amplitude sine in their work π and the AGM on page 60 ff:

sn⁡(u;k)=4K¯(k)q¯(k)2sin⁡[u÷K¯(k)]∑n=1∞q(k)n−1[1+q(k)2n−1]1−2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

sn⁡(u;k)=2q¯(k)sin⁡[u÷K¯(k)]∏n=1∞1−2q(k)2ncos⁡[2u÷K¯(k)]+q(k)4n1−2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

This defining formula, which results directly from the inner substitution z->K(k)−z, applies analogously to the cd function:

cd⁡(u;k)=4K¯(k)q¯(k)2cos⁡[u÷K¯(k)]∑n=1∞q(k)n−1[1+q(k)2n−1]1+2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

cd⁡(u;k)=2q¯(k)cos⁡[u÷K¯(k)]∏n=1∞1+2q(k)2ncos⁡[2u÷K¯(k)]+q(k)4n1+2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

These formulas are based on Whittaker and Watson's definition of theta non-zero value functions.

These formulas[10] apply to the cosine amplitude:

cn⁡(u;k)=4K¯(k)q¯(k)2cos⁡[u÷K¯(k)]∑n=1∞(−1)n+1q(k)n−1[1−q(k)2n−1]1−2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

cn⁡(u;k)=21−k24q¯(k)cos⁡[u÷K¯(k)]∏n=1∞1+2q(k)2ncos⁡[2u÷K¯(k)]+q(k)4n1−2q(k)2n−1cos⁡[2u÷K¯(k)]+q(k)4n−2

According to the Whittaker-Watson product formulas, this formula also applies to the delta amplitude function:

dn⁡(u;k)=1−k24∏n=1∞1+2cos⁡[2u÷K¯(k)]q(k)2n−1+q(k)4n−21−2cos⁡[2u÷K¯(k)]q(k)2n−1+q(k)4n−2

With a Hyperbolic secant sum is a definition[11] possible for the Delta Amplitudinis:

dn⁡(z;k)=π2K(1−k2)∑n=−∞∞sech⁡{πK(1−k2)−1[K(k)n+12z]}

## The elliptic nome and its series

### Elliptic integral and elliptic nome

Since the Jacobi functions are defined in terms of the elliptic modulus k(τ), we need to invert this and find τ in terms of k. We start from k′=1−k2, the complementary modulus. As a function of τ it is

k′(τ)=1−k2={ϑ01[q(k)]ϑ00[q(k)]}2

Let us define the elliptic nome and the complete elliptic integral of the first kind:

q(k)=exp⁡[−πK(1−k2)K(k)]

These are two identical definitions of the complete elliptic integral of the first kind:

K(k)=∫0π/211−k2sin⁡(φ)2∂φ

K(k)=π2∑a=0∞[(2a)!]216a(a!)4k2a

An identical definition of the nome function can me produced by using a series. Following function has this identity:

1−1−k241+1−k24=ϑ00[q(k)]−ϑ01[q(k)]ϑ00[q(k)]+ϑ01[q(k)]=[∑n=1∞2q(k)(2n−1)2][1+∑n=1∞2q(k)4n2]−1

Since we may reduce to the case where the imaginary part of τ is greater than or equal to 3/2 (see Modular group), we can assume the absolute value of q is less than or equal to exp⁡(−π3/2)≈0.0658; for values this small the above series converges very rapidly and easily allows us to find the appropriate value for q. By solving this function after q we get this[12][13][14] result:

q(k)=∑n=1∞Sc(n)24n−3(1−1−k241+1−k24)4n−3=k2{12+[∑n=1∞Sc(n+1)24n+1k2n]}4

This table shows numbers of the Schwarz integer sequence A002103 accurately:

Sc(1)  | Sc(2)  | Sc(3)  | Sc(4)  | Sc(5)  | Sc(6)  | Sc(7)  | Sc(8)

|
|
|
|
|
|
|

1  | 2  | 15  | 150  | 1707  | 20910  | 268616  | 3567400

### Kneser integer sequence

The German mathematician Adolf Kneser researched on the integer sequence of the elliptic period ratio in his essay Neue Untersuchung einer Reihe aus der Theorie der elliptischen Funktionen and showed that the generating function of this sequence is an elliptic function. Also a further mathematician with the name Robert Fricke analyzed this integer sequence in his essay Die elliptischen Funktionen und ihre Anwendungen and described the accurate computing methods by using this mentioned sequence. The Kneser integer sequence Kn(n) can be constructed in this way:

Kn(2n)=24n−3(4n2n)+∑m=1n42n−2m(4n2n−2m)Kn(m)

Kn(2n+1)=24n−1(4n+22n+1)+∑m=1n42n−2m+1(4n+22n−2m+1)Kn(m)

Executed examples:

Kn(2)=2×6+1×1=13

Kn(3)=8×20+24×1=184

Kn(4)=32×70+448×1+1×13=2701

Kn(5)=128×252+7680×1+40×13=40456

Kn(6)=512×924+126720×1+1056×13+1×184=613720

Kn(7)=2048×3432+2050048×1+23296×13+56×184=9391936

The Kneser sequence appears in the Taylor series of the period ratio (half period ratio):

14ln⁡(16x2)−πK′(x)4K(x)=∑n=1∞Kn(n)24n−1nx2n

### Schellbach Schwarz integer sequence

The mathematician Karl Heinrich Schellbach discovered the integer number sequence that appears in the MacLaurin series of the Elliptic Nome function. This scientist[15] constructed this sequence A002103 in his work Die Lehre von den elliptischen Integralen und den Thetafunktionen in detail. Especially on page 60 of this work a synthesis route of this sequence is written down in his work. Also the Silesian German mathematician Hermann Amandus Schwarz wrote in his work Formeln und Lehrsätze zum Gebrauche der elliptischen Funktionen in the chapter Berechnung der Grösse k on pages 54 to 56 that integer number sequence down. This Schellbach Schwarz number sequence Sc(n) (OEIS: A002103) was also analyzed by the mathematicians Karl Theodor Wilhelm Weierstrass and Louis Melville Milne-Thomson in the 20th century. The mathematician Adolf Kneser determined a synthesis method for this sequence based on the following pattern:

Sc(n+1)=2n∑m=1nSc(m)Kn(n+1−m)

The Schellbach Schwarz sequence Sc(n) is entered in the online encyclopedia of number sequences under the number A002103 and the Kneser sequence Kn(n) is entered under the number A227503. Adolf Kneser researched on this integer sequence in his essay Neue Untersuchung einer Reihe aus der Theorie der elliptischen Funktionen and showed that the generating function of this sequence is an elliptic function. Also Robert Fricke analyzed this integer sequence in his essay Die elliptischen Funktionen und ihre Anwendungen and described accurate computing methods by using this sequence. Following table[16][17] contains the Kneser numbers and the Schellbach Schwarz numbers:

Constructed sequences Kneser and Schellbach Schwarz  Index n  | Kn(n) (A227503)  | Sc(n) (A002103)

|
|

1  | 1  | 1
2  | 13  | 2
3  | 184  | 15
4  | 2701  | 150
5  | 40456  | 1707
6  | 613720  | 20910
7  | 9391936  | 268616
8  | 144644749  | 3567400

In the following, it will be shown as an example how the Schellbach Schwarz numbers are built up successively. For this, the examples with the numbers Sc(4) = 150, Sc(5) = 1707 and Sc(6) = 20910 are used:

Sc(4)=23∑m=13Sc(m)Kn(4−m)=23[Sc(1)Kn(3)+Sc(2)Kn(2)+Sc(3)Kn(1)]

Sc(4)=23(1×184+2×13+15×1)=150

Sc(5)=24∑m=14Sc(m)Kn(5−m)=24[Sc(1)Kn(4)+Sc(2)Kn(3)+Sc(3)Kn(2)+Sc(4)Kn(1)]

Sc(5)=24(1×2701+2×184+15×13+150×1)=1707

Sc(6)=25∑m=15Sc(m)Kn(6−m)=25[Sc(1)Kn(5)+Sc(2)Kn(4)+Sc(3)Kn(3)+Sc(4)Kn(2)+Sc(5)Kn(1)]

Sc(6)=25(1×40456+2×2701+15×184+150×13+1707×1)=20910

And this sequence creates the MacLaurin series of the elliptic nome in exactly this way mentioned above:

q(x)=∑n=1∞Sc(n)24n−3(1−1−x241+1−x24)4n−3=x2{12+[∑n=1∞Sc(n+1)24n+1x2n]}4

## Definition in terms of Neville theta functions

The Jacobi elliptic functions can be defined very simply using the Neville theta functions:[18]

pq⁡(u,m)=θp(u,m)θq(u,m)

Simplifications of complicated products of the Jacobi elliptic functions are often made easier using these identities.

## Jacobi transformations

### The Jacobi imaginary transformations

Plot of the degenerate Jacobi curve (x2 + y2/b2 = 1, b = ∞) and the twelve Jacobi Elliptic functions pq(u,1) for a particular value of angle φ. The solid curve is the degenerate ellipse (x2 = 1) with m = 1 and u = F(φ,1) where F(*,*) is the elliptic integral of the first kind. The dotted curve is the unit circle. Since these are the Jacobi functions for m = 0 (circular trigonometric functions) but with imaginary arguments, they correspond to the six hyperbolic trigonometric functions.

The Jacobi imaginary transformations relate various functions of the imaginary variable i u or, equivalently, relations between various values of the m parameter. In terms of the major functions:[19]: 506

cn⁡(u,m)=nc⁡(iu,1−m)
sn⁡(u,m)=−isc⁡(iu,1−m)
dn⁡(u,m)=dc⁡(iu,1−m)

Using the multiplication rule, all other functions may be expressed in terms of the above three. The transformations may be generally written as pq⁡(u,m)=γpqpq′⁡(iu,1−m). The following table gives the γpqpq′⁡(iu,1−m) for the specified pq(u,m).[18] (The arguments (iu,1−m) are suppressed)

Jacobi Imaginary transformations γpqpq′⁡(iu,1−m) | q

|

c  | s  | n  | d

c  | 1 | i ns | nc | nd
s  | −i sn | 1 | −i sc | −i sd
n  | cn | i cs | 1 | cd
d  | dn | i ds | dc | 1

Since the hyperbolic trigonometric functions are proportional to the circular trigonometric functions with imaginary arguments, it follows that the Jacobi functions will yield the hyperbolic functions for m=1.[5]: 249 In the figure, the Jacobi curve has degenerated to two vertical lines at x = 1 and x = −1.

### The Jacobi real transformations

The Jacobi real transformations[5]: 308 yield expressions for the elliptic functions in terms with alternate values of m. The transformations may be generally written as pq⁡(u,m)=γpqpq′⁡(ku,1/m). The following table gives the γpqpq′⁡(ku,1/m) for the specified pq(u,m).[18] (The arguments (ku,1/m) are suppressed)

Jacobi real transformations γpqpq′⁡(ku,1/m) | q

|

c  | s  | n  | d

c  | 1 | kds | dn | dc
s  | 1ksd | 1 | 1ksn | 1ksc
n  | nd | kns | 1 | nc
d  | cd | kcs | cn | 1

### Other Jacobi transformations

Jacobi's real and imaginary transformations can be combined in various ways to yield three more simple transformations .[5]: 214 The real and imaginary transformations are two transformations in a group (D3 or anharmonic group) of six transformations. If

μR(m)=1/m

is the transformation for the m parameter in the real transformation, and

μI(m)=1−m=m′

is the transformation of m in the imaginary transformation, then the other transformations can be built up by successive application of these two basic transformations, yielding only three more possibilities:

μIR(m)=μI(μR(m))=−m′/mμRI(m)=μR(μI(m))=1/m′μRIR(m)=μR(μI(μR(m)))=−m/m′

These five transformations, along with the identity transformation (μU(m) = m) yield the six-element group. With regard to the Jacobi elliptic functions, the general transformation can be expressed using just three functions:

cs⁡(u,m)=γics'⁡(γiu,μi(m))
ns⁡(u,m)=γins'⁡(γiu,μi(m))
ds⁡(u,m)=γids'⁡(γiu,μi(m))

where i = U, I, IR, R, RI, or RIR, identifying the transformation, γi is a multiplication factor common to these three functions, and the prime indicates the transformed function. The other nine transformed functions can be built up from the above three. The reason the cs, ns, ds functions were chosen to represent the transformation is that the other functions will be ratios of these three (except for their inverses) and the multiplication factors will cancel.

The following table lists the multiplication factors for the three ps functions, the transformed m's, and the transformed function names for each of the six transformations.[5]: 214 (As usual, k2 = m, 1 − k2 = k12 = m′ and the arguments (γiu,μi(m)) are suppressed)

Parameters for the six transformations  Transformation i | γi | μi(m) | cs' | ns' | ds'

|
|
|
|
|

U  | 1 | m | cs | ns | ds
I  | i | m′ | ns | cs | ds
IR  | i k | −m′/m | ds | cs | ns
R  | k | 1/m | ds | ns | cs
RI  | i k1 | 1/m′ | ns | ds | cs
RIR  | k1 | −m/m′ | cs | ds | ns

Thus, for example, we may build the following table for the RIR transformation.[18] The transformation is generally written pq⁡(u,m)=γpqpq'⁡(k′u,−m/m′) (The arguments (k′u,−m/m′) are suppressed)

The RIR transformation γpqpq'⁡(k′u,−m/m′) | q

|

c  | s  | n  | d

c  | 1 | k' cs | cd | cn
s  | 1k′ sc | 1 | 1k′ sd | 1k′ sn
n  | dc | k′ ds | 1 | dn
d  | nc | k′ ns | nd | 1

The value of the Jacobi transformations is that any set of Jacobi elliptic functions with any complex-valued parameter m can be converted into another set for which 0 ≤ m ≤ 1 and, for real values of u, the function values will be real.[5]: p. 215

## The Jacobi hyperbola

Plot of the Jacobi hyperbola (x2 + y2/b2 = 1, b imaginary) and the twelve Jacobi Elliptic functions pq(u,m) for particular values of angle φ and parameter b. The solid curve is the hyperbola, with m = 1 − 1/b2 and u = F(φ,m) where F(*,*) is the elliptic integral of the first kind. The dotted curve is the unit circle. For the ds-dc triangle, σ = sinφ cos φ.

Introducing complex numbers, our ellipse has an associated hyperbola:

x2−y2b2=1

from applying Jacobi's imaginary transformation[18] to the elliptic functions in the above equation for x and y.

x=1dn⁡(u,1−m),y=sn⁡(u,1−m)dn⁡(u,1−m)

It follows that we can put x=dn⁡(u,1−m),y=sn⁡(u,1−m). So our ellipse has a dual ellipse with m replaced by 1-m. This leads to the complex torus mentioned in the Introduction.[20] Generally, m may be a complex number, but when m is real and m<0, the curve is an ellipse with major axis in the x direction. At m=0 the curve is a circle, and for 0<m<1, the curve is an ellipse with major axis in the y direction. At m = 1, the curve degenerates into two vertical lines at x = ±1. For m > 1, the curve is a hyperbola. When m is complex but not real, x or y or both are complex and the curve cannot be described on a real x-y diagram.

## Minor functions

Reversing the order of the two letters of the function name results in the reciprocals of the three functions above:

ns⁡(u)=1sn⁡(u),nc⁡(u)=1cn⁡(u),nd⁡(u)=1dn⁡(u).

Similarly, the ratios of the three primary functions correspond to the first letter of the numerator followed by the first letter of the denominator:

sc⁡(u)=sn⁡(u)cn⁡(u),sd⁡(u)=sn⁡(u)dn⁡(u),dc⁡(u)=dn⁡(u)cn⁡(u),ds⁡(u)=dn⁡(u)sn⁡(u),cs⁡(u)=cn⁡(u)sn⁡(u),cd⁡(u)=cn⁡(u)dn⁡(u).

More compactly, we have

pq⁡(u)=pn⁡(u)qn⁡(u)

where p and q are any of the letters s, c, d.

## Periodicity, poles, and residues

Plots of the phase for the twelve Jacobi Elliptic functions pq(u,m) as a function complex argument u, with poles and zeroes indicated. The plots are over one full cycle in the real and imaginary directions with the colored portion indicating phase according to the color wheel at the lower right (which replaces the trivial dd function). Regions with absolute value below 1/3 are colored black, roughly indicating the location of a zero, while regions with absolute value above 3 are colored white, roughly indicating the position of a pole. All plots use m = 2/3 with K = K(m), K′ = K(1 − m), K(*) being the complete elliptic integral of the first kind. Arrows at the poles point in direction of zero phase. Right and left arrows imply positive and negative real residues respectively. Up and down arrows imply positive and negative imaginary residues respectively.

In the complex plane of the argument u, the Jacobi elliptic functions form a repeating pattern of poles (and zeroes). The residues of the poles all have the same absolute value, differing only in sign. Each function pq(u,m) has an "inverse function" (in the multiplicative sense) qp(u,m) in which the positions of the poles and zeroes are exchanged. The periods of repetition are generally different in the real and imaginary directions, hence the use of the term "doubly periodic" to describe them.

The Jacobi amplitude and the Jacobi epsilon function are quasi-periodic:

am⁡(u+2K,m)=am⁡(u,m)+π,
ℰ(u+2K,m)=ℰ(u,m)+2E

where E(m) is the complete elliptic integral of the second kind with parameter m.

Also

zn⁡(u+2iK′,m)=zn⁡(u,m)−πiK.

The double periodicity of the Jacobi elliptic functions may be expressed as:

pq⁡(u+2αK(m)+2iβK(1−m),m)=(−1)γpq⁡(u,m)

where α and β are any pair of integers. K(*) is the complete elliptic integral of the first kind, also known as the quarter period. The power of negative unity (γ) is given in the following table:

γ | q

|

c  | s  | n  | d

c  | 0 | β | α + β | α
s  | β | 0 | α | α + β
n  | α + β | α | 0 | β
d  | α | α + β | β | 0

When the factor (−1)γ is equal to −1, the equation expresses quasi-periodicity. When it is equal to unity, it expresses full periodicity. It can be seen, for example, that for the entries containing only α when α is even, full periodicity is expressed by the above equation, and the function has full periods of 4K(m) and 2iK(1 − m). Likewise, functions with entries containing only β have full periods of 2K(m) and 4iK(1 − m), while those with α \+ β have full periods of 4K(m) and 4iK(1 − m).

In the diagram on the right, which plots one repeating unit for each function, indicating phase along with the location of poles and zeroes, a number of regularities can be noted: The inverse of each function is opposite the diagonal, and has the same size unit cell, with poles and zeroes exchanged. The pole and zero arrangement in the auxiliary rectangle formed by (0,0), (K,0), (0,K′) and (K,K′) are in accordance with the description of the pole and zero placement described in the introduction above. Also, the size of the white ovals indicating poles are a rough measure of the absolute value of the residue for that pole. The residues of the poles closest to the origin in the figure (i.e. in the auxiliary rectangle) are listed in the following table:

Residues of Jacobi elliptic functions  | q

|

c  | s  | n  | d

c  |  | 1 | −ik | −1k
s  | −1k′ |  | 1k | −ikk′
n  | −1k′ | 1 |  | −ik′
d  | -1 | 1 | −i |

When applicable, poles displaced above by 2K or displaced to the right by 2K′ have the same value but with signs reversed, while those diagonally opposite have the same value. Note that poles and zeroes on the left and lower edges are considered part of the unit cell, while those on the upper and right edges are not.

## Special values

### Lemniscatic values

The Values along with the modulus k=122=λ*(1) are called lemniscatic values:

Values for the thirds of the integral K:

sn⁡[13K(122);122]=12(124−3+1)

cn⁡[23K(122);122]=12(3−124+1)

cn⁡[13K(122);122]=23−34

sn⁡[23K(122);122]=12(274+34−2)23−34

Values for the fifths of the integral K:

sn⁡[15K(122);122]=12(tan⁡15π+1−2tan⁡320π)

sn⁡[35K(122);122]=12(tan⁡15π−1+2cot⁡320π)

cn⁡[25K(122);122]=12(cot⁡110π−1−2tan⁡120π)

cn⁡[45K(122);122]=12(cot⁡110π+1−2cot⁡120π)

For the corresponding Pythagorean opposites these formulas are valid:

cn⁡[15K(122);122]=25−24sin⁡320πcos⁡120π

cn⁡[35K(122);122]=25−24sin⁡120πcos⁡320π

sn⁡[25K(122);122]=122tan⁡120π4(cot⁡110π+1−cot⁡110π−3)

sn⁡[45K(122);122]=122cot⁡120π4(cot⁡110π+3−cot⁡110π−1)

Values for the sevenths of the integral K:

cn⁡[27K(122);122]=tanh⁡[12arcoth⁡(122sin⁡17πcot⁡328π+sin⁡114π)]

cn⁡[47K(122);122]=tanh⁡[12arcoth⁡(122cos⁡114πtan⁡528π+sin⁡314π)]

cn⁡[67K(122);122]=tanh⁡[12arcoth⁡(122cos⁡314πcot⁡128π+cos⁡17π)]

For the corresponding Pythagorean opposites these formulas are valid:

sn⁡[27K(122);122]=sech⁡[12arcoth⁡(122sin⁡17πcot⁡328π+sin⁡114π)]

sn⁡[47K(122);122]=sech⁡[12arcoth⁡(122cos⁡114πtan⁡528π+sin⁡314π)]

sn⁡[67K(122);122]=sech⁡[12arcoth⁡(122cos⁡314πcot⁡128π+cos⁡17π)]

Important identities:

Setting k=122 gives the lemniscate elliptic functions sl and cl in this way:

sl⁡(u)=122sd⁡(2u;122)=122sn⁡(2u;122)dn⁡(2u;122),cl⁡(u)=cn⁡(2u;122).

Setting k=i gives the lemniscate elliptic functions sl and cl in the following way:

sl⁡(u)=sn⁡(u;i),cl⁡(u)=cd⁡(u;i)=cn⁡(u;i)dn⁡(u;i).

### Equianharmonic values

The Values along with the modulus k=sin⁡112π=λ*(3) and k=cos⁡112π=λ*(13) are called equianharmonic values or also values of the Equianharmonic Case:

cn⁡[25K(sin⁡112π);sin⁡112π]=tan⁡[arctan⁡(103tan⁡110π+133(103−1))−112π]

cn⁡[45K(sin⁡112π);sin⁡112π]=tan⁡[arctan⁡(103tan⁡110π−133(103−1))+112π]

cn⁡[25K(cos⁡112π);cos⁡112π]=cot⁡[arctan⁡(103cot⁡15π+133(103−1))−112π]

cn⁡[45K(cos⁡112π);cos⁡112π]=cot⁡[arctan⁡(103cot⁡15π−133(103−1))+112π]

The flower marks the Elliptic lambda star function!

### Elementary functions

When k=0 or k=1, the Jacobi elliptic functions are reduced to non-elliptic functions:

Function  | k = 0  | k = 1

|
|

sn⁡(u;k) | sin⁡u | tanh⁡u
cn⁡(u;k) | cos⁡u | sech⁡u
dn⁡(u;k) | 1 | sech⁡u
ns⁡(u;k) | csc⁡u | coth⁡u
nc⁡(u;k) | sec⁡u | cosh⁡u
nd⁡(u;k) | 1 | cosh⁡u
sd⁡(u;k) | sin⁡u | sinh⁡u
cd⁡(u;k) | cos⁡u | 1
cs⁡(u;k) | cot⁡u | csch⁡u
ds⁡(u;k) | csc⁡u | csch⁡u
dc⁡(u;k) | sec⁡u | 1
sc⁡(u;k) | tan⁡u | sinh⁡u

For the Jacobi amplitude, am⁡(u;0)=u and am⁡(u;1)=gd⁡u where gd is the Gudermannian function.

In general if neither of p, q is d then pq⁡(u;1)=pq⁡(gd⁡(u);0).

## Identities

### Addition theorems

The functions satisfy the two square relations in dependency of the modulus k in Legendre form: cn2(u;k)+sn2(u;k)=1,

dn2(u;k)+k2sn2(u;k)=1.

From this we see that (cn, sn, dn) parametrizes an elliptic curve which is the intersection of the two quadrics defined by the above two equations. We now may define a group law for points on this curve by the addition formulas for the Jacobi functions[3][21]

sn⁡(x+y;k)=sn⁡(x;k)cn⁡(y;k)dn⁡(y;k)+sn⁡(y;k)cn⁡(x;k)dn⁡(x;k)1−k2sn2(x;k)sn2(y;k),cn⁡(x+y;k)=cn⁡(x;k)cn⁡(y;k)−sn⁡(x;k)sn⁡(y;k)dn⁡(x;k)dn⁡(y;k)1−k2sn2(x;k)sn2(y;k),dn⁡(x+y;k)=dn⁡(x;k)dn⁡(y;k)−k2sn⁡(x;k)sn⁡(y;k)cn⁡(x;k)cn⁡(y;k)1−k2sn2(x;k)sn2(y;k).

By only treating the functions sn and cd, the following pair of addition theorems can be established in which the two theorems are antisymmetric to each other:

sn(x+y;k)=sn(x;k)cd(y;k)+cd(x;k)sn(y;k)1+k2sn(x;k)cd(x;k)sn(y;k)cd(y;k)

cd(x+y;k)=cd(x;k)cd(y;k)−sn(x;k)sn(y;k)1−k2sn(x;k)cd(x;k)sn(y;k)cd(y;k)

The Jacobi epsilon and zn functions satisfy a quasi-addition theorem: ℰ(x+y;k)=ℰ(x;k)+ℰ(y;k)−k2sn⁡(x;k)sn⁡(y;k)sn⁡(x+y;k),zn⁡(x+y;k)=zn⁡(x;k)+zn⁡(y;k)−k2sn⁡(x;k)sn⁡(y;k)sn⁡(x+y;k).

Double angle formulae can be easily derived from the above equations by setting x = y.[3] Half angle formulae[18][3] are all of the form:

pq⁡(12u;k)2=fp/fq

where: fc=cn⁡(u;k)+dn⁡(u;k) fs=1−cn⁡(u;k) fn=1+dn⁡(u;k) fd=(1+dn⁡(u;k))−m(1−cn⁡(u;k))

### Half Angle formula

These three formulas describe the halving theorem:

sn⁡(12u;k)2=1−cn⁡(u;k)1+dn⁡(u;k) cn⁡(12u;k)2=cn⁡(u;k)+dn⁡(u;k)1+dn⁡(u;k) dn⁡(12u;k)2=1−k2+dn⁡(u;k)+k2cn⁡(u;k)1+dn⁡(u;k)

And the theorem of the arithmetic mean is described by this formula:

sn⁡(12x+12y;k)2=1+sn⁡(x;k)sn⁡(y;k)−cn⁡(x;k)cn⁡(y;k)1+k2sn⁡(x;k)sn⁡(y;k)+dn⁡(x;k)dn⁡(y;k)

### K formulas

Half K formula

sn⁡[12K(k);k]=21+k+1−k

cn⁡[12K(k);k]=21−k241+k+1−k

dn⁡[12K(k);k]=1−k24

Third K formula

Also this equation[21] leads to the sn-value of the third of K:

k2s4−2k2s3+2s−1=0

s=sn⁡[13K(k);k]

These equations lead to the other values of the Jacobi functions:

cn⁡[23K(k);k]=1−sn⁡[13K(k);k]

dn⁡[23K(k);k]=1/sn⁡[13K(k);k]−1

On the basis of that now mentioned quartic equation we can a simplified formula by using the solution algorithm for the General case of the quartic equation. This parametrized formula can be generated in that way:

sn⁡{13K[tan⁡(12arctan⁡x3)];tan⁡(12arctan⁡x3)}=tanh⁡[12ln⁡(2x4−x2+1−x2+2+x2+1)]

The Pythagorean opposite gives that formula successively:

cn⁡{13K[tan⁡(12arctan⁡x3);tan⁡(12arctan⁡x3)}=sech⁡[12ln⁡(2x4−x2+1−x2+2+x2+1)]

These two now mentioned formulas are valid for all real values of x, the criterion x∈ℝ is fulfilled in both equations.

To get x, we take the tangent of twice the arctangent of the modulus and then we take the cube root so that x appears.

And then we insert the generated value x into the right side of the balance of the shown parametrized equations.

In relation to this calculation algorithm, three examples are given below:

First example:

122=tan⁡[12arctan⁡(x3)](x=2)

If we insert the value x=2 into the fraction showed above, we get this result:

sn⁡[13K(122);122]=tanh⁡[12ln⁡(124+3)]=12(124−3+1)

cn⁡[13K(122);122]=sech⁡[12ln⁡(124+3)]=23−34

Second example:

12=tan⁡[12arctan⁡(x3)](x=4/33)

If we insert the value x=433 into the fraction showed above, we get this result:

sn⁡[13K(12);12]=tanh⁡[12ln⁡(2(43)4/3−(43)2/3+1−(43)2/3+2+(43)2/3+1)]

cn⁡[13K(12);12]=sech⁡[12ln⁡(2(43)4/3−(43)2/3+1−(43)2/3+2+(43)2/3+1)]

Third example:

14=tan⁡[12arctan⁡(x3)](x=8/153)

If we insert the value x=8153 into the fraction showed above, we get this result:

sn⁡[13K(14);14]=tanh⁡[12ln⁡(2(815)4/3−(815)2/3+1−(815)2/3+2+(815)2/3+1)]

cn⁡[13K(14);14]=sech⁡[12ln⁡(2(815)4/3−(815)2/3+1−(815)2/3+2+(815)2/3+1)]

Fifth K formula

Following equation has following solution:

4k2t6+8k2t5+2(1−k2)2t−(1−k2)2=0

t=12−12k2sn⁡[25K(k);k]2sn⁡[45K(k);k]2=sn⁡[45K(k);k]2−sn⁡[25K(k);k]22sn⁡[25K(k);k]sn⁡[45K(k);k]

To get the sn-values, we put the solution x into following expressions:

sn⁡[25K(k);k]=(1+k2)−1/22(1−t−t2)(t2+1−tt2+1)

sn⁡[45K(k);k]=(1+k2)−1/22(1−t−t2)(t2+1+tt2+1)

The t solution between zero and a half is this solution which has to be used for these expressions.

In relation to the now mentioned calculation algorithm, an important example are given below:

This equation of sixth degree is given:

4k2t6+8k2t5+2(1−k2)2t−(1−k2)2=0

The value k=122 shall be entered in this example.

By entering this modulus value and multiplying the whole equation with the factor 4 we get that equation:

8t6+16t5+2t−1=0

This is the solution of that equation for t that is located between zero and a half:

t=14(5+1)(54−1)

This solution t and the corresponding modulus k will be inserted into the formulas for the amplitude sine values:

sn⁡[25K(122);122]=(1+k2)−1/22(1−t−t2)(t2+1−tt2+1)[k=122∩t=14(5+1)(54−1)]

sn⁡[45K(122);122]=(1+k2)−1/22(1−t−t2)(t2+1+tt2+1)[k=122∩t=14(5+1)(54−1)]

By doing this, following results appear, that are mentioned already in the section of the special values:

sn⁡[25K(122);122]=122tan⁡120π4(cot⁡110π+1−cot⁡110π−3)

sn⁡[45K(122);122]=122cot⁡120π4(cot⁡110π+3−cot⁡110π−1)

By using the hyperbolic functions, these equivalent expressions can be produced:

sn⁡[25K(122);122]=sech⁡{12arcoth⁡[12cot⁡110π−12]}

sn⁡[45K(122);122]=sech⁡{12arcoth⁡[12cot⁡110π+12]}

According to Abel–Ruffini theorem, it is completely impossible to set up elementary expressions for the amplitude sine values of the fifths of the complete first kind elliptic integral for the regular case of the elliptic modulus.

### Relations between squares of the functions

Relations between squares of the functions can be derived from two basic relationships (Arguments (u,m) suppressed): cn2+sn2=1 cn2+m′sn2=dn2 where m + m' = 1\. Multiplying by any function of the form nq yields more general equations:

cq2+sq2=nq2 cq2+m′sq2=dq2

With q = d, these correspond trigonometrically to the equations for the unit circle (x2+y2=r2) and the unit ellipse (x2+m′y2=1), with x = cd, y = sd and r = nd. Using the multiplication rule, other relationships may be derived. For example:

−dn2+m′=−mcn2=msn2−m

−m′nd2+m′=−mm′sd2=mcd2−m

m′sc2+m′=m′nc2=dc2−m

cs2+m′=ds2=ns2−m

## Representations of the function values via Theta functions

Given are the product definitions of the Jacobi Theta Nullwert Functions as the mathematicians Edmund Taylor Whittaker and George Neville Watson[22][23][24] set them up:

ϑ00(w)=∏n=1∞(1−w2n)(1+w2n−1)2
ϑ01(w)=∏n=1∞(1−w2n)(1−w2n−1)2

Identities of thirds of the integral K:

With the so-called theta zero-value functions of the elliptic noun of the module, many Jacobian function values can be represented:

sc⁡[23K(k);k]=3ϑ01[q(k)6]1−k2ϑ01[q(k)2]
sn⁡[13K(k);k]=2ϑ00[q(k)]23ϑ00[q(k)3]2+ϑ00[q(k)]2=3ϑ01[q(k)3]2−ϑ01[q(k)]23ϑ01[q(k)3]2+ϑ01[q(k)]2
cn⁡[23K(k);k]=3ϑ00[q(k)3]2−ϑ00[q(k)]23ϑ00[q(k)3]2+ϑ00[q(k)]2=2ϑ01[q(k)]23ϑ01[q(k)3]2+ϑ01[q(k)]2

Identities of the fifths of the integral K:

sn⁡[15K(k);k]={5ϑ01[q(k)5]ϑ01[q(k)]−1}{5ϑ01[q(k)10]2ϑ01[q(k)2]2−1}−1
sn⁡[35K(k);k]={5ϑ01[q(k)5]ϑ01[q(k)]+1}{5ϑ01[q(k)10]2ϑ01[q(k)2]2−1}−1
cn⁡[25K(k);k]={5ϑ00[q(k)5]ϑ00[q(k)]+1}{5ϑ01[q(k)10]2ϑ01[q(k)2]2−1}−1
cn⁡[45K(k);k]={5ϑ00[q(k)5]ϑ00[q(k)]−1}{5ϑ01[q(k)10]2ϑ01[q(k)2]2−1}−1

The elementary combinations of theta zero-value functions and the elliptic nome are not sufficient for the representation of the Jacobian function values of left-hand parenthesis entries beyond rationally broken K-integrals. This requires the theta non-zero functions of the pattern described above.

## Jacobi elliptic functions as solutions of nonlinear ordinary differential equations

The derivatives of the three basic Jacobi elliptic functions are:

ddzsn⁡(z)=cn⁡(z)dn⁡(z), ddzcn⁡(z)=−sn⁡(z)dn⁡(z), ddzdn⁡(z)=−msn⁡(z)cn⁡(z).

These can be used to derive the derivatives of all other functions as shown in the table below (arguments (u,m) suppressed):

Derivatives ddupq⁡(u,m) | q

|

c  | s  | n  | d

c  | 0 | −ds ns | −dn sn | −m' nd sd
s  | dc nc | 0 | cn dn | cd nd
n  | dc sc | −cs ds | 0 | m cd sd
d  | m' nc sc | −cs ns | −m cn sn | 0

With the addition theorems above and for a given m with 0 < m < 1 the major functions are therefore solutions to the following nonlinear ordinary differential equations:

* sn⁡(x) solves the differential equations d2ydx2+(1+m)y−2my3=0 and (dydx)2=(1−y2)(1−my2)
* cn⁡(x) solves the differential equations d2ydx2+(1−2m)y+2my3=0 and (dydx)2=(1−y2)(1−m+my2)
* dn⁡(x) solves the differential equations d2ydx2−(2−m)y+2y3=0 and (dydx)2=(y2−1)(1−m−y2)

The Jacobi amplitude provides a non-trivial solution of the differential equation describing the exact motion of a simple pendulum. In particular,

θ=2am⁡(t2c2,2)->d2θdt2+csin⁡θ=0.

A function which solves the above pendulum differential equation with initial angle θ0 is

θ=2arcsin⁡(kcd⁡(ct,k2)),k=sin⁡θ02.

## Lambert series expansion in terms of the nome

Let the nome be q=exp⁡(−πK′(m)/K(m))=eiπτ, Im⁡(τ)>0, m=k2 and let v=πu/(2K(m)). Then the functions have expansions as Lambert series

sn⁡(u,m)=2πkK(m)∑n=0∞qn+1/21−q2n+1sin⁡((2n+1)v),

cn⁡(u,m)=2πkK(m)∑n=0∞qn+1/21+q2n+1cos⁡((2n+1)v),

dn⁡(u,m)=π2K(m)+2πK(m)∑n=1∞qn1+q2ncos⁡(2nv)

zn⁡(u,m)=2πK(m)∑n=1∞qn1−q2nsin⁡(2nv)

when qexp⁡(2|Im⁡(v)|)<1.

For the Jacobi amplitude,

am⁡(u,m)=πu2K(m)+2∑n=1∞qnn(1+q2n)sin⁡(2nv)

where 0<m<1 and u∈ℝ.

Bivariate power series expansions have been published by Schett.[25]

## Fast computation

The theta function ratios provide an efficient way of computing the Jacobi elliptic functions. There is an alternative method, based on the arithmetic-geometric mean and Landen's transformations:[6]

Initialize

a0=1,b0=1−m,c0=1−b02

where 0<m<1. Define

an=an−1+bn−12,bn=an−1bn−1,cn=an−1−bn−12

where n≥1. Then define

φN=2NaNu

for u∈ℝ and a fixed N∈ℕ. If

φn−1=12(φn+arcsin⁡(cnansin⁡φn))

for n≥1, then

am⁡(u,m)=φ0

as N->∞. This is notable for its rapid convergence. It is then trivial to compute all Jacobi elliptic functions from the Jacobi amplitude am on the real line.

In conjunction with the addition theorems for elliptic functions (which hold for complex numbers in general) and the Jacobi transformations, the method of computation described above can be used to compute all Jacobi elliptic functions in the whole complex plane.

Another method of fast computation of the Jacobi elliptic functions via the arithmetic–geometric mean, avoiding the computation of the Jacobi amplitude, is due to Herbert E. Salzer:[26]

Let

0≤m≤1,0≤u≤K(m),a0=1,b0=1−m,
an+1=an+bn2,bn+1=anbn,cn+1=an−bn2.

Set

yN=aNsin⁡(aNu)yN−1=yN+aNcNyNyN−2=yN−1+aN−1cN−1yN−1⋮y0=y1+m4y1.

Then

sn⁡(u,m)=1y0cn⁡(u,m)=1−1y02dn⁡(u,m)=1−my02

as N->∞.

## Approximation in terms of hyperbolic functions

The Jacobi elliptic functions can be expanded in terms of the hyperbolic functions. When m is close to unity, such that m'2 and higher powers of m′ can be neglected, we have:[27][28]

* sn(u): sn⁡(u,m)≈tanh⁡(u)+14m′(sinh⁡(u)cosh⁡(u)−u)sech2(u).
* cn(u): cn⁡(u,m)≈sech⁡(u)−14m′(sinh⁡(u)cosh⁡(u)−u)tanh⁡(u)sech⁡(u).
* dn(u): dn⁡(u,m)≈sech⁡(u)+14m′(sinh⁡(u)cosh⁡(u)+u)tanh⁡(u)sech⁡(u).

For the Jacobi amplitude, am⁡(u,m)≈gd⁡(u)+14m′(sinh⁡(u)cosh⁡(u)−u)sech⁡(u).

## Continued fractions

Assuming real numbers a,p with 0<a<p and the nome q=eπiτ, Im⁡(τ)>0 with elliptic modulus k(τ)=1−k′(τ)2=(ϑ10(0;τ)/ϑ00(0;τ))2. If K[τ]=K(k(τ)), where K(x)=π/2⋅2F1(1/2,1/2;1;x2) is the complete elliptic integral of the first kind, then holds the following continued fraction expansion[29]

dn((p/2−a)τK[pτ2];k(pτ2))k′(pτ2)=∑n=−∞∞qp/2n2+(p/2−a)n∑n=−∞∞(−1)nqp/2n2+(p/2−a)n=−1+21−qa+qp−a1−qp+(qa+q2p−a)(qa+p+qp−a)1−q3p+qp(qa+q3p−a)(qa+2p+qp−a)1−q5p+q2p(qa+q4p−a)(qa+3p+qp−a)1−q7p+⋯

Known continued fractions involving sn(t),cn(t) and dn(t) with elliptic modulus k are

For z∈ℂ, |k|<1:[30] pg. 374

∫0∞sn(t)e−tzdt=112(1+k2)+z2−1⋅22⋅3k232(1+k2)+z2−3⋅42⋅5k252(1+k2)+z2−⋯

For z∈ℂ∖{0}, |k|<1:[30] pg. 375

∫0∞sn2(t)e−tzdt=2z−122(1+k2)+z2−2⋅32⋅4k242(1+k2)+z2−4⋅52⋅6k262(1+k2)+z2−⋯

For z∈ℂ∖{0}, |k|<1:[31] pg. 220

∫0∞cn(t)e−tzdt=1z+12z+22k2z+32z+42k2z+52z+⋯

For z∈ℂ∖{0}, |k|<1:[30] pg. 374

∫0∞dn(t)e−tzdt=1z+12k2z+22z+32k2z+42z+52k2z+⋯

For z∈ℂ, |k|<1:[30] pg. 375

∫0∞sn(t)cn(t)dn(t)e−tzdt=12⋅12(2−k2)+z2−1⋅22⋅3k42⋅32(2−k2)+z2−3⋅42⋅5k42⋅52(2−k2)+z2−⋯

## Inverse functions

The inverses of the Jacobi elliptic functions can be defined similarly to the inverse trigonometric functions; if x=sn⁡(ξ,m), ξ=arcsn⁡(x,m). They can be represented as elliptic integrals,[32][33][34] and power series representations have been found.[35][3]

* arcsn⁡(x,m)=∫0xdt(1−t2)(1−mt2)
* arccn⁡(x,m)=∫x1dt(1−t2)(1−m+mt2)
* arcdn⁡(x,m)=∫x1dt(1−t2)(t2+m−1)

## Map projection

The Peirce quincuncial projection is a map projection based on Jacobian elliptic functions.

## See also

* Elliptic curve
* Schwarz–Christoffel mapping
* Carlson symmetric form
* Jacobi theta function
* Ramanujan theta function
* Dixon elliptic functions
* Abel elliptic functions
* Weierstrass elliptic functions
* Lemniscate elliptic functions

## Notes

1. ↑ If u∈ℝ and m is restricted to [0,1], then dn⁡(u,m) can be also written as 1−msin2am⁡(u,m).

## References

* Abramowitz, Milton; Stegun, Irene Ann, eds (1983). "Chapter 16". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. pp. 569. LCCN 65-12253. ISBN 978-0-486-61272-0. http://www.math.sfu.ca/~cbm/aands/page_569.htm.
* Jonathan M. Borwein, Peter B. Borwein: Pi and the AGM. A study in analytical Number Theory and Computational Complexity. John Wiley & Sons, 1987.
* N. I. Akhiezer, Elements of the Theory of Elliptic Functions (1970) Moscow, translated into English as AMS Translations of Mathematical Monographs Volume 79 (1990) AMS, Rhode Island ISBN 0-8218-4532-2
* A. C. Dixon The elementary properties of the elliptic functions, with examples (Macmillan, 1894)
* Alfred George Greenhill The applications of elliptic functions (London, New York, Macmillan, 1892)
* Edmund T. Whittaker, George Neville Watson: A Course in Modern Analysis. 4th ed. Cambridge, England: Cambridge University Press, 1990. pp. 469–470.
* H. Hancock Lectures on the theory of elliptic functions (New York, J. Wiley & sons, 1910)
* Jacobi, C. G. J. (1829) (in la), Fundamenta nova theoriae functionum ellipticarum, Königsberg, Reprinted by Cambridge University Press 2012, ISBN 978-1-108-05200-9, https://archive.org/details/fundamentanovat00jacogoog
* Reinhardt, William P.; Walker, Peter L. (2010), "Jacobian Elliptic Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/22
* (in French) P. Appell and E. Lacour Principes de la théorie des fonctions elliptiques et applications (Paris, Gauthier Villars, 1897)
* (in French) G. H. Halphen Traité des fonctions elliptiques et de leurs applications (vol. 1) (Paris, Gauthier-Villars, 1886–1891)
* (in French) G. H. Halphen Traité des fonctions elliptiques et de leurs applications (vol. 2) (Paris, Gauthier-Villars, 1886–1891)
* (in French) G. H. Halphen Traité des fonctions elliptiques et de leurs applications (vol. 3) (Paris, Gauthier-Villars, 1886–1891)
* (in French) J. Tannery and J. Molk  Eléments de la théorie des fonctions elliptiques. Tome I, Introduction. Calcul différentiel. Ire partie (Paris : Gauthier-Villars et fils, 1893)
* (in French) J. Tannery and J. Molk  Eléments de la théorie des fonctions elliptiques. Tome II, Calcul différentiel. IIe partie (Paris : Gauthier-Villars et fils, 1893)
* (in French) J. Tannery and J. Molk  Eléments de la théorie des fonctions elliptiques. Tome III, Calcul intégral. Ire partie, Théorèmes généraux. Inversion (Paris : Gauthier-Villars et fils, 1893)
* (in French) J. Tannery and J. Molk  Eléments de la théorie des fonctions elliptiques. Tome IV, Calcul intégral. IIe partie, Applications (Paris : Gauthier-Villars et fils, 1893)
* (in French) C. Briot and J. C. Bouquet Théorie des fonctions elliptiques (Paris : Gauthier-Villars, 1875)
* Toshio Fukushima: Fast Computation of Complete Elliptic Integrals and Jacobian Elliptic Functions. 2012, National Astronomical Observatory of Japan (国立天文台)
* Lowan, Blanch und Horenstein: On the Inversion of the q-Series Associated with Jacobian Elliptic Functions. Bull. Amer. Math. Soc. 48, 1942
* H. Ferguson, D. E. Nielsen, G. Cook: A partition formula for the integer coefficients of the theta function nome. Mathematics of computation, Volume 29, Nummer 131, July 1975
* J. D. Fenton and R. S. Gardiner-Garden: Rapidly-convergent methods for evaluating elliptic integrals and theta and elliptic functions. J. Austral. Math. Soc. (Series B) 24, 1982, p. 57
* Adolf Kneser: Neue Untersuchung einer Reihe aus der Theorie der elliptischen Funktionen. J. reine u. angew. Math. 157, 1927. pages 209 – 218

## External links

* Hazewinkel, Michiel, ed. (2001), "Jacobi elliptic functions", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4, https://www.encyclopediaofmath.org/index.php?title=p/j054050
* Weisstein, Eric W.. "Jacobi Elliptic Functions". http://mathworld.wolfram.com/JacobiEllipticFunctions.html.

1. ↑ Armitage, J. V.; Eberlein, W. F. (2006). Elliptic Functions (First ed.). Cambridge University Press. ISBN 978-0-521-78078-0. p. 48
2. ↑ 2.0 2.1 Walker, Peter (2003). "The Analyticity of Jacobian Functions with Respect to the Parameter k". Proceedings of the Royal Society 459 (2038): 2569–2574. doi:10.1098/rspa.2003.1157. Bibcode: 2003RSPSA.459.2569W. https://www.jstor.org/stable/3560143.
3. ↑ 3.0 3.1 3.2 3.3 3.4 Olver, F. W. J., ed (2017-12-22). "NIST Digital Library of Mathematical Functions (Release 1.0.17)". National Institute of Standards and Technology. http://dlmf.nist.gov/22.
4. ↑ "cplot, Python package for plotting complex-valued functions". https://github.com/nschloe/cplot.
5. ↑ 5.0 5.1 5.2 5.3 5.4 5.5 Neville, Eric Harold (1944). Jacobian Elliptic Functions. Oxford: Oxford University Press. https://archive.org/details/jacobianelliptic00neviuoft.
6. ↑ 6.0 6.1 Sala, Kenneth L. (November 1989). "Transformations of the Jacobian Amplitude Function and Its Calculation via the Arithmetic-Geometric Mean". SIAM Journal on Mathematical Analysis 20 (6): 1514–1528. doi:10.1137/0520100. https://epubs.siam.org/doi/abs/10.1137/0520100.
7. ↑ Reinhardt, W. P.; Walker, P. L. (2010), "Jacobian Elliptic Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/22.16.E17
8. ↑ Carlson, B. C. (2010), "Elliptic Integrals", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/19.8.E13
9. ↑ Whittaker and Watson (1990) pp. 469–470
10. ↑ Alvaro H. Salas, Lorenzo J. H. Martinez, David L. R. Ocampo R. (2021-10-11), "Approximation of Elliptic Functions by Means of Trigonometric Functions with Applications", Mathematical Problems in Engineering 2021: e5546666, doi:10.1155/2021/5546666, ISSN 1024-123X
11. ↑ "Table of Infinite Products Infinite Sums Infinite Series Elliptic Theta". http://www-elsa.physik.uni-bonn.de/~dieckman/InfProd/InfProd.html#specialxvaluesxofxEllipticKxandxEllipticE.
12. ↑ "A002103 - OEIS". https://oeis.org/A002103.
13. ↑ "Series Expansion of EllipticNomeQ differs from older Mathematica Version" (in en). https://mathematica.stackexchange.com/questions/269455/series-expansion-of-ellipticnomeq-differs-from-older-mathematica-version.
14. ↑ R. B. King, E. R. Canfield (1992-08-01), "Icosahedral symmetry and the quintic equation", Computers & Mathematics with Applications 24 (3): 13–28, doi:10.1016/0898-1221(92)90210-9, ISSN 0898-1221
15. ↑ Karl Heinrich Schellbach (1864), Die Lehre von den Elliptischen Integralen und den ThetaFunctionen, G. Reimer, http://archive.org/details/bub_gb_WiYPAAAAQAAJ, retrieved 2023-06-11
16. ↑ Adolf Kneser (1927), "Neue Untersuchung einer Reihe aus der Theorie der elliptischen Funktionen.", Journal für die reine und angewandte Mathematik 158: 209–218, ISSN 0075-4102, https://eudml.org/doc/149645, retrieved 2023-06-11
17. ↑ D. K. Lee (1989-03-01), Application of theta functions for numerical evaluation of complete elliptic integrals of the first and second kinds, Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), doi:10.2172/6137964, https://www.osti.gov/biblio/6137964, retrieved 2023-06-11
18. ↑ 18.0 18.1 18.2 18.3 18.4 18.5 "Introduction to the Jacobi elliptic functions". Wolfram Research, Inc.. 2018. http://functions.wolfram.com/EllipticFunctions/JacobiAmplitude/introductions/JacobiPQs/ShowAll.html.
19. ↑ Whittaker, E.T.; Watson, G.N. (1940). A Course in Modern Analysis. New York, USA: The MacMillan Co.. ISBN 978-0-521-58807-2. https://archive.org/details/courseofmodernan00whit.
20. ↑ "Elliptic Functions: Complex Variables". https://paramanands.blogspot.co.uk/2011/01/elliptic-functions-complex-variables.html#.WlHhTbp2t9A.
21. ↑ 21.0 21.1 Viktor Prasolov and Yuri Soloviev (1997), "Elliptic Functions and Elliptic Integrals", arXiv: General Mathematics, https://staff.math.su.se/mleites/books/prasolov-soloviev-1997-elliptic.pdf, retrieved 2023-06-24
22. ↑ Weisstein, Eric W.. "Jacobi Theta Functions". http://mathworld.wolfram.com/JacobiThetaFunctions.html.
23. ↑ Borwein, Jonathan M.; Borwein, Peter B. (1998). Pi and the AGM: a study in analytic number theory and computational complexity. New York: Wiley. ISBN 978-0471315155. http://wayback.cecm.sfu.ca/~pborwein/TEMP_PROTECTED/pi-agm.pdf.
24. ↑ "DLMF: 20.5 Infinite Products and Related Results". https://dlmf.nist.gov/20.5.
25. ↑ Schett, Alois (1976). "Properties of the Taylor series expansion coefficients of the Jacobian Elliptic Functions". Math. Comp. 30 (133): 143–147. doi:10.1090/S0025-5718-1976-0391477-3.
26. ↑ Salzer, Herbert E. (July 1962). "Quick calculation of Jacobian elliptic functions". Communications of the ACM 5 (7): 399. doi:10.1145/368273.368573.
27. ↑ Reinhardt, W. P.; Walker, P. L. (2010), "Jacobian Elliptic Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/22.10
28. ↑ Reinhardt, W. P.; Walker, P. L. (2010), "Jacobian Elliptic Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/22.16.E8
29. ↑ N.Bagis.(2020)."Evaluations of series related to Jacobi elliptic functions". preprint https://www.researchgate.net/publication/331370071_Evaluations_of_Series_Related_to_Jacobi_Elliptic_Functions
30. ↑ 30.0 30.1 30.2 30.3 H.S. Wall. (1948). "Analytic Theory of Continued Fractions", Van Nostrand, New York.
31. ↑ Perron, O. (1957). "Die Lehre von den Kettenbruchen", Band II, B.G. Teubner, Stuttgart.
32. ↑ Reinhardt, W. P.; Walker, P. L. (2010), "§22.15 Inverse Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, http://dlmf.nist.gov/22.15
33. ↑ Ehrhardt, Wolfgang. "The AMath and DAMath Special Functions: Reference Manual and Implementation Notes". p. 42. http://www.wolfgang-ehrhardt.de/specialfunctions.pdf.
34. ↑ Byrd, P.F.; Friedman, M.D. (1971). Handbook of Elliptic Integrals for Engineers and Scientists (2nd ed.). Berlin: Springer-Verlag.
35. ↑ Carlson, B. C. (2008). "Power series for inverse Jacobian elliptic functions". Mathematics of Computation 77 (263): 1615–1621. doi:10.1090/s0025-5718-07-02049-2. Bibcode: 2008MaCom..77.1615C. https://www.ams.org/journals/mcom/2008-77-263/S0025-5718-07-02049-2/S0025-5718-07-02049-2.pdf. Retrieved 17 July 2013.

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Jacobi elliptic functions. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Jacobi_elliptic_functions&oldid=337090"