Lemniscate elliptic functions

From testwiki
Jump to navigation Jump to search

Template:Short description

The lemniscate sine (red) and lemniscate cosine (purple) applied to a real argument, in comparison with the trigonometric sine Template:Math (pale dashed red).

In mathematics, the lemniscate elliptic functions are elliptic functions related to the arc length of the lemniscate of Bernoulli. They were first studied by Giulio Fagnano in 1718 and later by Leonhard Euler and Carl Friedrich Gauss, among others.[1]

The lemniscate sine and lemniscate cosine functions, usually written with the symbols Template:Math and Template:Math (sometimes the symbols Template:Math and Template:Math or Template:Math and Template:Math are used instead),[2] are analogous to the trigonometric functions sine and cosine. While the trigonometric sine relates the arc length to the chord length in a unit-diameter circle x2+y2=x,[3] the lemniscate sine relates the arc length to the chord length of a lemniscate (x2+y2)2=x2−y2.

The lemniscate functions have periods related to a number Template:Math called the lemniscate constant, the ratio of a lemniscate's perimeter to its diameter. This number is a quartic analog of the (quadratic) Template:Math, ratio of perimeter to diameter of a circle.

As complex functions, Template:Math and Template:Math have a square period lattice (a multiple of the Gaussian integers) with fundamental periods {(1+i)ϖ,(1−i)ϖ},[4] and are a special case of two Jacobi elliptic functions on that lattice, sl⁡z=sn⁡(z;−1), cl⁡z=cd⁡(z;−1).

Similarly, the hyperbolic lemniscate sine Template:Math and hyperbolic lemniscate cosine Template:Math have a square period lattice with fundamental periods {2ϖ,2ϖi}.

The lemniscate functions and the hyperbolic lemniscate functions are related to the Weierstrass elliptic function ℘(z;a,0).

Lemniscate sine and cosine functions

Definitions

The lemniscate functions Template:Math and Template:Math can be defined as the solution to the initial value problem:[5]

ddzsl⁡z=(1+sl2z)cl⁡z, ddzcl⁡z=−(1+cl2z)sl⁡z, sl⁡0=0, cl⁡0=1,

or equivalently as the inverses of an elliptic integral, the Schwarz–Christoffel map from the complex unit disk to a square with corners {12ϖ,12ϖi,−12ϖ,−12ϖi}:[6]

z=∫0sl⁡zdt1−t4=∫cl⁡z1dt1−t4.

Beyond that square, the functions can be extended to the complex plane via analytic continuation by successive reflections.

By comparison, the circular sine and cosine can be defined as the solution to the initial value problem:

ddzsin⁡z=cos⁡z, ddzcos⁡z=−sin⁡z, sin⁡0=0, cos⁡0=1,

or as inverses of a map from the upper half-plane to a half-infinite strip with real part between −12π,12π and positive imaginary part:

z=∫0sin⁡zdt1−t2=∫cos⁡z1dt1−t2.

Relation to the lemniscate constant

Template:Main

The lemniscate sine function and hyperbolic lemniscate sine functions are defined as inverses of elliptic integrals. The complete integrals are related to the lemniscate constant Template:Mvar.

The lemniscate functions have minimal real period Template:Tmath, minimal imaginary period Template:Tmath and fundamental complex periods (1+i)ϖ and (1−i)ϖ for a constant Template:Tmath called the lemniscate constant,[7]

ϖ=2∫01dt1−t4=2π322(Γ(34))2=2.62205…

The lemniscate functions satisfy the basic relation cl⁡z=sl(12ϖ−z), analogous to the relation cos⁡z=sin(12π−z).

The lemniscate constant Template:Tmath is a close analog of the [[pi|circle constant Template:Tmath]], and many identities involving Template:Tmath have analogues involving Template:Tmath, as identities involving the trigonometric functions have analogues involving the lemniscate functions. For example, Viète's formula for Template:Tmath can be written:[8]

2π=12⋅12+1212⋅12+1212+1212⋯

An analogous formula for Template:Tmath is:[9]

2ϖ=12⋅12+12/12⋅12+12/12+12/12⋯

The Machin formula for Template:Tmath is 14π=4arctan⁡15−arctan⁡1239, and several similar formulas for Template:Tmath can be developed using trigonometric angle sum identities, e.g. Euler's formula 14π=arctan⁡12+arctan⁡13. Analogous formulas can be developed for Template:Tmath, including the following found by Gauss: 12ϖ=2arcsl⁡12+arcsl⁡723.[10]

The lemniscate and circle constants were found by Gauss to be related to each-other by the arithmetic-geometric mean Template:Tmath:[11]

πϖ=M(1,2)

Argument identities

Zeros, poles and symmetries

sl in the complex plane.[12] In the picture, it can be seen that the fundamental periods (1+i)ϖ and (1−i)ϖ are "minimal" in the sense that they have the smallest absolute value of all periods whose real part is non-negative.

The lemniscate functions Template:Math and Template:Math are even and odd functions, respectively,

cl⁡(−z)=cl⁡zsl⁡(−z)=−sl⁡z

At translations of 12ϖ, Template:Math and Template:Math are exchanged, and at translations of 12iϖ they are additionally rotated and reciprocated:[13]

cl(z±12ϖ)=∓sl⁡z,cl(z±12iϖ)=∓isl⁡zsl(z±12ϖ)=±cl⁡z,sl(z±12iϖ)=±icl⁡z

Doubling these to translations by a unit-Gaussian-integer multiple of ϖ (that is, ±ϖ or ±iϖ), negates each function, an involution:

cl⁡(z+ϖ)=cl⁡(z+iϖ)=−cl⁡zsl⁡(z+ϖ)=sl⁡(z+iϖ)=−sl⁡z

As a result, both functions are invariant under translation by an even-Gaussian-integer multiple of ϖ.[14] That is, a displacement (a+bi)ϖ, with a+b=2k for integers Template:Tmath, Template:Tmath, and Template:Tmath.

cl(z+(1+i)ϖ)=cl(z+(1−i)ϖ)=cl⁡zsl(z+(1+i)ϖ)=sl(z+(1−i)ϖ)=sl⁡z

This makes them elliptic functions (doubly periodic meromorphic functions in the complex plane) with a diagonal square period lattice of fundamental periods (1+i)ϖ and (1−i)ϖ.[15] Elliptic functions with a square period lattice are more symmetrical than arbitrary elliptic functions, following the symmetries of the square.

Reflections and quarter-turn rotations of lemniscate function arguments have simple expressions:

cl⁡z¯=cl⁡z‾sl⁡z¯=sl⁡z‾cl⁡iz=1cl⁡zsl⁡iz=isl⁡z

The Template:Math function has simple zeros at Gaussian integer multiples of Template:Tmath, complex numbers of the form aϖ+bϖi for integers Template:Tmath and Template:Tmath. It has simple poles at Gaussian half-integer multiples of Template:Tmath, complex numbers of the form (a+12)ϖ+(b+12)ϖi, with residues (−1)a−b+1i. The Template:Math function is reflected and offset from the Template:Math function, cl⁡z=sl(12ϖ−z). It has zeros for arguments (a+12)ϖ+bϖi and poles for arguments aϖ+(b+12)ϖi, with residues (−1)a−bi.

Also

sl⁡z=sl⁡w↔z=(−1)m+nw+(m+ni)ϖ

for some m,n∈ℤ and

sl⁡((1±i)z)=(1±i)sl⁡zsl′⁡z.

The last formula is a special case of complex multiplication. Analogous formulas can be given for sl⁡((n+mi)z) where n+mi is any Gaussian integer – the function sl has complex multiplication by ℤ[i].[16]

There are also infinite series reflecting the distribution of the zeros and poles of Template:Math:[17][18]

1sl⁡z=∑(n,k)∈ℤ2(−1)n+kz+nϖ+kϖi
sl⁡z=−i∑(n,k)∈ℤ2(−1)n+kz+(n+1/2)ϖ+(k+1/2)ϖi.

Pythagorean-like identity

Curves Template:Math for various values of a. Negative Template:Math in green, positive Template:Math in blue, Template:Math in red, Template:Math in black.

The lemniscate functions satisfy a Pythagorean-like identity:

cl2⁡z+sl2⁡z+cl2⁡zsl2⁡z=1

