Lemniscate elliptic functions – Wikipedia

Mathematical functions

The lemniscate sine (red) and lemniscate cosine (purple) applied to a real argument, in comparison with the trigonometric sine y = sin(πx/ϖ) (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 sl and cl (sometimes the symbols sinlem and coslem or sin lemn and cos lemn 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,{displaystyle x^{2}+y^{2}=x,}

[3] the lemniscate sine relates the arc length to the chord length of a lemniscate

(x2+y2)2=x2y2.{displaystyle {bigl (}x^{2}+y^{2}{bigr )}{}^{2}=x^{2}-y^{2}.}

The lemniscate functions have periods related to a number

ϖ={displaystyle varpi =}

2.622057… called the lemniscate constant, the ratio of a lemniscate’s perimeter to its diameter. This number is a quartic analog of the (quadratic)

π={displaystyle pi =}

3.141592…, ratio of perimeter to diameter of a circle.

As complex functions, sl and cl have a square period lattice (a multiple of the Gaussian integers) with fundamental periods

{(1+i)ϖ,(1i)ϖ},{displaystyle {(1+i)varpi ,(1-i)varpi },}

[4] and are a special case of two Jacobi elliptic functions on that lattice,

slz=sn(z;i),{displaystyle operatorname {sl} z=operatorname {sn} (z;i),}

clz=cd(z;i){displaystyle operatorname {cl} z=operatorname {cd} (z;i)}

.

Similarly, the hyperbolic lemniscate sine slh and hyperbolic lemniscate cosine clh have a square period lattice with fundamental periods

{2ϖ,2ϖi}.{displaystyle {bigl {}{sqrt {2}}varpi ,{sqrt {2}}varpi i{bigr }}.}

The lemniscate functions and the hyperbolic lemniscate functions are related to the Weierstrass elliptic function

(z;a,0){displaystyle wp (z;a,0)}

.

Lemniscate sine and cosine functions[edit]

Definitions[edit]

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

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}:{displaystyle {big {}{tfrac {1}{2}}varpi ,{tfrac {1}{2}}varpi i,-{tfrac {1}{2}}varpi ,-{tfrac {1}{2}}varpi i{big }}colon }

[6]

Beyond that square, the functions can be analytically continued to the whole complex plane by a series of reflections.

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

or as inverses of a map from the upper half-plane to a half-infinite strip with real part between

12π,12π{displaystyle -{tfrac {1}{2}}pi ,{tfrac {1}{2}}pi }

and positive imaginary part:

Arc length of Bernoulli’s lemniscate[edit]

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 half-width 1 is the locus of points in the plane such that the product of their distances from the two focal points

F1=(12,0){displaystyle F_{1}={bigl (}{-{tfrac {1}{sqrt {2}}}},0{bigr )}}

and

F2=(12,0){displaystyle F_{2}={bigl (}{tfrac {1}{sqrt {2}}},0{bigr )}}

is the constant

12{displaystyle {tfrac {1}{2}}}

. This is a quartic curve satisfying the polar equation

r2=cos2θ{displaystyle r^{2}=cos 2theta }

or the Cartesian equation

(x2+y2)2=x2y2.{displaystyle {bigl (}x^{2}+y^{2}{bigr )}{}^{2}=x^{2}-y^{2}.}

The points on the lemniscate at distance

r{displaystyle r}

from the origin are the intersections of the circle

x2+y2=r2{displaystyle x^{2}+y^{2}=r^{2}}

and the hyperbola

x2y2=r4{displaystyle x^{2}-y^{2}=r^{4}}

. The intersection in the positive quadrant has Cartesian coordinates:

Using this parametrization with

r[0,1]{displaystyle rin [0,1]}

for a quarter of the lemniscate, the arc length from the origin to a point

(x(r),y(r)){displaystyle {big (}x(r),y(r){big )}}

is:[7]

Likewise, the arc length from

(1,0){displaystyle (1,0)}

to

(x(r),y(r)){displaystyle {big (}x(r),y(r){big )}}

is:

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){displaystyle (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θ{displaystyle r=cos theta }

or Cartesian equation

x2+y2=x,{displaystyle x^{2}+y^{2}=x,}

using the same argument above but with the parametrization:

Alternatively, just as the unit circle

x2+y2=1{displaystyle x^{2}+y^{2}=1}

is parametrized in terms of the arc length

s{displaystyle s}

from the point

(1,0){displaystyle (1,0)}

by

the lemniscate is parametrized in terms of the arc length

s{displaystyle s}

from the point

(1,0){displaystyle (1,0)}

by[8]

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

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 n sections of equal arc length using only straightedge and compass if and only if n is of the form

n=2kp1p2pm{displaystyle n=2^{k}p_{1}p_{2}cdots p_{m}}

where k is a non-negative integer and each pi (if any) is a distinct Fermat prime.[10] 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.[11] Equivalently, the lemniscate can be divided into n sections of equal arc length using only straightedge and compass if and only if

φ(n){displaystyle varphi (n)}

is a power of two (where

φ{displaystyle varphi }

is Euler’s totient function). The lemniscate is not assumed to be already drawn; the theorem refers to constructing the division points only.

Let

rj=sl2jϖn{displaystyle r_{j}=operatorname {sl} {dfrac {2jvarpi }{n}}}

. Then the n-division points for the lemniscate

(x2+y2)2=x2y2{displaystyle (x^{2}+y^{2})^{2}=x^{2}-y^{2}}

are the points

where

{displaystyle lfloor cdot rfloor }

is the floor function. See below for some specific values of

sl2ϖn{displaystyle operatorname {sl} {dfrac {2varpi }{n}}}

.

Arc length of rectangular elastica[edit]

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

The inverse lemniscate sine also describes the arc length s relative to the x coordinate of the rectangular elastica.[12] This curve has y coordinate and arc length:

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[edit]

The lemniscate elliptic functions and an ellipse

Let

C{displaystyle C}

be a point on the ellipse

x2+2y2=1{displaystyle x^{2}+2y^{2}=1}

in the first quadrant and let

D{displaystyle D}

be the projection of

C{displaystyle C}

on the unit circle

x2+y2=1{displaystyle x^{2}+y^{2}=1}

. The distance

r{displaystyle r}

between the origin

A{displaystyle A}

and the point

C{displaystyle C}

is a function of

φ{displaystyle varphi }

(the angle

BAC{displaystyle BAC}

where

B=(1,0){displaystyle B=(1,0)}

; equivalently the length of the circular arc

BD{displaystyle BD}

). The parameter

u{displaystyle u}

is given by

If

E{displaystyle E}

is the projection of

D{displaystyle D}

on the x-axis and if

F{displaystyle F}

is the projection of

C{displaystyle C}

on the x-axis, then the lemniscate elliptic functions are given by

Relation to the lemniscate constant[edit]

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 ϖ.

The lemniscate functions have minimal real period 2ϖ[13] and fundamental complex periods

(1+i)ϖ{displaystyle (1+i)varpi }

and

(1i)ϖ{displaystyle (1-i)varpi }

for a constant ϖ called the lemniscate constant,[14]

The lemniscate functions satisfy the basic relation

clz=sl(12ϖz),{displaystyle operatorname {cl} z={operatorname {sl} }{bigl (}{tfrac {1}{2}}varpi -z{bigr )},}

analogous to the relation

cosz=sin(12πz).{displaystyle cos z={sin }{bigl (}{tfrac {1}{2}}pi -z{bigr )}.}

The lemniscate constant ϖ is a close analog of the circle constant π, and many identities involving π have analogues involving ϖ, as identities involving the trigonometric functions have analogues involving the lemniscate functions. For example, Viète’s formula for π can be written:

An analogous formula for ϖ is:[15]

The Machin formula for π is

14π=4arctan15arctan1239,{textstyle {tfrac {1}{4}}pi =4arctan {tfrac {1}{5}}-arctan {tfrac {1}{239}},}

and several similar formulas for π can be developed using trigonometric angle sum identities, e.g. Euler’s formula

14π=arctan12+arctan13{textstyle {tfrac {1}{4}}pi =arctan {tfrac {1}{2}}+arctan {tfrac {1}{3}}}

. Analogous formulas can be developed for ϖ, including the following found by Gauss:

12ϖ=2arcsl12+arcsl723.{displaystyle {tfrac {1}{2}}varpi =2operatorname {arcsl} {tfrac {1}{2}}+operatorname {arcsl} {tfrac {7}{23}}.}

[16]

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

Zeros, poles and symmetries[edit]

The lemniscate functions cl and sl are even and odd functions, respectively,

At translations of

12ϖ,{displaystyle {tfrac {1}{2}}varpi ,}

cl and sl are exchanged, and at translations of

12iϖ{displaystyle {tfrac {1}{2}}ivarpi }

they are additionally rotated and reciprocated:[19]

Doubling these to translations by a unit-Gaussian-integer multiple of

ϖ{displaystyle varpi }

(that is,

±ϖ{displaystyle pm varpi }

or

±iϖ{displaystyle pm ivarpi }

), negates each function, an involution:

As a result, both functions are invariant under translation by an even-Gaussian-integer multiple of

ϖ{displaystyle varpi }

.[20] That is, a displacement

(a+bi)ϖ,{displaystyle (a+bi)varpi ,}

with

a+b=2k{displaystyle a+b=2k}

for integers a, b, and k.

This makes them elliptic functions (doubly periodic meromorphic functions in the complex plane) with a diagonal square period lattice of fundamental periods

(1+i)ϖ{displaystyle (1+i)varpi }

and

(1i)ϖ{displaystyle (1-i)varpi }

.[21] 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:

The sl function has simple zeros at Gaussian integer multiples of ϖ, complex numbers of the form

aϖ+bϖi{displaystyle avarpi +bvarpi i}

for integers a and b. It has simple poles at Gaussian half-integer multiples of ϖ, complex numbers of the form

(a+12)ϖ+(b+12)ϖi{displaystyle {bigl (}a+{tfrac {1}{2}}{bigr )}varpi +{bigl (}b+{tfrac {1}{2}}{bigr )}varpi i}

, with residues

(1)ab+1i{displaystyle (-1)^{a-b+1}i}

. The cl function is reflected and offset from the sl function,

clz=sl(12ϖz){displaystyle operatorname {cl} z={operatorname {sl} }{bigl (}{tfrac {1}{2}}varpi -z{bigr )}}

. It has zeros for arguments

(a+12)ϖ+bϖi{displaystyle {bigl (}a+{tfrac {1}{2}}{bigr )}varpi +bvarpi i}

and poles for arguments

aϖ+(b+12)ϖi,{displaystyle avarpi +{bigl (}b+{tfrac {1}{2}}{bigr )}varpi i,}

with residues

(1)abi.{displaystyle (-1)^{a-b}i.}

Also

for some

m,nZ{displaystyle m,nin mathbb {Z} }

and

The last formula is a special case of complex multiplication. Analogous formulas can be given for

sl((n+mi)z){displaystyle operatorname {sl} ((n+mi)z)}

where

n+mi{displaystyle n+mi}

is any Gaussian integer – the function

sl{displaystyle operatorname {sl} }

has complex multiplication by

Z[i]{displaystyle mathbb {Z} [i]}

.[22]

There are also infinite series reflecting the distribution of the zeros and poles of sl:[23][24]

The lemniscate sine as a ratio of entire functions[edit]

The

The

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 sl has the following product expansion, reflecting the distribution of its zeros and poles:[25]

where

Here,

α{displaystyle alpha }

and

β{displaystyle beta }

denote, respectively, the zeros and poles of sl which are in the quadrant

Rez>0,Imz0{displaystyle operatorname {Re} z>0,operatorname {Im} zgeq 0}

lnN(ϖ)=π/2{displaystyle ln N(varpi )=pi /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”.[26] Gauss expanded the products for

M{displaystyle M}

and

N{displaystyle N}

as infinite series. He also discovered several identities involving the functions

M{displaystyle M}

and

N{displaystyle N}

, such as

and

Since the functions

M{displaystyle M}

and

N{displaystyle N}

are entire, their power series expansions converge everywhere in the complex plane:[27][28][29]

Pythagorean-like identity[edit]

Curves x² ⊕ y² = a for various values of a. Negative a in green, positive a in blue, a = ±1 in red, a = ∞ in black.

The lemniscate functions satisfy a Pythagorean-like identity:

As a result, the parametric equation

(x,y)=(clt,slt){displaystyle (x,y)=(operatorname {cl} t,operatorname {sl} t)}

parametrizes the quartic curve

x2+y2+x2y2=1.{displaystyle x^{2}+y^{2}+x^{2}y^{2}=1.}

This identity can alternately be rewritten:[30]

Defining a tangent-sum operator as

ab:=tan(arctana+arctanb),{displaystyle aoplus bmathrel {:=} tan(arctan a+arctan b),}

gives:

The functions

cl~{displaystyle {tilde {operatorname {cl} }}}

and

sl~{displaystyle {tilde {operatorname {sl} }}}

satisfy another Pythagorean-like identity:

Derivatives and integrals[edit]

The derivatives are as follows:

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

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

Argument sum and multiple identities[edit]

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:[31]

The derivative and Pythagorean-like identities can be used to rework the identity used by Fagano in terms of sl and cl. Defining a tangent-sum operator

ab:=tan(arctana+arctanb){displaystyle aoplus bmathrel {:=} tan(arctan a+arctan b)}

and tangent-difference operator

ab:=a(b),{displaystyle aominus bmathrel {:=} aoplus (-b),}

the argument sum and difference identities can be expressed as:[32]

These resemble their trigonometric analogs:

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

Bisection formulas:

Duplication formulas:[33]

Triplication formulas:[33]

Note the “reverse symmetry” of the coefficients of numerator and denominator of

sl3x{displaystyle operatorname {sl} 3x}

. This phenomenon can be observed in multiplication formulas for

slβx{displaystyle operatorname {sl} beta x}

where

β=m+ni{displaystyle beta =m+ni}

whenever

m,nZ{displaystyle m,nin mathbb {Z} }

and

m+n{displaystyle m+n}

is odd.[22]

Lemnatomic polynomials[edit]

Let

L{displaystyle L}

be the lattice

Furthermore, let

K=Q(i){displaystyle K=mathbb {Q} (i)}

,

O=Z[i]{displaystyle {mathcal {O}}=mathbb {Z} [i]}

,

zC{displaystyle zin mathbb {C} }

,

β=m+in{displaystyle beta =m+in}

,

γ=m+in{displaystyle gamma =m’+in’}

(where

m,n,m,nZ{displaystyle m,n,m’,n’in mathbb {Z} }

),

m+n{displaystyle m+n}

be odd,

m+n{displaystyle m’+n’}

be odd,

γ1mod2(1+i){displaystyle gamma equiv 1,operatorname {mod} ,2(1+i)}

and

slβz=Mβ(slz){displaystyle operatorname {sl} beta z=M_{beta }(operatorname {sl} z)}

. Then

for some coprime polynomials

Pβ(x),Qβ(x)O[x]{displaystyle P_{beta }(x),Q_{beta }(x)in {mathcal {O}}[x]}


and some

ε{0,1,2,3}{displaystyle varepsilon in {0,1,2,3}}

[34] where

and

where

δβ{displaystyle delta _{beta }}

is any

β{displaystyle beta }

-torsion generator (i.e.

δβ(1/β)L{displaystyle delta _{beta }in (1/beta )L}

and

[δβ](1/β)L/L{displaystyle [delta _{beta }]in (1/beta )L/L}

generates

(1/β)L/L{displaystyle (1/beta )L/L}

as an

O{displaystyle {mathcal {O}}}

-module). Examples of

β{displaystyle beta }

-torsion generators include

2ϖ/β{displaystyle 2varpi /beta }

and

(1+i)ϖ/β{displaystyle (1+i)varpi /beta }

. The polynomial

Λβ(x)O[x]{displaystyle Lambda _{beta }(x)in {mathcal {O}}[x]}

is called the

β{displaystyle beta }

-th lemnatomic polynomial. It is monic and is irreducible over

K{displaystyle K}

. The lemnatomic polynomials are the “lemniscate analogs” of the cyclotomic polynomials,[35]

The

β{displaystyle beta }

-th lemnatomic polynomial

Λβ(x){displaystyle Lambda _{beta }(x)}

is the minimal polynomial of

slδβ{displaystyle operatorname {sl} delta _{beta }}

in

K[x]{displaystyle K[x]}

. For convenience, let

ωβ=sl(2ϖ/β){displaystyle omega _{beta }=operatorname {sl} (2varpi /beta )}

and

ω~β=sl((1+i)ϖ/β){displaystyle {tilde {omega }}_{beta }=operatorname {sl} ((1+i)varpi /beta )}

. So for example, the minimal polynomial of

ω5{displaystyle omega _{5}}

(and also of

ω~5{displaystyle {tilde {omega }}_{5}}

) in

K[x]{displaystyle K[x]}

is

and[36]

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

which is the minimal polynomial of

ω1+2i{displaystyle omega _{-1+2i}}

(and also of

ω~1+2i{displaystyle {tilde {omega }}_{-1+2i}}

) in

K[x].{displaystyle K[x].}

If

p{displaystyle p}

is prime and

β{displaystyle beta }

is positive and odd,[38] then[39]

which can be compared to the cyclotomic analog

Specific values[edit]

Just as for the trigonometric functions, values of the lemniscate functions can be computed for divisions of the lemniscate into n parts of equal length, using only basic arithmetic and square roots, if and only if n is of the form

n=2kp1p2pm{displaystyle n=2^{k}p_{1}p_{2}cdots p_{m}}

where k is a non-negative integer and each pi (if any) is a distinct Fermat prime.[40] The expressions become unwieldy as n grows. Below are the expressions for dividing the lemniscate

(x2+y2)2=x2y2{displaystyle (x^{2}+y^{2})^{2}=x^{2}-y^{2}}

into n parts of equal length for some n ≤ 20.

Power series[edit]

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

where the coefficients

an{displaystyle a_{n}}

are determined as follows:

where

i+j+k=n{displaystyle i+j+k=n}

stands for all three-term compositions of

n{displaystyle n}

. For example, to evaluate

a13{displaystyle a_{13}}

, it can be seen that there are only six compositions of

132=11{displaystyle 13-2=11}

that give a nonzero contribution to the sum:

11=9+1+1=1+9+1=1+1+9{displaystyle 11=9+1+1=1+9+1=1+1+9}

and

11=5+5+1=5+1+5=1+5+5{displaystyle 11=5+5+1=5+1+5=1+5+5}

, so

The expansion can be equivalently written as[42]

where

The power series expansion of

sl~{displaystyle {tilde {operatorname {sl} }}}

at the origin is

where

αn=0{displaystyle alpha _{n}=0}

if

n{displaystyle n}

is even and[43]

if

n{displaystyle n}

is odd.

The expansion can be equivalently written as[44]

where

For the lemniscate cosine,[45]

where

Relation to Weierstrass and Jacobi elliptic functions[edit]

The lemniscate functions are closely related to the Weierstrass elliptic function

(z;1,0){displaystyle wp (z;1,0)}

(the “lemniscatic case”), with invariants g2 = 1 and g3 = 0. This lattice has fundamental periods

ω1=2ϖ,{displaystyle omega _{1}={sqrt {2}}varpi ,}

and

ω2=iω1{displaystyle omega _{2}=iomega _{1}}

. The associated constants of the Weierstrass function are

e1=12, e2=0, e3=12.{displaystyle e_{1}={tfrac {1}{2}}, e_{2}=0, e_{3}=-{tfrac {1}{2}}.}

The related case of a Weierstrass elliptic function with g2 = a, g3 = 0 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: a > 0 and a < 0. The period parallelogram is either a square or a rhombus. The Weierstrass elliptic function

(z;1,0){displaystyle wp (z;-1,0)}

is called the “pseudolemniscatic case”.[46]

The square of the lemniscate sine can be represented as

where the second and third argument of

{displaystyle wp }

denote the lattice invariants g2 and g3. Another representation is

where the second argument of

{displaystyle wp }

denotes the period ratio

τ{displaystyle tau }

.[47] The lemniscate sine is a rational function in the Weierstrass elliptic function and its derivative:[48]

where the second and third argument of

{displaystyle wp }

denote the lattice invariants g2 and g3. In terms of the period ratio

τ{displaystyle tau }

, this becomes

The lemniscate functions can also be written in terms of Jacobi elliptic functions. The Jacobi elliptic functions

sn{displaystyle operatorname {sn} }

and

cd{displaystyle operatorname {cd} }

with positive real elliptic modulus have an “upright” rectangular lattice aligned with real and imaginary axes. Alternately, the functions

sn{displaystyle operatorname {sn} }

and

cd{displaystyle operatorname {cd} }

with modulus i (and

sd{displaystyle operatorname {sd} }

and

cn{displaystyle operatorname {cn} }

with modulus

1/2{displaystyle 1/{sqrt {2}}}

) have a square period lattice rotated 1/8 turn.[49][50]

where the second arguments denote the elliptic modulus

k{displaystyle k}

.

The functions

sl~{displaystyle {tilde {operatorname {sl} }}}

and

cl~{displaystyle {tilde {operatorname {cl} }}}

can also be expressed in terms of Jacobi elliptic functions:

Relation to the modular lambda function[edit]

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

For example:

Ramanujan’s cos/cosh identity[edit]

Ramanujan’s famous cos/cosh identity states that if

then[43]

There is a close relation between the lemniscate functions and

R(s){displaystyle R(s)}

. Indeed,[43][51]

and

Continued fractions[edit]

For

zC{0}{displaystyle zin mathbb {C} setminus {0}}

:[52]

Methods of computation[edit]

A fast algorithm, returning approximations to

slx{displaystyle operatorname {sl} x}

(which get closer to

slx{displaystyle operatorname {sl} x}

with increasing

N{displaystyle N}

), is the following:[53]

This is effectively using the arithmetic-geometric mean and is based on Landen’s transformations.[54]

Several methods of computing

slx{displaystyle operatorname {sl} x}

involve first making the change of variables

πx=ϖx~{displaystyle pi x=varpi {tilde {x}}}

and then computing

sl(ϖx~/π).{displaystyle operatorname {sl} (varpi {tilde {x}}/pi ).}

A hyperbolic series method:[55][56][57]

Fourier series method:[58]

The lemniscate functions can be computed more rapidly by

where

are the Jacobi theta functions.[59]

Two other fast computation methods use the following sum and product series:

where

f(x)=tan(2arctanx)=2x/(1x2).{displaystyle f(x)=tan(2arctan x)=2x/(1-x^{2}).}

Fourier series for the logarithm of the lemniscate sine:

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

The functions

sl~{displaystyle {tilde {operatorname {sl} }}}

and

cl~{displaystyle {tilde {operatorname {cl} }}}

analogous to

sin{displaystyle sin }

and

cos{displaystyle cos }

on the unit circle have the following Fourier and hyperbolic series expansions:[43][51][61]

Inverse functions[edit]

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

It can also be represented by the hypergeometric function:

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

For x in the interval

1x1{displaystyle -1leq xleq 1}

,

slarcslx=x{displaystyle operatorname {sl} operatorname {arcsl} x=x}

and

clarcclx=x{displaystyle operatorname {cl} operatorname {arccl} x=x}

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

Expression using elliptic integrals[edit]

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:

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):

The lemniscate arccosine has this expression:

Use in integration[edit]

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

Hyperbolic lemniscate functions[edit]

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ϖ{displaystyle sigma ={sqrt {2}}varpi }

.

σ{displaystyle sigma }

is the “squircular” analog of

π{displaystyle pi }

(see below). The decimal expansion of

σ{displaystyle sigma }

(i.e.

3.7081{displaystyle 3.7081ldots }

[62]) appears in entry 34e of chapter 11 of Ramanujan’s second notebook.[63]

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

where in

(){displaystyle (*)}

,

z{displaystyle z}

is in the square with corners

{σ/2,σi/2,σ/2,σi/2}{displaystyle {sigma /2,sigma i/2,-sigma /2,-sigma 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:

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

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

The functions

slh{displaystyle operatorname {slh} }

and

clh{displaystyle operatorname {clh} }

have a square period lattice with fundamental periods

{σ,σi}{displaystyle {sigma ,sigma i}}

.

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

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

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

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

In a quartic Fermat curve

x4+y4=1{displaystyle x^{4}+y^{4}=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{displaystyle x^{2}+y^{2}=1}

(the quadratic Fermat curve). If the origin and a point on the curve are connected to each other by a line L, the hyperbolic lemniscate sine of twice the enclosed area between this line and the x-axis is the y-coordinate of the intersection of L with the line

x=1{displaystyle x=1}

.[65] Just as

π{displaystyle pi }

is the area enclosed by the circle

x2+y2=1{displaystyle x^{2}+y^{2}=1}

, the area enclosed by the squircle

x4+y4=1{displaystyle x^{4}+y^{4}=1}

is

σ{displaystyle sigma }

. Moreover,

where

M{displaystyle M}

is the arithmetic–geometric mean.

The hyperbolic lemniscate sine satisfies the argument addition identity:

When

x{displaystyle x}

is real, the derivative can be expressed in this way:

Number theory[edit]

In algebraic number theory, every finite abelian extension of the Gaussian rationals

Q(i){displaystyle mathbb {Q} (i)}

is a subfield of

Q(i,ωn){displaystyle mathbb {Q} (i,omega _{n})}

for some positive integer

n{displaystyle n}

.[35][66] This is analogous to the Kronecker–Weber theorem for the rational numbers

Q{displaystyle mathbb {Q} }

which is based on division of the circle – in particular, every finite abelian extension of

Q{displaystyle mathbb {Q} }

is a subfield of

Q(ζn){displaystyle mathbb {Q} (zeta _{n})}

for some positive integer

n{displaystyle n}

. Both are special cases of Kronecker’s Jugendtraum, which became Hilbert’s twelfth problem.

The field

Q(i,sl(ϖ/n)){displaystyle mathbb {Q} (i,operatorname {sl} (varpi /n))}

(for positive odd

n{displaystyle n}

) is the extension of

Q(i){displaystyle mathbb {Q} (i)}

generated by the

x{displaystyle x}

– and

y{displaystyle y}

-coordinates of the

(1+i)n{displaystyle (1+i)n}

-torsion points on the elliptic curve

y2=4x3+x{displaystyle y^{2}=4x^{3}+x}

.[66]

Hurwitz numbers[edit]

The Bernoulli numbers

Bn{displaystyle mathrm {B} _{n}}

can be defined by

and appear in

where

ζ{displaystyle zeta }

is the Riemann zeta function.

The Hurwitz numbers

Hn,{displaystyle mathrm {H} _{n},}

named after Adolf Hurwitz, are the “lemniscate analogs” of the Bernoulli numbers. They can be defined by[67][68]

where

ζ(;1/4,0){displaystyle zeta (cdot ;1/4,0)}

is the Weierstrass zeta function with lattice invariants

1/4{displaystyle 1/4}

and

0{displaystyle 0}

. They appear in

where

Z[i]{displaystyle mathbb {Z} [i]}

are the Gaussian integers and

G4n{displaystyle G_{4n}}

are the Eisenstein series of weight

4n{displaystyle 4n}

, and in

The Hurwitz numbers can also be determined as follows:

H4=1/10{displaystyle mathrm {H} _{4}=1/10}

,

and

Hn=0{displaystyle mathrm {H} _{n}=0}

if

n{displaystyle n}

is not a multiple of

4{displaystyle 4}

.[69] This yields[67]

Also[70]

where

pP{displaystyle pin mathbb {P} }

such that

p1(mod4),{displaystyle pequiv 1,({text{mod}},4),}


just as

where

pP{displaystyle pin mathbb {P} }

(by the von Staudt–Clausen theorem).

In fact, the von Staudt–Clausen theorem states that

(sequence A000146 in the OEIS) where

p{displaystyle p}

is any prime, and an analogous theorem holds for the Hurwitz numbers: suppose that

aZ{displaystyle ain mathbb {Z} }

is odd,

bZ{displaystyle bin mathbb {Z} }

is even,

p{displaystyle p}

is a prime such that

p1(mod4){displaystyle pequiv 1,(mathrm {mod} ,4)}

,

p=a2+b2{displaystyle p=a^{2}+b^{2}}

(see Fermat’s theorem on sums of two squares) and

ab+1(mod4){displaystyle aequiv b+1,(mathrm {mod} ,4)}

. Then for any given

p{displaystyle p}

,

a=ap{displaystyle a=a_{p}}

is uniquely determined and[67]

The sequence of the integers

Gn{displaystyle mathrm {G} _{n}}

starts with

0,1,5,253,.{displaystyle 0,-1,5,253,ldots .}

[67]

Let

n2{displaystyle ngeq 2}

. If

4n+1{displaystyle 4n+1}

is a prime, then

Gn1(mod4){displaystyle mathrm {G} _{n}equiv 1,(mathrm {mod} ,4)}

. If

4n+1{displaystyle 4n+1}

is not a prime, then

Gn3(mod4){displaystyle mathrm {G} _{n}equiv 3,(mathrm {mod} ,4)}

.[71]

Some authors instead define the Hurwitz numbers as

Hn=H4n{displaystyle mathrm {H} _{n}’=mathrm {H} _{4n}}

.

Appearances in Laurent series[edit]

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

Analogously, in terms of the Bernoulli numbers:

World map projections[edit]

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).[73]

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.[74] 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.[75] Because many partial differential equations can be effectively solved by conformal mapping, this map from sphere to cube is convenient for atmospheric modeling.[76]

See also[edit]

  1. ^ Fagnano (1718–1723); Euler (1761); Gauss (1917)
  2. ^ Gauss (1917) p. 199 used the symbols sl and cl for the lemniscate sine and cosine, respectively, and this notation is most common today: see e.g. Cox (1984) p. 316, Eymard & Lafon (2004) p. 204, and Lemmermeyer (2000) p. 240. Ayoub (1984) uses sinlem and coslem. Whittaker & Watson (1920) use the symbols sin lemn and cos lemn. Some sources use the generic letters s and c. Prasolov & Solovyev (1997) use the letter φ for the lemniscate sine and φ′ for its derivative.
  3. ^ The circle
  4. ^ The fundamental periods
  5. ^ Robinson (2019a) 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 Schwarz (1869) p. 113.
  7. ^ Euler (1761); Siegel (1969). Prasolov & Solovyev (1997) use the polar-coordinate representation of the Lemniscate to derive differential arc length, but the result is the same.
  8. ^ Reinhardt & Walker (2010a) §22.18.E6
  9. ^ Siegel (1969); Schappacher (1997)
  10. ^ Such numbers are OEIS sequence A003401.
  11. ^ Abel (1827–1828); Rosen (1981); Prasolov & Solovyev (1997)
  12. ^ Euler (1786); Sridharan (2004); Levien (2008)
  13. ^ 2ϖi is also a period.
  14. ^ Schappacher (1997). OEIS sequence A062539 lists the lemniscate constant’s decimal digits.
  15. ^ Levin (2006)
  16. ^ Todd (1975)
  17. ^ Cox (1984)
  18. ^ Dark areas represent zeros, and bright areas represent poles. As the argument of
  19. ^ Combining the first and fourth identity gives
  20. ^ The even Gaussian integers are the residue class of 0, modulo 1 + i, the black squares on a checkerboard.
  21. ^ Prasolov & Solovyev (1997); Robinson (2019a)
  22. ^ a b Cox (2012)
  23. ^ Reinhardt & Walker (2010a) §22.12.6, §22.12.12
  24. ^ Analogously,
  25. ^ Eymard & Lafon (2004) p. 227.
  26. ^ Bottazzini & Gray (2013) p. 58
  27. ^ Gauss, C. F. (1866). Werke (Band III) (in Latin and German). Herausgegeben der Königlichen Gesellschaft der Wissenschaften zu Göttingen. p. 405; there’s an error on the page: the coefficient of
  28. ^ The
  29. ^ If
  30. ^ Lindqvist & Peetre (2001) generalizes the first of these forms.
  31. ^ Ayoub (1984); Prasolov & Solovyev (1997)
  32. ^ Euler (1761) §44 p. 79, §47 pp. 80–81
  33. ^ a b Euler (1761) §46 p. 80
  34. ^ In fact,
  35. ^ a b c Cox & Hyde (2014)
  36. ^ Gómez-Molleda & Lario (2019)
  37. ^ a b The fourth root with the least positive principal argument is chosen.
  38. ^ The restriction to positive and odd
  39. ^ Cox (2013) p. 142, Example 7.29(c)
  40. ^ Rosen (1981)
  41. ^ “A104203”. The On-Line Encyclopedia of Integer Sequences.
  42. ^ Lomont, J.S.; Brillhart, John (2001). Elliptic Polynomials. CRC Press. pp. 12, 44. ISBN 1-58488-210-7.
  43. ^ a b c d “A193543 – Oeis”.
  44. ^ Lomont, J.S.; Brillhart, John (2001). Elliptic Polynomials. CRC Press. ISBN 1-58488-210-7. p. 79, eq. 5.36
  45. ^ Lomont, J.S.; Brillhart, John (2001). Elliptic Polynomials. CRC Press. ISBN 1-58488-210-7. p. 79, eq. 5. 36 and p. 78, eq. 5.33
  46. ^ Robinson (2019a)
  47. ^
  48. ^ Eymard & Lafon (2004) p. 234
  49. ^ Armitage, J. V.; Eberlein, W. F. (2006). Elliptic Functions. Cambridge University Press. p. 49. ISBN 978-0-521-78563-1.
  50. ^ The identity
  51. ^ a b “A289695 – Oeis”.
  52. ^ Wall, H. S. (1948). Analytic Theory of Continued Fractions. Chelsea Publishing Company. pp. 374–375.
  53. ^ Reinhardt & Walker (2010a) §22.20(ii)
  54. ^ Carlson (2010) §19.8
  55. ^ Dieckmann, Andreas. “Collection of Infinite Products and Series”.
  56. ^ Reinhardt & Walker (2010a) §22.12.12; Vigren & Dieckmann (2020) p. 7
  57. ^ In general,
  58. ^ Reinhardt & Walker (2010a) §22.11
  59. ^ Reinhardt & Walker (2010a) §22.2.E7
  60. ^ Berndt (1994) p. 247, 248, 253
  61. ^ Reinhardt & Walker (2010a) §22.11.E1
  62. ^ http://oeis.org/A175576
  63. ^ Berndt, Bruce C. (1989). Ramanujan’s Notebooks Part II. Springer. ISBN 978-1-4612-4530-8. p. 96
  64. ^ Levin (2006) p. 515
  65. ^ Levin (2006); Robinson (2019b)
  66. ^ a b Cox (2012) p. 508, 509
  67. ^ a b c d Arakawa, Tsuneo; Ibukiyama, Tomoyoshi; Kaneko, Masanobu (2014). Bernoulli Numbers and Zeta Functions. Springer. ISBN 978-4-431-54918-5. p. 203—206
  68. ^ Equivalently,
  69. ^ The Bernoulli numbers can be determined by an analogous recurrence:
  70. ^ Katz, Nicholas M. (1975). “The congruences of Clausen — von Staudt and Kummer for Bernoulli-Hurwitz numbers”. Mathematische Annalen. 216 (1): 1–4. See eq. (9)
  71. ^ Hurwitz, Adolf (1963). Mathematische Werke: Band II (in German). Springer Basel AG. p. 370
  72. ^ Arakawa et al. (2014) define
  73. ^ Peirce (1879). Guyou (1887) and Adams (1925) introduced transverse and oblique aspects of the same projection, respectively. Also see Lee (1976). These authors write their projection formulas in terms of Jacobi elliptic functions, with a square lattice.
  74. ^ Adams (1925)
  75. ^ Adams (1925); Lee (1976).
  76. ^ Rančić, Purser & Mesinger (1996); McGregor (2005).

External links[edit]

References[edit]