As a result, the parametric equation (x,y)=(cl⁡t,sl⁡t) parametrizes the quartic curve x2+y2+x2y2=1.

This identity can alternately be rewritten:[19]

(1+cl2⁡z)(1+sl2⁡z)=2
cl2⁡z=1−sl2⁡z1+sl2⁡z,sl2⁡z=1−cl2⁡z1+cl2⁡z

Defining a tangent-sum operator as a⊕b:=tan⁡(arctan⁡a+arctan⁡b)=a+b1−ab, gives:

cl2⁡z⊕sl2⁡z=1.

Derivatives and integrals

The derivatives are as follows:

ddzcl⁡z=cl'⁡z=−(1+cl2⁡z)sl⁡z=−2sl⁡zsl2z+1cl'2⁡z=1−cl4⁡zddzsl⁡z=sl'⁡z=(1+sl2⁡z)cl⁡z=2cl⁡zcl2z+1sl'2⁡z=1−sl4⁡z
ddzcl~z=−2sl~zcl⁡z−sl~zcl⁡zddzsl~z=2cl~zcl⁡z−cl~zcl⁡z

The second derivatives of lemniscate sine and lemniscate cosine are their negative duplicated cubes:

d2dz2cl⁡z=−2cl3⁡z
d2dz2sl⁡z=−2sl3⁡z

The lemniscate functions can be integrated using the inverse tangent function:

∫cl⁡zdz=arctan⁡sl⁡z+C∫sl⁡zdz=−arctan⁡cl⁡z+C∫cl~zdz=sl~zcl⁡z+C∫sl~zdz=−cl~zcl⁡z+C

Argument sum and multiple identities

Like the trigonometric functions, the lemniscate functions satisfy argument sum and difference identities. The original identity used by Fagnano for bisection of the lemniscate was:[20]

sl⁡(u+v)=sl⁡usl'⁡v+sl⁡vsl'⁡u1+sl2⁡usl2⁡v

The derivative and Pythagorean-like identities can be used to rework the identity used by Fagano in terms of Template:Math and Template:Math. Defining a tangent-sum operator a⊕b:=tan⁡(arctan⁡a+arctan⁡b) and tangent-difference operator a⊖b:=a⊕(−b), the argument sum and difference identities can be expressed as:[21]

cl⁡(u+v)=cl⁡ucl⁡v⊖sl⁡usl⁡v=cl⁡ucl⁡v−sl⁡usl⁡v1+sl⁡ucl⁡usl⁡vcl⁡vcl⁡(u−v)=cl⁡ucl⁡v⊕sl⁡usl⁡vsl⁡(u+v)=sl⁡ucl⁡v⊕cl⁡usl⁡v=sl⁡ucl⁡v+cl⁡usl⁡v1−sl⁡ucl⁡usl⁡vcl⁡vsl⁡(u−v)=sl⁡ucl⁡v⊖cl⁡usl⁡v

These resemble their trigonometric analogs:

cos⁡(u±v)=cos⁡ucos⁡v∓sin⁡usin⁡vsin⁡(u±v)=sin⁡ucos⁡v±cos⁡usin⁡v

In particular, to compute the complex-valued functions in real components,

cl⁡(x+iy)=cl⁡x−isl⁡xsl⁡ycl⁡ycl⁡y+isl⁡xcl⁡xsl⁡y=cl⁡xcl⁡y(1−sl2xsl2y)cl2y+sl2xcl2xsl2y−isl⁡xsl⁡y(cl2x+cl2y)cl2y+sl2xcl2xsl2ysl⁡(x+iy)=sl⁡x+icl⁡xsl⁡ycl⁡ycl⁡y−isl⁡xcl⁡xsl⁡y=sl⁡xcl⁡y(1−cl2xsl2y)cl2y+sl2xcl2xsl2y+icl⁡xsl⁡y(sl2x+cl2y)cl2y+sl2xcl2xsl2y

Gauss discovered that

sl⁡(u−v)sl⁡(u+v)=sl⁡((1+i)u)−sl⁡((1+i)v)sl⁡((1+i)u)+sl⁡((1+i)v)

where u,v∈ℂ such that both sides are well-defined.

Also

sl⁡(u+v)sl⁡(u−v)=sl2u−sl2v1+sl2usl2v

where u,v∈ℂ such that both sides are well-defined; this resembles the trigonometric analog

sin⁡(u+v)sin⁡(u−v)=sin2u−sin2v.

Bisection formulas:

cl212x=1+cl⁡x1+sl2x1+1+sl2x
sl212x=1−cl⁡x1+sl2x1+1+sl2x

Duplication formulas:[22]

cl⁡2x=−1+2cl2x+cl4x1+2cl2x−cl4x
sl⁡2x=2sl⁡xcl⁡x1+sl2x1+sl4x

Triplication formulas:[22]

cl⁡3x=−3cl⁡x+6cl5x+cl9x1+6cl4x−3cl8x
sl⁡3x=3sl⁡x−6sl5x−1sl9x1+6sl4x−3sl8x

Note the "reverse symmetry" of the coefficients of numerator and denominator of sl⁡3x. This phenomenon can be observed in multiplication formulas for sl⁡βx where β=m+ni whenever m,n∈ℤ and m+n is odd.[16]

Lemnatomic polynomials

Let L be the lattice

L=ℤ(1+i)ϖ+ℤ(1−i)ϖ.

Furthermore, let K=ℚ(i), 𝒪=ℤ[i], z∈ℂ, β=m+in, γ=m′+in′ (where m,n,m′,n′∈ℤ), m+n be odd, m′+n′ be odd, γ≡1mod⁡2(1+i) and sl⁡βz=Mβ(sl⁡z). Then

Mβ(x)=iεxPβ(x4)Qβ(x4)

for some coprime polynomials Pβ(x),Qβ(x)∈𝒪[x] and some ε∈{0,1,2,3}[23] where

xPβ(x4)=∏γ|βΛγ(x)

and

Λβ(x)=∏[α]∈(𝒪/β𝒪)×(x−sl⁡αδβ)

where δβ is any β-torsion generator (i.e. δβ∈(1/β)L and [δβ]∈(1/β)L/L generates (1/β)L/L as an 𝒪-module). Examples of β-torsion generators include 2ϖ/β and (1+i)ϖ/β. The polynomial Λβ(x)∈𝒪[x] is called the β-th lemnatomic polynomial. It is monic and is irreducible over K. The lemnatomic polynomials are the "lemniscate analogs" of the cyclotomic polynomials,[24]

Φk(x)=∏[a]∈(ℤ/kℤ)×(x−ζka).

The β-th lemnatomic polynomial Λβ(x) is the minimal polynomial of sl⁡δβ in K[x]. For convenience, let ωβ=sl⁡(2ϖ/β) and ω~β=sl⁡((1+i)ϖ/β). So for example, the minimal polynomial of ω5 (and also of ω~5) in K[x] is

Λ5(x)=x16+52x12−26x8−12x4+1,

and[25]

ω5=−13+65+285−3854
ω~5=−13−65+285+3854[26]

(an equivalent expression is given in the table below). Another example is[24]

Λ−1+2i(x)=x4−1+2i

which is the minimal polynomial of ω−1+2i (and also of ω~−1+2i) in K[x].

If p is prime and β is positive and odd,[27] then[28]

deg⁡Λβ=β2∏p|β(1−1p)(1−(−1)(p−1)/2p)

which can be compared to the cyclotomic analog

deg⁡Φk=k∏p|k(1−1p).

Specific values

Just as for the trigonometric functions, values of the lemniscate functions can be computed for divisions of the lemniscate into Template:Tmath parts of equal length, using only basic arithmetic and square roots, if and only if Template:Tmath is of the form n=2kp1p2⋯pm where Template:Tmath is a non-negative integer and each Template:Tmath (if any) is a distinct Fermat prime.[29]

n cl⁡nϖ sl⁡nϖ
1 −1 0
56 −23−34 12(3+1−124)
34 −2−1 2−1
23 −12(3+1−124) 23−34
12 0 1
13 12(3+1−124) 23−34
14 2−1 2−1
16 23−34 12(3+1−124)

Relation to geometric shapes

Arc length of Bernoulli's lemniscate

The lemniscate sine and cosine relate the arc length of an arc of the lemniscate to the distance of one endpoint from the origin.
The trigonometric sine and cosine analogously relate the arc length of an arc of a unit-diameter circle to the distance of one endpoint from the origin.

ℒ, the lemniscate of Bernoulli with unit distance from its center to its furthest point (i.e. with unit "half-width"), is essential in the theory of the lemniscate elliptic functions. It can be characterized in at least three ways:

Angular characterization: Given two points A and B which are unit distance apart, let B′ be the reflection of B about A. Then ℒ is the closure of the locus of the points P such that |APB−APB′| is a right angle.[30]

Focal characterization: ℒ is the locus of points in the plane such that the product of their distances from the two focal points F1=(−12,0) and F2=(12,0) is the constant 12.

Explicit coordinate characterization: ℒ is a quartic curve satisfying the polar equation r2=cos⁡2θ or the Cartesian equation (x2+y2)2=x2−y2.

The perimeter of ℒ is 2ϖ.[31]

The points on ℒ at distance r from the origin are the intersections of the circle x2+y2=r2 and the hyperbola x2−y2=r4. The intersection in the positive quadrant has Cartesian coordinates:

(x(r),y(r))=(12r2(1+r2),12r2(1−r2)).

Using this parametrization with r∈[0,1] for a quarter of ℒ, the arc length from the origin to a point (x(r),y(r)) is:[32]

∫0rx′(t)2+y′(t)2dt=∫0r(1+2t2)22(1+t2)+(1−2t2)22(1−t2)dt=∫0rdt1−t4=arcsl⁡r.

Likewise, the arc length from (1,0) to (x(r),y(r)) is:

∫r1x′(t)2+y′(t)2dt=∫r1dt1−t4=arccl⁡r=12ϖ−arcsl⁡r.

Or in the inverse direction, the lemniscate sine and cosine functions give the distance from the origin as functions of arc length from the origin and the point (1,0), respectively.

Analogously, the circular sine and cosine functions relate the chord length to the arc length for the unit diameter circle with polar equation r=cos⁡θ or Cartesian equation x2+y2=x, using the same argument above but with the parametrization:

(x(r),y(r))=(r2,r2(1−r2)).

Alternatively, just as the unit circle x2+y2=1 is parametrized in terms of the arc length s from the point (1,0) by

(x(s),y(s))=(cos⁡s,sin⁡s),

ℒ is parametrized in terms of the arc length s from the point (1,0) by[33]

(x(s),y(s))=(cl⁡s1+sl2s,sl⁡scl⁡s1+sl2s)=(cl~s,sl~s).

The notation cl~,sl~ is used solely for the purposes of this article; in references, notation for general Jacobi elliptic functions is used instead.

The lemniscate integral and lemniscate functions satisfy an argument duplication identity discovered by Fagnano in 1718:[34]

∫0zdt1−t4=2∫0udt1−t4,if z=2u1−u41+u4 and 0≤u≤2−1.
A lemniscate divided into 15 sections of equal arclength (red curves). Because the prime factors of 15 (3 and 5) are both Fermat primes, this polygon (in black) is constructible using a straightedge and compass.

Later mathematicians generalized this result. Analogously to the constructible polygons in the circle, the lemniscate can be divided into Template:Tmath sections of equal arc length using only straightedge and compass if and only if Template:Tmath is of the form n=2kp1p2⋯pm where Template:Tmath is a non-negative integer and each Template:Tmath (if any) is a distinct Fermat prime.[35] The "if" part of the theorem was proved by Niels Abel in 1827–1828, and the "only if" part was proved by Michael Rosen in 1981.[36] Equivalently, the lemniscate can be divided into Template:Tmath sections of equal arc length using only straightedge and compass if and only if φ(n) is a power of two (where φ is Euler's totient function). The lemniscate is not assumed to be already drawn, as that would go against the rules of straightedge and compass constructions; instead, it is assumed that we are given only two points by which the lemniscate is defined, such as its center and radial point (one of the two points on the lemniscate such that their distance from the center is maximal) or its two foci.

Let rj=sl⁡2jϖn. Then the Template:Tmath-division points for ℒ are the points

(rj12(1+rj2), (−1)⌊4j/n⌋12rj2(1−rj2)),j∈{1,2,…,n}

where ⌊⋅⌋ is the floor function. See below for some specific values of sl⁡2ϖn.

Arc length of rectangular elastica

The lemniscate sine relates the arc length to the x coordinate in the rectangular elastica.

The inverse lemniscate sine also describes the arc length Template:Tmath relative to the Template:Tmath coordinate of the rectangular elastica.[37] This curve has Template:Tmath coordinate and arc length:

y=∫x1t2dt1−t4,s=arcsl⁡x=∫0xdt1−t4

The rectangular elastica solves a problem posed by Jacob Bernoulli, in 1691, to describe the shape of an idealized flexible rod fixed in a vertical orientation at the bottom end, and pulled down by a weight from the far end until it has been bent horizontal. Bernoulli's proposed solution established Euler–Bernoulli beam theory, further developed by Euler in the 18th century.

Elliptic characterization

The lemniscate elliptic functions and an ellipse

Let C be a point on the ellipse x2+2y2=1 in the first quadrant and let D be the projection of C on the unit circle x2+y2=1. The distance r between the origin A and the point C is a function of φ (the angle BAC where B=(1,0); equivalently the length of the circular arc BD). The parameter u is given by

u=∫0φr(θ)dθ=∫0φdθ1+sin2θ.

If E is the projection of D on the x-axis and if F is the projection of C on the x-axis, then the lemniscate elliptic functions are given by

cl⁡u=AF‾,sl⁡u=DE‾,
cl~u=AF‾AC‾,sl~u=AF‾FC‾.

Series Identities

Power series

The power series expansion of the lemniscate sine at the origin is[38]

sl⁡z=∑n=0∞anzn=z−12z55!+3024z99!−4390848z1313!+⋯,|z|<ϖ2

where the coefficients an are determined as follows:

n≢1(mod4)⟹an=0,
a1=1,∀n∈ℕ0:an+2=−2(n+1)(n+2)∑i+j+k=naiajak

where i+j+k=n stands for all three-term compositions of n. For example, to evaluate a13, it can be seen that there are only six compositions of 13−2=11 that give a nonzero contribution to the sum: 11=9+1+1=1+9+1=1+1+9 and 11=5+5+1=5+1+5=1+5+5, so

a13=−212⋅13(a9a1a1+a1a9a1+a1a1a9+a5a5a1+a5a1a5+a1a5a5)=−1115600.

The expansion can be equivalently written as[39]

sl⁡z=∑n=0∞p2nz4n+1(4n+1)!,|z|<ϖ2

where

pn+2=−12∑j=0n(2n+22j+2)pn−j∑k=0j(2j+12k+1)pkpj−k,p0=1,p1=0.

The power series expansion of sl~ at the origin is

sl~z=∑n=0∞αnzn=z−9z33!+153z55!−4977z77!+⋯,|z|<ϖ2

where αn=0 if n is even and[40]

αn=2πϖ(−1)(n−1)/2n!∑k=1∞(2kπ/ϖ)n+1cosh⁡kπ,|αn|∼2n+5/2n+1ϖn+2

if n is odd.

The expansion can be equivalently written as[41]

sl~z=∑n=0∞(−1)n2n+1(∑l=0n2l(2n+22l+1)sltn−l)z2n+1(2n+1)!,|z|<ϖ2

where

sn+2=3sn+1+24∑j=0n(2n+22j+2)sn−j∑k=0j(2j+12k+1)sksj−k,s0=1,s1=3,
tn+2=3tn+1+3∑j=0n(2n+22j+2)tn−j∑k=0j(2j+12k+1)tktj−k,t0=1,t1=3.

For the lemniscate cosine,[42]

cl⁡z=1−∑n=0∞(−1)n(∑l=0n2l(2n+22l+1)qlrn−l)z2n+2(2n+2)!=1−2z22!+12z44!−216z66!+⋯,|z|<ϖ2,
cl~z=∑n=0∞(−1)n2nqnz2n(2n)!=1−3z22!+33z44!−819z66!+⋯,|z|<ϖ2

where

rn+2=3∑j=0n(2n+22j+2)rn−j∑k=0j(2j+12k+1)rkrj−k,r0=1,r1=0,
qn+2=32qn+1+6∑j=0n(2n+22j+2)qn−j∑k=0j(2j+12k+1)qkqj−k,q0=1,q1=32.

Ramanujan's cos/cosh identity

Ramanujan's famous cos/cosh identity states that if

R(s)=πϖ2∑n∈ℤcos⁡(2nπs/ϖ)cosh⁡nπ,

then[40]

R(s)−2+R(is)−2=2,|Re⁡s|<ϖ2,|Im⁡s|<ϖ2.

There is a close relation between the lemniscate functions and R(s). Indeed,[40][43]

sl~s=−ddsR(s)|Im⁡s|<ϖ2
cl~s=dds1−R(s)2,|Re⁡s−ϖ2|<ϖ2,|Im⁡s|<ϖ2

and

R(s)=11+sl2s,|Im⁡s|<ϖ2.

Continued fractions

For z∈ℂ∖{0}:[44]

∫0∞e−tz2cl⁡tdt=1/2z+a1z+a2z+a3z+⋱,an=n24((−1)n+1+3)
∫0∞e−tz2sl⁡tcl⁡tdt=1/2z2+b1−a1z2+b2−a2z2+b3−⋱,an=n2(4n2−1),bn=3(2n−1)2

Methods of computation

Template:Quote box

Several methods of computing sl⁡x involve first making the change of variables πx=ϖx~ and then computing sl⁡(ϖx~/π).

A hyperbolic series method:[45][46]

sl⁡(ϖπx)=πϖ∑n∈ℤ(−1)ncosh⁡(x−(n+1/2)π),x∈ℂ
1sl⁡(ϖx/π)=πϖ∑n∈ℤ(−1)nsinh(x−nπ)=πϖ∑n∈ℤ(−1)nsin⁡(x−nπi),x∈ℂ

Fourier series method:[47]

sl⁡(ϖπx)=2πϖ∑n=0∞(−1)nsin⁡((2n+1)x)cosh⁡((n+1/2)π),|Im⁡x|<π2
cl⁡(ϖπx)=2πϖ∑n=0∞cos⁡((2n+1)x)cosh⁡((n+1/2)π),|Im⁡x|<π2
1sl⁡(ϖx/π)=πϖ(1sin⁡x−4∑n=0∞sin⁡((2n+1)x)e(2n+1)π+1),|Im⁡x|<π

The lemniscate functions can be computed more rapidly by

sl⁡(ϖπx)=θ1(x,e−π)θ3(x,e−π),x∈ℂcl⁡(ϖπx)=θ2(x,e−π)θ4(x,e−π),x∈ℂ

where

θ1(x,e−π)=∑n∈ℤ(−1)n+1e−π(n+1/2+x/π)2=∑n∈ℤ(−1)ne−π(n+1/2)2sin⁡((2n+1)x),θ2(x,e−π)=∑n∈ℤ(−1)ne−π(n+x/π)2=∑n∈ℤe−π(n+1/2)2cos⁡((2n+1)x),θ3(x,e−π)=∑n∈ℤe−π(n+x/π)2=∑n∈ℤe−πn2cos⁡2nx,θ4(x,e−π)=∑n∈ℤe−π(n+1/2+x/π)2=∑n∈ℤ(−1)ne−πn2cos⁡2nx

are the Jacobi theta functions.[48]

Fourier series for the logarithm of the lemniscate sine:

ln⁡sl⁡(ϖπx)=ln⁡2−π4+ln⁡sin⁡x+2∑n=1∞(−1)ncos⁡2nxn(enπ+(−1)n),|Im⁡x|<π2

The following series identities were discovered by Ramanujan:[49]

ϖ2π2sl2(ϖx/π)=1sin2x−1π−8∑n=1∞ncos⁡2nxe2nπ−1,|Im⁡x|<π
arctan⁡sl⁡(ϖπx)=2∑n=0∞sin⁡((2n+1)x)(2n+1)cosh⁡((n+1/2)π),|Im⁡x|<π2

The functions sl~ and cl~ analogous to sin and cos on the unit circle have the following Fourier and hyperbolic series expansions:[40][43][50]

sl~s=22π2ϖ2∑n=1∞nsin⁡(2nπs/ϖ)cosh⁡nπ,|Im⁡s|<ϖ2
cl~s=2π2ϖ2∑n=0∞(2n+1)cos⁡((2n+1)πs/ϖ)sinh⁡((n+1/2)π),|Im⁡s|<ϖ2
sl~s=π2ϖ22∑n∈ℤsinh⁡(π(n+s/ϖ))cosh2(π(n+s/ϖ)),s∈ℂ
cl~s=π2ϖ22∑n∈ℤ(−1)ncosh2(π(n+s/ϖ)),s∈ℂ

The following identities come from product representations of the theta functions:[51]

sl(ϖπx)=2e−π/4sin⁡x∏n=1∞1−2e−2nπcos⁡2x+e−4nπ1+2e−(2n−1)πcos⁡2x+e−(4n−2)π,x∈ℂ
cl(ϖπx)=2e−π/4cos⁡x∏n=1∞1+2e−2nπcos⁡2x+e−4nπ1−2e−(2n−1)πcos⁡2x+e−(4n−2)π,x∈ℂ

A similar formula involving the sn function can be given.[52]

The lemniscate functions as a ratio of entire functions

Since the lemniscate sine is a meromorphic function in the whole complex plane, it can be written as a ratio of entire functions. Gauss showed that Template:Math has the following product expansion, reflecting the distribution of its zeros and poles:[53]

sl⁡z=M(z)N(z)

where

M(z)=z∏α(1−z4α4),N(z)=∏β(1−z4β4).

Here, α and β denote, respectively, the zeros and poles of Template:Math which are in the quadrant Re⁡z>0,Im⁡z≥0. A proof can be found in.[53][54] Importantly, the infinite products converge to the same value for all possible orders in which their terms can be multiplied, as a consequence of uniform convergence.[55] Template:Collapse top Proof by logarithmic differentiation

It can be easily seen (using uniform and absolute convergence arguments to justify interchanging of limiting operations) that

M′(z)M(z)=−∑n=0∞24nH4nz4n−1(4n)!,|z|<ϖ

(where Hn are the Hurwitz numbers defined in Lemniscate elliptic functions § Hurwitz numbers) and

N′(z)N(z)=(1+i)M′((1+i)z)M((1+i)z)−M′(z)M(z).

Therefore

N′(z)N(z)=∑n=0∞24n(1−(−1)n22n)H4nz4n−1(4n)!,|z|<ϖ2.

It is known that

1sl2z=∑n=0∞24n(4n−1)H4nz4n−2(4n)!,|z|<ϖ.

Then from

ddzsl′⁡zsl⁡z=−1sl2z−sl2z

and

sl2z=1sl2z−(1+i)2sl2((1+i)z)

we get

sl′⁡zsl⁡z=−∑n=0∞24n(2−(−1)n22n)H4nz4n−1(4n)!,|z|<ϖ2.

Hence

sl′⁡zsl⁡z=M′(z)M(z)−N′(z)N(z),|z|<ϖ2.

Therefore

sl⁡z=CM(z)N(z)

for some constant C for |z|<ϖ/2 but this result holds for all z∈ℂ by analytic continuation. Using

limz→0sl⁡zz=1

gives C=1 which completes the proof. ◼

Proof by Liouville's theorem

Let

f(z)=M(z)N(z)=(1+i)M(z)2M((1+i)z),

with patches at removable singularities. The shifting formulas

M(z+2ϖ)=e2πϖ(z+ϖ)M(z),M(z+2ϖi)=e−2πϖ(iz−ϖ)M(z)

imply that f is an elliptic function with periods 2ϖ and 2ϖi, just as sl. It follows that the function g defined by

g(z)=sl⁡zf(z),

when patched, is an elliptic function without poles. By Liouville's theorem, it is a constant. By using sl⁡z=z+O⁡(z5), M(z)=z+O⁡(z5) and N(z)=1+O⁡(z4), this constant is 1, which proves the theorem. ◼ Template:Collapse bottom Gauss conjectured that ln⁡N(ϖ)=π/2 (this later turned out to be true) and commented that this “is most remarkable and a proof of this property promises the most serious increase in analysis”.[56] Gauss expanded the products for M and N as infinite series (see below). He also discovered several identities involving the functions M and N, such as

The M function in the complex plane. The complex argument is represented by varying hue.
The N function in the complex plane. The complex argument is represented by varying hue.
N(z)=M((1+i)z)(1+i)M(z),z∉ϖℤ[i]

and

N(2z)=M(z)4+N(z)4.

Thanks to a certain theorem[57] on splitting limits, we are allowed to multiply out the infinite products and collect like powers of z. Doing so gives the following power series expansions that are convergent everywhere in the complex plane:[58][59][60][61][62]

M(z)=z−2z55!−36z99!+552z1313!+⋯,z∈ℂ
N(z)=1+2z44!−4z88!+408z1212!+⋯,z∈ℂ.

This can be contrasted with the power series of sl which has only finite radius of convergence (because it is not entire).

We define S and T by

S(z)=N(z1+i)2−iM(z1+i)2,T(z)=S(iz).

Then the lemniscate cosine can be written as

cl⁡z=S(z)T(z)

where[63]

S(z)=1−z22!−z44!−3z66!+17z88!−9z1010!+111z1212!+⋯,z∈ℂ
T(z)=1+z22!−z44!+3z66!+17z88!+9z1010!+111z1212!+⋯,z∈ℂ.

Furthermore, the identities

M(2z)=2M(z)N(z)S(z)T(z),
S(2z)=S(z)4−2M(z)4,
T(2z)=T(z)4−2M(z)4

and the Pythagorean-like identities

M(z)2+S(z)2=N(z)2,
M(z)2+N(z)2=T(z)2

hold for all z∈ℂ.

The quasi-addition formulas

M(z+w)M(z−w)=M(z)2N(w)2−N(z)2M(w)2,
N(z+w)N(z−w)=N(z)2N(w)2+M(z)2M(w)2

(where z,w∈ℂ) imply further multiplication formulas for M and N by recursion.[64]

Gauss' M and N satisfy the following system of differential equations:

M(z)M″(z)=M′(z)2−N(z)2,
N(z)N″(z)=N′(z)2+M(z)2

where z∈ℂ. Both M and N satisfy the differential equation[65]

X(z)X⁗(z)=4X′(z)X‴(z)−3X″(z)2+2X(z)2,z∈ℂ.

The functions can be also expressed by integrals involving elliptic functions:

M(z)=zexp⁡(−∫0z∫0w(1sl2v−1v2)dvdw),
N(z)=exp⁡(∫0z∫0wsl2vdvdw)

where the contours do not cross the poles; while the innermost integrals are path-independent, the outermost ones are path-dependent; however, the path dependence cancels out with the non-injectivity of the complex exponential function.

An alternative way of expressing the lemniscate functions as a ratio of entire functions involves the theta functions (see Lemniscate elliptic functions § Methods of computation); the relation between M,N and θ1,θ3 is

M(z)=2−1/4eπz2/(2ϖ2)πϖθ1(πzϖ,e−π),
N(z)=2−1/4eπz2/(2ϖ2)πϖθ3(πzϖ,e−π)

where z∈ℂ.

Relation to other functions

Relation to Weierstrass and Jacobi elliptic functions

The lemniscate functions are closely related to the Weierstrass elliptic function ℘(z;1,0) (the "lemniscatic case"), with invariants Template:Tmath and Template:Tmath. This lattice has fundamental periods ω1=2ϖ, and ω2=iω1. The associated constants of the Weierstrass function are e1=12, e2=0, e3=−12.

The related case of a Weierstrass elliptic function with Template:Tmath, Template:Tmath may be handled by a scaling transformation. However, this may involve complex numbers. If it is desired to remain within real numbers, there are two cases to consider: Template:Tmath and Template:Tmath. The period parallelogram is either a square or a rhombus. The Weierstrass elliptic function ℘(z;−1,0) is called the "pseudolemniscatic case".[66]

The square of the lemniscate sine can be represented as

sl2z=1℘(z;4,0)=i2℘((1−i)z;−1,0)=−2℘(2z+(i−1)ϖ2;1,0)

where the second and third argument of ℘ denote the lattice invariants Template:Tmath and Template:Tmath. The lemniscate sine is a rational function in the Weierstrass elliptic function and its derivative:[67]

sl⁡z=−2℘(z;−1,0)℘′(z;−1,0).

The lemniscate functions can also be written in terms of Jacobi elliptic functions. The Jacobi elliptic functions sn and cd with positive real elliptic modulus have an "upright" rectangular lattice aligned with real and imaginary axes. Alternately, the functions sn and cd with modulus Template:Tmath (and sd and cn with modulus 1/2) have a square period lattice rotated 1/8 turn.[68][69]

sl⁡z=sn⁡(z;i)=sc⁡(z;2)=12sd(2z;12)
cl⁡z=cd⁡(z;i)=dn⁡(z;2)=cn(2z;12)

where the second arguments denote the elliptic modulus k.

The functions sl~ and cl~ can also be expressed in terms of Jacobi elliptic functions:

sl~z=cd⁡(z;i)sd⁡(z;i)=dn⁡(z;2)sn⁡(z;2)=12cn⁡(2z;12)sn⁡(2z;12),
cl~z=cd⁡(z;i)nd⁡(z;i)=dn⁡(z;2)cn⁡(z;2)=cn⁡(2z;12)dn⁡(2z;12).

Relation to the modular lambda function

Template:Unreferenced section

The lemniscate sine can be used for the computation of values of the modular lambda function:

∏k=1nsl(2k−12n+1ϖ2)=λ((2n+1)i)1−λ((2n+1)i)8

For example:

sl(114ϖ)sl(314ϖ)sl(514ϖ)=λ(7i)1−λ(7i)8=tan(12arccsc(1287+21+127+1))=22+7+21+87+214+67+455+1727sl(118ϖ)sl(318ϖ)sl(518ϖ)sl(718ϖ)=λ(9i)1−λ(9i)8=tan⁡(||π4−arctan⁡(||223−23−22−33+3−1124||))

Inverse functions

The inverse function of the lemniscate sine is the lemniscate arcsine, defined as[70]

arcsl⁡x=∫0xdt1−t4.

It can also be represented by the hypergeometric function:

arcsl⁡x=x2F1(12,14;54;x4)

which can be easily seen by using the binomial series.

The inverse function of the lemniscate cosine is the lemniscate arccosine. This function is defined by following expression:

arccl⁡x=∫x1dt1−t4=12ϖ−arcsl⁡x

For Template:Tmath in the interval −1≤x≤1, sl⁡arcsl⁡x=x and cl⁡arccl⁡x=x

For the halving of the lemniscate arc length these formulas are valid:Template:Cn

sl(12arcsl⁡x)=sin(12arcsin⁡x)sech(12arsinh⁡x)sl(12arcsl⁡x)2=tan(14arcsin⁡x2)

Furthermore there are the so called Hyperbolic lemniscate area functions:Template:Cn

aslh⁡(x)=∫0x1y4+1dy=12F(2arctan⁡x;12)
aclh⁡(x)=∫x∞1y4+1dy=12F(2arccot⁡x;12)
aclh⁡(x)=ϖ2−aslh⁡(x)
aslh⁡(x)=2arcsl⁡(x/1+x4+1)
arcsl⁡(x)=2aslh⁡(x/1+1−x4)

Expression using elliptic integrals

The lemniscate arcsine and the lemniscate arccosine can also be expressed by the Legendre-Form:

These functions can be displayed directly by using the incomplete elliptic integral of the first kind:Template:Cn

arcsl⁡x=12F(arcsin2x1+x2;12)
arcsl⁡x=2(2−1)F(arcsin(2+1)x1+x2+1;(2−1)2)

The arc lengths of the lemniscate can also be expressed by only using the arc lengths of ellipses (calculated by elliptic integrals of the second kind):Template:Cn

arcsl⁡x=2+22E(arcsin(2+1)x1+x2+1;(2−1)2)  −E(arcsin2x1+x2;12)+x1−x22(1+x2+1+x2)

The lemniscate arccosine has this expression:Template:Cn

arccl⁡x=12F(arccos⁡x;12)

Use in integration

The lemniscate arcsine can be used to integrate many functions. Here is a list of important integrals (the constants of integration are omitted):

∫11−x4dx=arcsl⁡x
∫1(x2+1)(2x2+1)dx=arcslxx2+1
∫1x4+6x2+1dx=arcsl2xx4+6x2+1+x2+1
∫1x4+1dx=2arcslxx4+1+1
∫1(1−x4)34dx=2arcslx1+1−x4
∫1(x4+1)34dx=arcslxx4+14
∫1(1−x2)34dx=2arcslx1+1−x2
∫1(x2+1)34dx=2arcslxx2+1+1
∫1(ax2+bx+c)34dx=224a2c−ab24arcsl2ax+b4a(ax2+bx+c)+4ac−b2
∫sech⁡xdx=2arcsltanh⁡12x
∫sec⁡xdx=2arcsltan⁡12x

Hyperbolic lemniscate functions

Fundamental information

The hyperbolic lemniscate sine (red) and hyperbolic lemniscate cosine (purple) applied to a real argument, in comparison with the trigonometric tangent (pale dashed red).
The hyperbolic lemniscate sine in the complex plane. Dark areas represent zeros and bright areas represent poles. The complex argument is represented by varying hue.

For convenience, let σ=2ϖ. σ is the "squircular" analog of π (see below). The decimal expansion of σ (i.e. 3.7081…[71]) appears in entry 34e of chapter 11 of Ramanujan's second notebook.[72]

The hyperbolic lemniscate sine (Template:Math) and cosine (Template:Math) can be defined as inverses of elliptic integrals as follows:

z=∗∫0slh⁡zdt1+t4=∫clh⁡z∞dt1+t4

where in (∗), z is in the square with corners {σ/2,σi/2,−σ/2,−σi/2}. Beyond that square, the functions can be analytically continued to meromorphic functions in the whole complex plane.

The complete integral has the value:

∫0∞dtt4+1=14B(14,14)=σ2=1.854074677301371…

Therefore, the two defined functions have following relation to each other:

slh⁡z=clh(σ2−z)

The product of hyperbolic lemniscate sine and hyperbolic lemniscate cosine is equal to one:

slh⁡zclh⁡z=1

The functions slh and clh have a square period lattice with fundamental periods {σ,σi}.

The hyperbolic lemniscate functions can be expressed in terms of lemniscate sine and lemniscate cosine:

slh⁡(2z)=(1+cl2z)sl⁡z2cl⁡z
clh⁡(2z)=(1+sl2z)cl⁡z2sl⁡z

But there is also a relation to the Jacobi elliptic functions with the elliptic modulus one by square root of two:

slh⁡z=sn⁡(z;1/2)cd⁡(z;1/2)
clh⁡z=cd⁡(z;1/2)sn⁡(z;1/2)

The hyperbolic lemniscate sine has following imaginary relation to the lemniscate sine:

slh⁡z=1−i2sl⁡(1+i2z)=sl⁡(−14z)−14

This is analogous to the relationship between hyperbolic and trigonometric sine:

sinh⁡z=−isin⁡(iz)=sin⁡(−12z)−12

Relation to quartic Fermat curve

Hyperbolic Lemniscate Tangent and Cotangent

This image shows the standardized superelliptic Fermat squircle curve of the fourth degree:

Superellipse with the relation x4+y4=1

In a quartic Fermat curve x4+y4=1 (sometimes called a squircle) the hyperbolic lemniscate sine and cosine are analogous to the tangent and cotangent functions in a unit circle x2+y2=1 (the quadratic Fermat curve). If the origin and a point on the curve are connected to each other by a line Template:Tmath, the hyperbolic lemniscate sine of twice the enclosed area between this line and the x-axis is the y-coordinate of the intersection of Template:Tmath with the line x=1.[73] Just as π is the area enclosed by the circle x2+y2=1, the area enclosed by the squircle x4+y4=1 is σ. Moreover,

M(1,1/2)=πσ

where M is the arithmetic–geometric mean.

The hyperbolic lemniscate sine satisfies the argument addition identity:

slh⁡(a+b)=slh⁡aslh′⁡b+slh⁡bslh′⁡a1−slh2aslh2b

When u is real, the derivative and the original antiderivative of slh and clh can be expressed in this way:

dduslh⁡(u)=1+slh⁡(u)4

dduclh⁡(u)=−1+clh⁡(u)4

ddu12arsinh⁡[slh⁡(u)2]=slh⁡(u)

ddu−12arsinh⁡[clh⁡(u)2]=clh⁡(u)

There are also the Hyperbolic lemniscate tangent and the Hyperbolic lemniscate cotangent as further functions:

The functions tlh and ctlh fulfill the identities described in the differential equation mentioned:

tlh(2u)=sin4(2u)=sl⁡(u)cl2u+1sl2u+cl2u
ctlh(2u)=cos4(2u)=cl⁡(u)sl2u+1sl2u+cl2u

The functional designation sl stands for the lemniscatic sine and the designation cl stands for the lemniscatic cosine. In addition, those relations to the Jacobi elliptic functions are valid:

tlh(u)=sn(u;122)cd(u;122)4+sn(u;122)44
ctlh(u)=cd(u;122)cd(u;122)4+sn(u;122)44

When u is real, the derivative and quarter period integral of tlh and ctlh can be expressed in this way:

ddutlh⁡(u)=ctlh⁡(u)3

dductlh⁡(u)=−tlh⁡(u)3

∫0ϖ/2tlh⁡(u)du=ϖ2

∫0ϖ/2ctlh⁡(u)du=ϖ2

Derivation of the Hyperbolic Lemniscate functions

With respect to the quartic Fermat curve x4+y4=1, the hyperbolic lemniscate sine is analogous to the trigonometric tangent function. Unlike slh and clh, the functions sin4 and cos4 cannot be analytically extended to meromorphic functions in the whole complex plane.[74]

The horizontal and vertical coordinates of this superellipse are dependent on twice the enclosed area w = 2A, so the following conditions must be met:

x(w)4+y(w)4=1
ddwx(w)=−y(w)3
ddwy(w)=x(w)3
x(w=0)=1
y(w=0)=0

The solutions to this system of equations are as follows:

x(w)=cl⁡(122w)[sl⁡(122w)2+1]1/2[sl⁡(122w)2+cl⁡(122w)2]−1/2
y(w)=sl⁡(122w)[cl⁡(122w)2+1]1/2[sl⁡(122w)2+cl⁡(122w)2]−1/2

The following therefore applies to the quotient:

y(w)x(w)=sl⁡(122w)[cl⁡(122w)2+1]1/2cl⁡(122w)[sl⁡(122w)2+1]1/2=slh⁡(w)

The functions x(w) and y(w) are called cotangent hyperbolic lemniscatus and hyperbolic tangent.

x(w)=ctlh(w)
y(w)=tlh(w)

The sketch also shows the fact that the derivation of the Areasinus hyperbolic lemniscatus function is equal to the reciprocal of the square root of the successor of the fourth power function.

Specific values

Template:Unreferenced section

This list shows the values of the Hyperbolic Lemniscate Sine accurately. Recall that,

∫0∞d⁡tt4+1=14B(14,14)=ϖ2=σ2=1.85407…

whereas 12B(12,12)=π2, so the values below such as slh(ϖ22)=slh(σ4)=1 are analogous to the trigonometric sin(π2)=1.

slh⁡(ϖ22)=1
slh⁡(ϖ32)=13423−34
slh⁡(2ϖ32)=23+34
slh⁡(ϖ42)=124(2+1−1)
slh⁡(3ϖ42)=124(2+1+1)
slh⁡(ϖ52)=1845−1204−5+1=25−24sin⁡(120π)sin⁡(320π)
slh⁡(2ϖ52)=1224(5+1)204−5+1=25+24sin⁡(120π)sin⁡(320π)
slh⁡(3ϖ52)=1845−1204+5+1=25−24cos⁡(120π)cos⁡(320π)
slh⁡(4ϖ52)=1224(5+1)204+5+1=25+24cos⁡(120π)cos⁡(320π)
slh⁡(ϖ62)=12(23+3+1)(1−23−34)
slh⁡(5ϖ62)=12(23+3+1)(1+23−34)

That table shows the most important values of the Hyperbolic Lemniscate Tangent and Cotangent functions:

z clh⁡z slh⁡z ctlh⁡z=cos4z tlh⁡z=sin4z
0 ∞ 0 1 0
14σ 1 1 1/24 1/24
12σ 0 ∞ 0 1
34σ −1 −1 −1/24 1/24
σ ∞ 0 −1 0

Number theory

In algebraic number theory, every finite abelian extension of the Gaussian rationals ℚ(i) is a subfield of ℚ(i,ωn) for some positive integer n.[24][75] This is analogous to the Kronecker–Weber theorem for the rational numbers ℚ which is based on division of the circle – in particular, every finite abelian extension of ℚ is a subfield of ℚ(ζn) for some positive integer n. Both are special cases of Kronecker's Jugendtraum, which became Hilbert's twelfth problem.

The field ℚ(i,sl⁡(ϖ/n)) (for positive odd n) is the extension of ℚ(i) generated by the x- and y-coordinates of the (1+i)n-torsion points on the elliptic curve y2=4x3+x.[75]

Hurwitz numbers

The Bernoulli numbers Bn can be defined by

Bn=limz→0dndznzez−1,n≥0

and appear in

∑k∈ℤ∖{0}1k2n=(−1)n−1B2n(2π)2n(2n)!=2ζ(2n),n≥1

where ζ is the Riemann zeta function.

The Hurwitz numbers Hn, named after Adolf Hurwitz, are the "lemniscate analogs" of the Bernoulli numbers. They can be defined by[76][77]

Hn=−limz→0dndznzζ(z;1/4,0),n≥0

where ζ(⋅;1/4,0) is the Weierstrass zeta function with lattice invariants 1/4 and 0. They appear in

∑z∈ℤ[i]∖{0}1z4n=H4n(2ϖ)4n(4n)!=G4n(i),n≥1

where ℤ[i] are the Gaussian integers and G4n are the Eisenstein series of weight 4n, and in

∑n=1∞nke2πn−1={124−18πif k=1Bk+12k+2if k≡1(mod4) and k≥5Bk+12k+2+Hk+12k+2(ϖπ)k+1if k≡3(mod4) and k≥3.

The Hurwitz numbers can also be determined as follows: H4=1/10,

H4n=3(2n−3)(16n2−1)∑k=1n−1(4n4k)(4k−1)(4(n−k)−1)H4kH4(n−k),n≥2

and Hn=0 if n is not a multiple of 4.[78] This yields[76]

H8=310,H12=567130,H16=43659170,…

Also[79]

denom⁡H4n=∏(p−1)|4np

where p∈ℙ such that p≢3(mod4), just as

denom⁡B2n=∏(p−1)|2np

where p∈ℙ (by the von Staudt–Clausen theorem).

In fact, the von Staudt–Clausen theorem determines the fractional part of the Bernoulli numbers:

B2n+∑(p−1)|2n1p∈ℤ,n≥1

Template:OEIS where p is any prime, and an analogous theorem holds for the Hurwitz numbers: suppose that a∈ℤ is odd, b∈ℤ is even, p is a prime such that p≡1(mod4), p=a2+b2 (see Fermat's theorem on sums of two squares) and a≡b+1(mod4). Then for any given p, 2a=ν(p) is uniquely determined; equivalently ν(p)=p−𝒩p where 𝒩p is the number of solutions of the congruence X3−X≡Y2(mod⁡p) in variables X,Y that are non-negative integers.[80] The Hurwitz theorem then determines the fractional part of the Hurwitz numbers:[76]

H4n−12−∑(p−1)|4nν(p)4n/(p−1)p=defGn∈ℤ,n≥1.

The sequence of the integers Gn starts with 0,−1,5,253,….[76]

Let n≥2. If 4n+1 is a prime, then Gn≡1(mod4). If 4n+1 is not a prime, then Gn≡3(mod4).[81]

Some authors instead define the Hurwitz numbers as Hn′=H4n.

Appearances in Laurent series

The Hurwitz numbers appear in several Laurent series expansions related to the lemniscate functions:[82]

sl2z=∑n=1∞24n(1−(−1)n22n)H4n4nz4n−2(4n−2)!,|z|<ϖ2sl′⁡zsl⁡z=1z−∑n=1∞24n(2−(−1)n22n)H4n4nz4n−1(4n−1)!,|z|<ϖ21sl⁡z=1z−∑n=1∞22n((−1)n2−22n)H4n4nz4n−1(4n−1)!,|z|<ϖ1sl2z=1z2+∑n=1∞24nH4n4nz4n−2(4n−2)!,|z|<ϖ

Analogously, in terms of the Bernoulli numbers:

1sinh2z=1z2−∑n=1∞22nB2n2nz2n−2(2n−2)!,|z|<π.

A quartic analog of the Legendre symbol

Let p be a prime such that p≡1(mod4). A quartic residue (mod p) is any number congruent to the fourth power of an integer. Define (ap)4 to be 1 if a is a quartic residue (mod p) and define it to be −1 if a is not a quartic residue (mod p).

If a and p are coprime, then there exist numbers p′∈ℤ[i] (see[83] for these numbers) such that[84]

(ap)4=∏p′sl⁡(2ϖap′/p)sl⁡(2ϖp′/p).

This theorem is analogous to

(ap)=∏n=1p−12sin⁡(2πan/p)sin⁡(2πn/p)

where (⋅⋅) is the Legendre symbol.

World map projections

"The World on a Quincuncial Projection", from Template:Harvp.

The Peirce quincuncial projection, designed by Charles Sanders Peirce of the US Coast Survey in the 1870s, is a world map projection based on the inverse lemniscate sine of stereographically projected points (treated as complex numbers).[85]

When lines of constant real or imaginary part are projected onto the complex plane via the hyperbolic lemniscate sine, and thence stereographically projected onto the sphere (see Riemann sphere), the resulting curves are spherical conics, the spherical analog of planar ellipses and hyperbolas.[86] Thus the lemniscate functions (and more generally, the Jacobi elliptic functions) provide a parametrization for spherical conics.

A conformal map projection from the globe onto the 6 square faces of a cube can also be defined using the lemniscate functions.[87] Because many partial differential equations can be effectively solved by conformal mapping, this map from sphere to cube is convenient for atmospheric modeling.[88]

See also

Notes

Template:Reflist

References

Template:Refbegin

Template:Refend


Template:Bots

  1. ↑ Template:Harvp; Template:Harvp; Template:Harvp
  2. ↑ Template:Harvp p. 199 used the symbols Template:Math and Template:Math for the lemniscate sine and cosine, respectively, and this notation is most common today: see e.g. Template:Harvp p. 316, Template:Harvp p. 204, and Template:Harvp p. 240. Template:Harvp uses Template:Math and Template:Math. Template:Harvp use the symbols Template:Math and Template:Math. Some sources use the generic letters Template:Math and Template:Math. Template:Harvp use the letter Template:Math for the lemniscate sine and Template:Math for its derivative.
  3. ↑ The circle x2+y2=x is the unit-diameter circle centered at (12,0) with polar equation r=cos⁡θ, the degree-2 clover under the definition from Template:Harvp. This is not the unit-radius circle x2+y2=1 centered at the origin. Notice that the lemniscate (x2+y2)2=x2−y2 is the degree-4 clover.
  4. ↑ The fundamental periods (1+i)ϖ and (1−i)ϖ are "minimal" in the sense that they have the smallest absolute value of all periods whose real part is non-negative.
  5. ↑ Template:Harvp starts from this definition and thence derives other properties of the lemniscate functions.
  6. ↑ This map was the first ever picture of a Schwarz–Christoffel mapping, in Template:Harvp p. 113.
  7. ↑ Template:Harvp. OEIS sequence A062539 lists the lemniscate constant's decimal digits.
  8. ↑ Template:Cite web
  9. ↑ Template:Harvp
  10. ↑ Template:Harvp
  11. ↑ Template:Harvp
  12. ↑ Dark areas represent zeros, and bright areas represent poles. As the argument of sl⁡z changes from −π (excluding −π) to π, the colors go through cyan, blue (Arg⁡≈−π/2), magneta, red (Arg⁡≈0), orange, yellow (Arg⁡≈π/2), green, and back to cyan (Arg⁡≈π).
  13. ↑ Combining the first and fourth identity gives sl⁡z=−i/sl⁡(z−(1+i)ϖ/2). This identity is (incorrectly) given in Template:Harvp p. 226, without the minus sign at the front of the right-hand side.
  14. ↑ The even Gaussian integers are the residue class of Template:Tmath, modulo Template:Tmath, the black squares on a checkerboard.
  15. ↑ Template:Harvp; Template:Harvp
  16. ↑ 16.0 16.1 Template:Harvp
  17. ↑ Template:Harvp §22.12.6, §22.12.12
  18. ↑ Analogously, 1sin⁡z=∑n∈ℤ(−1)nz+nπ.
  19. ↑ Template:Harvp generalizes the first of these forms.
  20. ↑ Template:Harvp; Template:Harvp
  21. ↑ Template:Harvp §44 p. 79, §47 pp. 80–81
  22. ↑ 22.0 22.1 Template:Harvp §46 p. 80
  23. ↑ In fact, iε=sl⁡βϖ2.
  24. ↑ 24.0 24.1 24.2 Template:Harvp
  25. ↑ Template:Harvp
  26. ↑ The fourth root with the least positive principal argument is chosen.
  27. ↑ The restriction to positive and odd β can be dropped in deg⁡Λβ=|(𝒪/β𝒪)×|.
  28. ↑ Template:Harvp p. 142, Example 7.29(c)
  29. ↑ Template:Harvp
  30. ↑ Template:Harvp p. 200
  31. ↑ And the area enclosed by ℒ is 1, which stands in stark contrast to the unit circle (whose enclosed area is a non-constructible number).
  32. ↑ Template:Harvp; Template:Harvp. Template:Harvp use the polar-coordinate representation of the Lemniscate to derive differential arc length, but the result is the same.
  33. ↑ Template:Harvp §22.18.E6
  34. ↑ Template:Harvp; Template:Harvp
  35. ↑ Such numbers are OEIS sequence A003401.
  36. ↑ Template:Harvp; Template:Harvp; Template:Harvp
  37. ↑ Template:Harvp; Template:Harvp; Template:Harvp
  38. ↑ Template:Cite web
  39. ↑ Template:Cite book
  40. ↑ 40.0 40.1 40.2 40.3 Template:Cite web
  41. ↑ Template:Cite book p. 79, eq. 5.36
  42. ↑ Template:Cite book p. 79, eq. 5. 36 and p. 78, eq. 5.33
  43. ↑ 43.0 43.1 Template:Cite web
  44. ↑ Template:Cite book
  45. ↑ Template:Harvp §22.12.12
  46. ↑ In general, sinh⁡(x−nπ) and sin⁡(x−nπi)=−isinh⁡(ix+nπ) are not equivalent, but the resulting infinite sum is the same.
  47. ↑ Template:Harvp §22.11
  48. ↑ Template:Harvp §22.2.E7
  49. ↑ Template:Harvp p. 247, 248, 253
  50. ↑ Template:Harvp §22.11.E1
  51. ↑ Template:Harvp
  52. ↑ Template:Harvp
  53. ↑ 53.0 53.1 Template:Harvp p. 227.
  54. ↑ Template:Cite book
  55. ↑ More precisely, suppose {an} is a sequence of bounded complex functions on a set S, such that ∑|an(z)| converges uniformly on S. If {n1,n2,n3,…} is any permutation of {1,2,3,…}, then ∏n=1∞(1+an(z))=∏k=1∞(1+ank(z)) for all z∈S. The theorem in question then follows from the fact that there exists a bijection between the natural numbers and α's (resp. β's).
  56. ↑ Template:Harvp p. 58
  57. ↑ More precisely, if for each k, limn→∞ak(n) exists and there is a convergent series ∑k=1∞Mk of nonnegative real numbers such that |ak(n)|≤Mk for all n∈ℕ and 1≤k≤n, then
    limn→∞∑k=1nak(n)=∑k=1∞limn→∞ak(n).
  58. ↑ Alternatively, it can be inferred that these expansions exist just from the analyticity of M and N. However, establishing the connection to "multiplying out and collecting like powers" reveals identities between sums of reciprocals and the coefficients of the power series, like ∑α1α4=−the coefficient ofz5 in the M series, and infinitely many others.
  59. ↑ Template:Cite book p. 405; there's an error on the page: the coefficient of φ17 should be 1077410154752000, not 107207484333056000.
  60. ↑ If M(z)=∑n=0∞anzn+1, then the coefficients an are given by the recurrence an+1=−1n+1∑k=0n2n−k+1akHn−k+1(n−k+1)! with a0=1 where Hn are the Hurwitz numbers defined in Lemniscate elliptic functions § Hurwitz numbers.
  61. ↑ The power series expansions of M and N are useful for finding a β-division polynomial for the β-division of the lemniscate ℒ (where β=m+ni where m,n∈ℤ such that m+n is odd). For example, suppose we want to find a 3-division polynomial. Given that
    M(3z)=d9M(z)9+d5M(z)5N(z)4+d1M(z)N(z)8
    for some constants d1,d5,d9, from
    3z−2(3z)55!−36(3z)99!+O⁡(z13)=d9x9+d5x5y4+d1xy8,
    where
    x=z−2z55!−36z99!+O⁡(z13),y=1+2z44!−4z88!+O⁡(z12),
    we have
    {d1,d5,d9}={3,−6,−1}.
    Therefore, a 3-division polynomial is
    −X9−6X5+3X
    (meaning one of its roots is sl⁡(2ϖ/3)). The equations arrived at by this process are the lemniscate analogs of
    Xn=1
    (so that e2πi/n is one of the solutions) which comes up when dividing the unit circle into n arcs of equal length. In the following note, the first few coefficients of the monic normalization of such β-division polynomials are described symbolically in terms of β.
  62. ↑ By utilizing the power series expansion of the N function, it can be proved that a polynomial having sl⁡(2ϖ/β) as one of its roots (with β from the previous note) is
    ∑n=0(ββ‾−1)/4a4n+1(β)Xββ‾−4n
    where
    a1(β)=1,a5(β)=β4−ββ‾12,a9(β)=−β8−70β5β‾+336β4+35β2β‾2−300ββ‾10080
    and so on.
  63. ↑ Template:Cite book
  64. ↑ For example, by the quasi-addition formulas, the duplication formulas and the Pythagorean-like identities, we have
    M(3z)=−M(z)9−6M(z)5N(z)4+3M(z)N(z)8,
    N(3z)=N(z)9+6M(z)4N(z)5−3M(z)8N(z),
    so
    sl⁡3z=−M(z)9−6M(z)5N(z)4+3M(z)N(z)8N(z)9+6M(z)4N(z)5−3M(z)8N(z).
    On dividing the numerator and the denominator by N(z)9, we obtain the triplication formula for sl:
    sl⁡3z=−sl9z−6sl5z+3sl⁡z1+6sl4z−3sl8z.
  65. ↑ Gauss (1866), p. 408
  66. ↑ Template:Harvp
  67. ↑ Template:Harvp p. 234
  68. ↑ Template:Cite book
  69. ↑ The identity cl⁡z=cn(2z;12) can be found in Template:Harvp p. 33.
  70. ↑ Template:Harvp
  71. ↑ Template:Cite OEIS
  72. ↑ Template:Cite book p. 96
  73. ↑ Template:Harvp; Template:Harvp
  74. ↑ Template:Harvp p. 515
  75. ↑ 75.0 75.1 Template:Harvp p. 508, 509
  76. ↑ 76.0 76.1 76.2 76.3 Template:Cite book p. 203—206
  77. ↑ Equivalently, Hn=−limz→0dndzn((1+i)z/2sl⁡((1+i)z/2)+z2ℰ(z2;i)) where n≥4 and ℰ(⋅;i) is the Jacobi epsilon function with modulus i.
  78. ↑ The Bernoulli numbers can be determined by an analogous recurrence: B2n=−12n+1∑k=1n−1(2n2k)B2kB2(n−k) where n≥2 and B2=1/6.
  79. ↑ Template:Cite journal See eq. (9)
  80. ↑ For more on the ν function, see Lemniscate constant.
  81. ↑ Template:Cite book p. 370
  82. ↑ Arakawa et al. (2014) define H4n by the expansion of 1/sl2.
  83. ↑ Template:Cite journal Eisenstein uses φ=sl and ω=2ϖ.
  84. ↑ Template:Harvp
  85. ↑ Template:Harvp. Template:Harvp and Template:Harvp introduced transverse and oblique aspects of the same projection, respectively. Also see Template:Harvp. These authors write their projection formulas in terms of Jacobi elliptic functions, with a square lattice.
  86. ↑ Template:Harvp
  87. ↑ Template:Harvp; Template:Harvp.
  88. ↑ Template:Harvp; Template:Harvp.