View source for Jacobi elliptic functions

In [[mathematics]], the '''Jacobi elliptic functions''' are a set of basic [[elliptic function]]s, and auxiliary [[theta function]]s, that are of historical importance.  They are found in the description of the motion of a [[pendulum]] (see also [[pendulum (mathematics)]]), as well as in the design of the electronic [[elliptic filter]]s. While [[trigonometry|trigonometric functions]] are defined with reference to a circle, the Jacobi elliptic functions are a generalization which refer to other [[conic section]]s, 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 {{harvs|txt|first=Carl Gustav Jakob |last=Jacobi|authorlink=Carl Gustav Jakob Jacobi|year=1829}}.

==Introduction==
[[Image:JacobiFunctionAbstract.png|width322px|thumb|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 the complex plane of the argument ''u'', the twelve functions form a repeating lattice of simple [[Zeros and poles|poles and zeroes]].<ref name="DLMF22">{{cite web|url=http://dlmf.nist.gov/22|title=NIST Digital Library of Mathematical Functions (Release 1.0.17)|editor-last=Olver|editor-first=F. W. J.|display-editors=et al |date=2017-12-22|publisher=National Institute of Standards and Technology|access-date=2018-02-26 }}</ref> 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 period]]s 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.

Mathematically, Jacobian elliptic functions are doubly periodic [[meromorphic]] functions on the [[complex plane]]. Since they are doubly periodic, 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 4''K'' and the second 4''K''&prime;, where ''K'' and ''K''&prime; are the [[quarter period]]s. Each function has two zeroes and two poles at opposite positions on the torus. Among the points {{nowrap|0, ''K'', ''K'' + ''iK''&prime;, ''iK''&prime;}} there is one zero and one pole.

The Jacobian elliptic functions are then the unique doubly periodic, [[meromorphic]] functions satisfying the following three properties:
* There is a simple zero at the corner p, and a simple pole at the corner&nbsp;q.
* The step from p to q is equal to half the period of the function pq&nbsp;''u''; that is, the function pq&nbsp;''u'' is periodic in the direction pq, with the period being twice the distance from p to q. The function pq&nbsp;''u'' is also periodic in the other two directions, with a period such that the distance from p to one of the other corners is a quarter period.
* If the function pq&nbsp;''u'' is expanded in terms of ''u'' at one of the corners, the leading term in the expansion has a coefficient of&nbsp;1. In other words, the leading term of the expansion of pq&nbsp;''u'' at the corner p is ''u''; the leading term of the expansion at the corner q is 1/''u'', and the leading term of an expansion at the other two corners is&nbsp;1.

 {{multiple image
 |align=center
 |footer=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.<ref>http://nbviewer.ipython.org/github/empet/Math/blob/master/DomainColoring.ipynb</ref> All have values of the ''k'' parameter equal to 0.8.
 | image1 = sn-k-08.png
 | alt1=Elliptic Jacobi function sn, ''k''&nbsp;=&nbsp;0.8
 | caption1=Jacobi elliptic function sn

 | image2 = cn-k-08.png
 | alt2=Elliptic Jacobi function cn, ''k''&nbsp;=&nbsp;0.8
 | caption2=Jacobi elliptic function cn

 | image3 = dn-k-08.png
 | alt3=Elliptic Jacobi function dn, ''k''&nbsp;=&nbsp;0.8
 | caption3=Jacobi elliptic function dn

 | image4 = sc-k-08.png
 | alt4=Elliptic Jacobi function sc, ''k''&nbsp;=&nbsp;0.8
 | caption4=Jacobi elliptic function sc
 }}
{{clear}}

==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 ''k''<sup>2</sup> = ''m'', or in terms of the '''[[modular angle]]''' &alpha;, where ''m''&nbsp;= &nbsp;sin<sup>2</sup>&nbsp;&alpha;. The complements of ''k'' and ''m'' are defined as ''m' ''=1-m and <math>k'=\sqrt{m'}</math>. 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 [[Christof Gudermann|Gudermann]] and [[James Whitbread Lee Glaisher|Glaisher]] and is not Jacobi's original notation.)

The parameter

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

:<math>\operatorname{pq}\,\,\,\,\operatorname{p'q'}= \operatorname{pq'}\,\,\,\operatorname{p'q}</math>

from which other commonly used relationships can be derived:

:<math>\operatorname{pr}/\operatorname{qr}=\operatorname{pq}</math>
:<math>\operatorname{pr}\,\,\operatorname{rq}=\operatorname{pq}</math>
:<math>1/\operatorname{qp}=\operatorname{pq}</math>

The multiplication rule follows immediately from the identification of the elliptic functions with the [[Neville theta function]]s<ref name="Neville1944">{{cite book |last=Neville |first=Eric Harold |date=1944 |title=Jacobian Elliptic Functions |url=https://archive.org/details/jacobianelliptic00neviuoft |location=Oxford |publisher=Oxford University Press |author-link=Eric Harold Neville}}</ref>
:

:<math>\operatorname{pq}(u,m)=\frac{\theta_p(u,m)}{\theta_q(u,m)}</math>

=={{anchor|sn|cn|dn|am}}Definition as inverses of elliptic integrals==

[[File:Modell der elliptischen Funktion φ=am (u, k) durch eine Fläche -Schilling V, 1 - 317-.jpg|thumb|Model of amplitude (measured along vertical axis) as a function of independent variables ''u'' and ''k'']]

The above definition, in terms of the unique meromorphic functions satisfying certain properties, is quite abstract.  There is a simpler, but completely equivalent definition, giving the elliptic functions as inverses of the incomplete [[elliptic integral]] of the first kind.   Let

:<math>u=\int_0^\varphi \frac{\mathrm d\theta} {\sqrt {1-m \sin^2 \theta}}\, . </math>

Then the ''elliptic sine'' sn&nbsp;''u'' (Latin: ''sinus amplitudinis'') is given by

:<math>\operatorname {sn} u = \sin \varphi\,</math>

and the ''elliptic cosine'' cn&nbsp;''u'' (Latin: ''cosinus amplitudinis'') is given by

:<math>\operatorname {cn} u = \cos \varphi</math>

and the ''delta amplitude'' dn&nbsp;''u'' (Latin: ''delta amplitudinis'')

:<math>\operatorname {dn} u = \sqrt {1-m\sin^2 \varphi}\, .</math>

Here, the angle <math>\varphi</math> is called the '''amplitude'''.  On occasion, dn&nbsp;''u''&nbsp;=&nbsp;&Delta;(''u'') is called the '''delta amplitude'''.  In the above, the value ''m'' is a free parameter, usually taken to be real, 0&nbsp;≤&nbsp;''m''&nbsp;≤&nbsp;1, and so the elliptic functions can be thought of as being given by two variables, the amplitude <math>\varphi</math> and the parameter&nbsp;''m''.

The remaining nine elliptic functions are easily built from the above three, and are given in a section below.

Note that when <math>\varphi=\pi/2</math>, that ''u'' then equals the [[quarter period]]&nbsp;''K''.

==Definition as trigonometry - the Jacobi ellipse==
[[File:Jacobi Elliptic Functions (on Jacobi Ellipse).svg|right|thumb|upright=1.5|Plot of the Jacobi ellipse (''x''<sup>2</sup>+''y''<sup>2</sup>/b<sup>2</sup>=1, ''b'' real) and the twelve Jacobi Elliptic functions ''pq(u,m)'' for particular values of angle &phi; and  parameter ''b''.  The solid curve is the ellipse, with  ''m''=1-1/b<sup>2</sup> and   ''u''=''F(&phi;,m)'' where ''F(.,.)'' is the [[elliptic integral]] of the first kind. 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.]]

<math> \cos \varphi, \sin \varphi </math> are defined on the unit circle, with radius ''r''&nbsp;=&nbsp;1 and angle <math>\varphi =</math> arc length of the unit circle measured from the positive ''x''-axis. Similarly, Jacobi elliptic functions are defined on the unit ellipse{{citation needed|reason= see talk of this section.  |date=July 2016}}, with ''a''&nbsp;=&nbsp;1.  Let

:<math>
\begin{align}
& x^2 + \frac{y^2}{b^2} = 1, \quad  b >  1, \\
& m = 1 - \frac{1}{b^2}, \quad 0 < m < 1, \\
& x = r \cos \varphi, \quad y = r \sin \varphi
\end{align}
</math>

then:

:<math> r( \varphi,m) =  \frac{1} {\sqrt {1-m \sin^2 \varphi}}\, . </math>

For each angle <math>\varphi</math> the angular component of the arc length, the
'angular arc length' is computed. An advantage of angular arc length is that we can calculate a total arc length
for hyperbolas. For the ellipse the angular arc length is:
:<math>u = u(\varphi,m)=\int_0^\varphi r(\theta,m) \, d\theta .</math>
Let <math>P=(x,y)</math> be the point on the ellipse with angular arc length <math>u(\varphi,m)</math> and let <math>P'=(x',y')</math> be the point on the unit circle with the angular arc length <math>\varphi</math>, (note that origin <math>O, P'</math> and <math>P</math> are on the same straight line). The familiar relations from the unit circle (<math>OP' = 1</math>):  
:<math> x' = \cos \varphi, \quad y' = \sin \varphi</math>
imply for the ellipse:
:<math>x' = \operatorname{cn}(u,m),\quad y' = \operatorname{sn}(u,m).</math>
So the projections of the intersection point <math>P'</math> of the line <math>OP</math> with the unit circle on the ''x''- and ''y''-axes are simply <math>\operatorname{cn}(u,m)</math> and <math>\operatorname{sn}(u,m)</math>. These projections may be interpreted as 'definition as trigonometry'. In short:

:<math> \operatorname{cn}(u,m) = \frac{x(u,m)}{r(u,m)}, \quad \operatorname{sn}(u,m) = \frac{y(u,m)}{r(u,m)}, \quad \operatorname{dn}(u,m) = \frac{1}{r(u,m)}. </math>

For the <math>x</math> and <math>y</math> value of the point <math>P</math> with 
<math>u</math> and parameter <math>m</math> we get, after inserting the relation:
:<math>r(\varphi,m) = \frac 1 {\operatorname{dn}(u,m)} </math>
into:<math>x = r(\varphi,m) \cos (\varphi), y = r(\varphi,m) \sin (\varphi)</math> that:
:<math> x = \frac{\operatorname{cn}(u,m)} {\operatorname{dn}(u,m)},\quad y = \frac{\operatorname{sn}(u,m)} {\operatorname{dn}(u,m)}.</math>
The latter relations for the ''x''- and ''y''-coordinates of points on the unit ellipse may be considered as generalization of the relations <math> x = \cos \varphi, y = \sin \varphi</math> 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 (''&phi;'',dn) with <math>r=\sqrt{x^2+y^2}</math>
{| class="wikitable" style="text-align:center;
|+ Jacobi Elliptic Functions pq[u,m] as functions of {x,y,r} and {&phi;,dn}
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! c
|1 || <math>x/y=\cot(\phi)</math> || <math>x/r=\cos(\phi)</math> ||  <math>x=\cos(\phi)/dn</math>
|-
! s
|<math>y/x=\tan(\phi)</math> || 1 ||<math>y/r=\sin(\phi)</math> || <math>y=\sin(\phi)/dn</math>
|-
! n
|<math>r/x=\sec(\phi)</math> || <math>r/y=\csc(\phi)</math> || 1 || <math>r=1/dn</math>
|-
! d
| <math>1/x=dn \sec(\phi)</math> || <math>1/y=dn \csc(\phi)</math> ||  <math>1/r=dn </math> || 1
|}

==Definition in terms of Jacobi theta functions==

Equivalently, Jacobi's elliptic functions can be defined in terms of his [[theta function]]s. If we abbreviate <math>\vartheta(0;\tau)</math> as <math>\vartheta</math>, and <math>\vartheta_{01}(0;\tau), \vartheta_{10}(0;\tau), \vartheta_{11}(0;\tau)</math> respectively as <math>\vartheta_{01}, \vartheta_{10}, \vartheta_{11}</math> (the ''theta constants'') then the [[elliptic modulus]] ''k'' is <math>k=\left({\vartheta_{10} \over \vartheta}\right)^2</math>. If we set <math>u = \pi \vartheta^2 z</math>, we have

: <math id="def.theta">
\begin{align}
\operatorname{sn}(u; k) & = -{\vartheta \vartheta_{11}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)} \\[7pt]
\operatorname{cn}(u; k) & = {\vartheta_{01} \vartheta_{10}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)} \\[7pt]
\operatorname{dn}(u; k) & = {\vartheta_{01} \vartheta(z;\tau) \over \vartheta \vartheta_{01}(z;\tau)}
\end{align}
</math>

Since the Jacobi functions are defined in terms of the elliptic modulus <math>k(\tau)</math>, we need to invert this and find <math>\tau</math> in terms of <math>k</math>. We start from <math>k' = \sqrt{1-k^2}</math>, the ''complementary modulus''. As a function of <math>\tau</math> it is

:<math>k'(\tau) = \left({\vartheta_{01} \over \vartheta}\right)^2.</math>

Let us first define

:<math>\ell = {1 \over 2} {1-\sqrt{k'} \over 1+\sqrt{k'}} =
{1 \over 2} {\vartheta - \vartheta_{01} \over \vartheta + \vartheta_{01}}.</math>

Then define the [[nome (mathematics)|nome]] as <math>q = \exp (\pi i \tau)</math> and expand <math>\ell</math> as a [[power series]] in the nome <math>q</math>, we obtain

:<math>\ell = {q + q^9 + q^{25} + \cdots \over 1 + 2q^4 + 2q^{16} + \cdots}.</math>

[[Reversion of series]] now gives

:<math>q = \ell + 2\ell^5 + 15\ell^9 + 150\ell^{13} + 1707\ell^{17} + 20910\ell^{21} + 268616\ell^{25} + \cdots.</math>

Since we may reduce to the case where the imaginary part of <math>\tau</math> is greater than or equal to <math>\sqrt{3}/2</math>, we can assume the absolute value of <math>q</math> is less than or equal to <math>\exp(-\pi\sqrt{3}/2) \approx 0.0658 </math>; for values this small the above series converges very rapidly and easily allows us to find the appropriate value for <math>q</math>.

==Definition in terms of Neville theta functions==

The Jacobi Elliptic functions can be defined very simply using the [[Neville theta functions]]:<ref name="WolframJE"/>

:<math>pq(u,m)=\frac{\theta_p(u,m)}{\theta_q(u,m)}</math>

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

== Jacobi Transformations==

===The Jacobi imaginary transformations===

[[File:JacobiElliptic.HT.svg|right|thumb|upright=1.5|Plot of the degenerate Jacobi curve (x<sup>2</sup>+y<sup>2</sup>/b<sup>2</sup>=1, b=infinity) and the twelve Jacobi Elliptic functions pq(u,1) for a particular value of angle &phi;.  The solid curve is the degenerate ellipse (x<sup>2</sup>=1) with m=1 and u=F(&phi;,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:<ref name="W&W">{{cite book |last1=Whittaker |first1=E.T. |last2=Watson |first2=G.N.|date=1940 |title=A Course in Modern Analysis |url=https://archive.org/details/courseofmodernan00whit |location=New York, USA |publisher=The MacMillan Co.|isbn=978-0-521-58807-2|author-link=A Course of Modern Analysis}}</ref>{{rp|506}}

:<math>\operatorname{cn}(u, m)=    \operatorname{nc}(i\,u,1\!-\!m)</math>
:<math>\operatorname{sn}(u, m)= -i \operatorname{sc}(i\,u,1\!-\!m)</math>
:<math>\operatorname{dn}(u, m)=    \operatorname{dc}(i\,u,1\!-\!m)</math>

Using the multiplication rule, all other functions may be expressed in terms of the above three. The transformations may be generally written as <math>pq(u,m)=\gamma_{pq} pq'(i\,u,1\!-\!m)</math>. The following table gives the <math>\gamma_{pq} pq'(i\,u,1\!-\!m)</math> for the specified pq(''u,m'').<ref name="WolframJE">{{cite web |url=http://functions.wolfram.com/EllipticFunctions/JacobiAmplitude/introductions/JacobiPQs/ShowAll.html |title=Introduction to the Jacobi elliptic functions |date=2018 |website=The Wolfram Functions Site |publisher=Wolfram Research, Inc. |access-date=January 7, 2018}}</ref> (The arguments <math>(i\,u,1\!-\!m)</math> are suppressed)
:{| class="wikitable" style="text-align:center;
|+ Jacobi Imaginary transformations <math>\gamma_{pq}\operatorname{pq}'(i\,u,1\!-\!m)</math>
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! 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 functions|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.<ref name="Neville1944"/>{{rp|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<ref name="Neville1944"/>{{rp|308}} yield expressions for the elliptic functions in terms with alternate values of ''m''. The transformations may be generally written as <math>pq(u,m)=\gamma_{pq} pq'(k\,u,1/m)</math>. The following table gives the <math>\gamma_{pq} pq'(k\,u,1/m)</math> for the specified pq(''u,m'').<ref name="WolframJE"/> (The arguments <math>(k\,u,1/m)</math> are suppressed)

:{| class="wikitable" style="text-align:center;
|+ Jacobi Real transformations <math>\gamma_{pq}\operatorname{pq}'(k\,u,1/m)</math>
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! c
|1||<math>k</math> ds || dn || dc
|-
! s
|<math>\frac{1}{k}</math> sd|| 1 || <math>\frac{1}{k}</math> sn || <math>\frac{1}{k}</math> sc
|-
! n
| nd || <math>k</math> ns || 1 || nc
|-
! d
| cd || <math>k</math> cs || cn || 1
|}

=== Other Jacobi transformations ===

Jacobi's real and imaginary transformations can be combined in various ways to yield three more simple transformations
.<ref name="Neville1944"/>{{rp|214}} The real and imaginary transformations are two transformations in a group ([[Dihedral group of order 6|D<sub>3</sub>]] or [[Anharmonic group]]) of six transformations. If

:<math>\mu_R(m) = 1/m</math>

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

:<math>\mu_I(m) = 1-m = m'</math>

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:

:<math>
\begin{align}
\mu_{IR}(m)&=&\mu_I(\mu_R(m))&=&-m'/m \\
\mu_{RI}(m)&=&\mu_R(\mu_I(m))&=&1/m' \\
\mu_{RIR}(m)&=&\mu_R(\mu_I(\mu_R(m)))&=&-m/m'
\end{align}
</math>

These five transformations, along with the identity transformation (&mu;<sub>U</sub>(m)=m) yield the 6 element group. With regard to the Jacobi elliptic functions, the general transformation can be expressed using just three functions:

:<math>\operatorname{cs}(u,m)=\gamma_i \operatorname{cs'}(\gamma_i u, \mu_i(m))</math>
:<math>\operatorname{ns}(u,m)=\gamma_i \operatorname{ns'}(\gamma_i u, \mu_i(m))</math>
:<math>\operatorname{ds}(u,m)=\gamma_i \operatorname{ds'}(\gamma_i u, \mu_i(m))</math>

where ''i'' = U, I, IR, R, RI, or RIR, identifying the transformation, &gamma;<sub>i</sub> 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.<ref name="Neville1944"/>{{rp|214}} (As usual, k<sup>2</sup>=m, 1-k<sup>2</sup>=k<sub>1</sub><sup>2</sup>=m' and the arguments (<math>\gamma_i u, \mu_i(m)</math>) are suppressed)

:{| class="wikitable" style="text-align:center;
|+ Parameters for the six transformations
!Transformation i||<math>\gamma_i</math>||<math>\mu_i(m)</math>||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 k<sub>1</sub>|| 1/m'  || ns || ds || cs
|-
! RIR
| k<sub>1</sub>  || -m/m' || cs || ds || ns
|-
|}

Thus, for example, we may build the following table for the RIR transformation.<ref name="WolframJE"/> The transformation is generally written <math>\operatorname{pq}(u,m)=\gamma_{pq}\,\operatorname{pq'}(k'\,u,-m/m')</math> (The arguments <math>(k'\,u,-m/m')</math> are suppressed)

:{| class="wikitable" style="text-align:center;
|+ The RIR transformation <math>\gamma_{pq}\,\operatorname{pq'}(k'\,u,-m/m')</math>
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! c
|1|| k' cs || cd || cn
|-
! s
|<math>\frac{1}{k'}</math> sc|| 1 ||<math>\frac{1}{k'}</math> sd ||<math> \frac{1}{k'}</math> sn
|-
! n
| dc || <math>k'</math> ds || 1 || dn
|-
! d
| nc || <math>k'</math> 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.<ref name="Neville1944"/>{{rp|p.215}}

==The Jacobi hyperbola==
[[File:Jacobi Elliptic Functions (on Jacobi Hyperbola).svg|right|thumb|upright=1.5|Plot of the Jacobi hyperbola (''x''<sup>2</sup>+''y''<sup>2</sup>/b<sup>2</sup>=1, ''b'' imaginary) and the twelve Jacobi Elliptic functions ''pq(u,m)'' for particular values of angle &phi; and  parameter ''b''.  The solid curve is the hyperbola, with  ''m''=1-1/b<sup>2</sup> and   ''u''=''F(&phi;,m)'' where ''F(.,.)'' is the [[elliptic integral]] of the first kind. The dotted curve is the unit circle. For the ds-dc triangle,''&sigma;''= ''Sin(&phi;)Cos(&phi;)''.]]

Introducing complex numbers, our ellipse has an associated hyperbola:

:<math> xx^2 - \frac{yy^2}{b^2} = 1 </math> 
from applying Jacobi's imaginary transformation<ref name="WolframJE"/> to the elliptic functions in the above equation for ''x'' and&nbsp;''y''.

:<math> xx = \frac{1} {\operatorname{dn}(u,1-m)},\quad yy = \frac{ \operatorname{sn}(u,1-m)} {\operatorname{dn}(u,1-m)}</math>

It follows that we can put <math> x=\operatorname{dn}(u,1-m), y=\operatorname{sn}(u,1-m)</math>. So our ellipse has a dual ellipse with m replaced by 1-m. This leads to the complex torus mentioned in the  Introduction.<ref>https://paramanands.blogspot.co.uk/2011/01/elliptic-functions-complex-variables.html#.WlHhTbp2t9A</ref> 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.

=={{anchor|ns|nc|nd|sc|sd|dc|ds|cs|cd|pg}}Minor functions==

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

:<math>
\begin{align}
\operatorname{ns}(u) = \frac{1}{\operatorname{sn}(u)}, \qquad \operatorname{nc}(u)  = \frac{1}{\operatorname{cn}(u)},  \qquad 
\operatorname{nd}(u)  = \frac{1}{\operatorname{dn}(u)}.
\end{align}
</math>

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

:<math>
\begin{align}
\operatorname{sc}(u)  = \frac{\operatorname{sn}(u)}{\operatorname{cn}(u)},  \qquad 
\operatorname{sd}(u)  = \frac{\operatorname{sn}(u)}{\operatorname{dn}(u)},  \qquad \operatorname{dc}(u)  = \frac{\operatorname{dn}(u)}{\operatorname{cn}(u)},  \qquad \operatorname{ds}(u)  = \frac{\operatorname{dn}(u)}{\operatorname{sn}(u)},  \qquad \operatorname{cs}(u)  = \frac{\operatorname{cn}(u)}{\operatorname{sn}(u)},  \qquad 
\operatorname{cd}(u) = \frac{\operatorname{cn}(u)}{\operatorname{dn}(u)}.
\end{align}
</math>
More compactly, we have

:<math>\operatorname{pq}(u)=\frac{\operatorname{pn}(u)}{\operatorname{qn}(u)}</math>

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

== Periodicity, poles, and residues ==

[[File:JacobiEllipticFunctions.svg|thumb|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 amplitude below 1/3 are colored black, roughly indicating the location of a zero, while regions with amplitude 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 amplitude, differing only in sign. Each function pq(u,m) has an inverse function 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 double periodicity of the Jacobi elliptic functions may be expressed as:

:<math>\operatorname{pq}(u + 2 \alpha K(m) + 2 i \beta K(1-m)\,,\,m)=(-1)^\gamma \operatorname{pq}(u,m)</math>

where &alpha; and &beta; 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 (&gamma;) is given in the following table:

:{| class="wikitable" style="text-align:center;
|+ <math>\gamma</math>
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! c
|0||&beta; || &alpha;+&beta; || &alpha;
|-
! s
|&beta; || 0 || &alpha; || &alpha;+&beta;
|-
! n
| &alpha;+&beta; || &alpha; || 0 || &beta;
|-
! d
| &alpha; || &alpha;+&beta; || &beta; || 0
|}

When the factor (-1)<sup>&gamma;</sup> 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 &alpha; when &alpha; 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 &beta; have full periods of 2K(m) and 4iK(1-m), while those with &alpha; + &beta; 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 amplitude 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:

:{| class="wikitable" style="text-align:center; width:200px”;
|+ Residues of Jacobi Elliptic Functions
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! width="40pt"|c
! width="40pt"|s
! width="40pt"|n
! width="40pt"|d
|-
!rowspan="6"|p
|-
! height="40pt" |c
|  ||1||<math>-\frac{i}{k}</math>||<math>-\frac{1}{k}</math>
|-
! height="40pt" |s
| <math>-\frac{1}{k'}</math>|| ||<math>\frac{1}{k}</math>||<math>-\frac{i}{k\,k'}</math>
|-
! height="40pt" |n
|<math>-\frac{1}{k'}</math>||1|| ||<math>-\frac{i}{k'}</math>
|-
! height="40pt" | d
| -1 || 1 || <math>-i</math> ||  
|-
|}

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.

==Relations between squares of the functions==

Relations between squares of the functions can be derived from two basic relationships (Arguments (''u'',''m'') suppressed):

:<math>\operatorname{cn}^2+\operatorname{sn}^2=1</math>
:<math>\operatorname{cn}^2+m' \operatorname{sn}^2=\operatorname{dn}^2</math>

where ''m&nbsp;+&nbsp;m'&nbsp;''=&nbsp;1 and ''m''&nbsp;=&nbsp;''k''<sup>2</sup>. Multiplying by any function of the form ''nq'' yields more general equations:

:<math>\operatorname{cq}^2+\operatorname{sq}^2=\operatorname{nq}^2</math>
:<math>\operatorname{cq}^2+m' \operatorname{sq}^2=\operatorname{dq}^2</math>

With ''q''=''d'', these correspond trigonometrically to the equations for the unit circle (<math>x^2+y^2=r^2</math>) and the unit ellipse  (<math>x^2+m' y^2=1</math>), with ''x=cd'', ''y=sd'' and ''r=nd''. Using the multiplication rule, other relationships may be derived. For example:

:<math>
-\operatorname{dn}^2+m'= -m\operatorname{cn}^2 = m\operatorname{sn}^2-m
</math>

:<math>
-m'\operatorname{nd}^2+m'= -mm'\operatorname{sd}^2 = m\operatorname{cd}^2-m
</math>

:<math>
m'\operatorname{sc}^2+m'= m'\operatorname{nc}^2 = \operatorname{dc}^2-m
</math>

:<math>
\operatorname{cs}^2+m'=\operatorname{ds}^2=\operatorname{ns}^2-m
</math>

==Addition theorems==

The functions satisfy the two square relations

:<math>\operatorname{cn}^2(u,k) + \operatorname{sn}^2(u,k) = 1,\,</math>

:<math>\operatorname{dn}^2(u,k) + k^2 \operatorname{sn}^2(u,k) = 1.\,</math>

From this we see that (cn, sn, dn) parametrizes an [[elliptic curve]] which is the intersection of the two [[quadric]]s 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<ref name="DLMF22"/>

:<math>
\begin{align}
\operatorname{cn}(x+y) & =
{\operatorname{cn}(x) \operatorname{cn}(y)
- \operatorname{sn}(x) \operatorname{sn}(y) \operatorname{dn}(x) \operatorname{dn}(y)
\over {1 - k^2 \operatorname{sn}^2 (x) \operatorname{sn}^2 (y)}}, \\[8pt]
\operatorname{sn}(x+y) & =
{\operatorname{sn}(x) \operatorname{cn}(y) \operatorname{dn}(y) +
\operatorname{sn}(y) \operatorname{cn}(x) \operatorname{dn}(x)
\over {1 - k^2 \operatorname{sn}^2 (x) \operatorname{sn}^2 (y)}}, \\[8pt]
\operatorname{dn}(x+y) & =
{\operatorname{dn}(x) \operatorname{dn}(y) - k^2 \operatorname{sn}(x) \operatorname{sn}(y) \operatorname{cn}(x) \operatorname{cn}(y)
\over {1 - k^2 \operatorname{sn}^2 (x) \operatorname{sn}^2 (y)}}.
\end{align}
</math>

Double angle formulae can be easily derived from the above equations by setting ''x''=''y''.<ref name="DLMF22"/> Half angle formulae<ref name="WolframJE"/><ref name="DLMF22"/> are all of the form:

:<math>\operatorname{pq}(\tfrac{1}{2}u,m)^2 = f_p/f_q</math>

where:
:<math>f_c = \operatorname{cn}(u,m)+\operatorname{dn}(u,m)</math>
:<math>f_s = 1-\operatorname{cn}(u,m)</math>
:<math>f_n = 1+\operatorname{dn}(u,m)</math>
:<math>f_d = (1+\operatorname{dn}(u,m))-m(1-\operatorname{cn}(u,m))</math>

==Expansion in terms of the nome==

Let the [[nome (mathematics)|nome]] be <math>q=\exp(-\pi K'/K)</math> and let the argument be <math>v=\pi u /(2K)</math>. Then the functions have expansions as [[Lambert series]]

:<math>\operatorname{sn}(u)=\frac{2\pi}{K\sqrt{m}}
\sum_{n=0}^\infty \frac{q^{n+1/2}}{1-q^{2n+1}} \sin ((2n+1)v),</math>

:<math>\operatorname{cn}(u)=\frac{2\pi}{K\sqrt{m}}
\sum_{n=0}^\infty \frac{q^{n+1/2}}{1+q^{2n+1}} \cos ((2n+1)v),</math>

:<math>\operatorname{dn}(u)=\frac{\pi}{2K} + \frac{2\pi}{K}
\sum_{n=1}^\infty \frac{q^{n}}{1+q^{2n}} \cos (2nv).</math>

==Jacobi elliptic functions as solutions of nonlinear ordinary differential equations==
The [[derivative]]s of the three basic Jacobi elliptic functions are:

:<math>
\frac{\mathrm{d}}{\mathrm{d}z} \operatorname{sn}(z) = \operatorname{cn}(z) \operatorname{dn}(z),</math>

<!-- extra blank line between two lines of "displayed" [[TeX]] for legibility -->

:<math>\frac{\mathrm{d}}{\mathrm{d}z} \operatorname{cn}(z) = -\operatorname{sn}(z) \operatorname{dn}(z),</math>

<!-- extra blank line between two lines of "displayed" [[TeX]] for legibility -->

:<math>\frac{\mathrm{d}}{\mathrm{d}z} \operatorname{dn}(z) = - k^2 \operatorname{sn}(z) \operatorname{cn}(z).
</math>

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

{| class="wikitable" style="text-align:center;
|+ Derivatives <math>\frac{d}{du} \operatorname{pq}(u,m)</math>
!colspan="2" rowspan="2"|
!colspan="4"|q
|-
! c
! s
! n
! d
|-
!rowspan="6"|p
|-
! 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|addition theorems above]] and for a given ''k'' with 0&nbsp;<&nbsp;''k''&nbsp;<&nbsp;1 the major functions are  therefore solutions to the following nonlinear [[ordinary differential equation]]s:

* <math>\operatorname{sn}(x)</math> solves the differential equations
:: <math>\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1+k^2) y - 2 k^2 y^3 = 0</math>
: and
:: <math> \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 y^2)</math>

* <math>\operatorname{cn}(x)</math> solves the differential equations
:: <math>\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1-2k^2) y + 2 k^2 y^3 = 0</math>
: and
:: <math> \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 + k^2 y^2)</math>

* <math>\operatorname{dn}(x)</math> solves the differential equations
:: <math>\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} - (2 - k^2) y + 2 y^3 = 0</math>
: and
:: <math> \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (y^2 - 1) (1 - k^2 - y^2)</math>

==Approximation in terms of hyperbolic functions==

The Jacobi elliptic functions can be expanded in terms of the hyperbolic functions.  When <math>m</math> is close to unity, such that <math>m'^2</math> and higher powers of  <math>m'</math> can be neglected, we have:
* sn(''u''):
::<math>\operatorname{sn} (u,m)\approx \tanh (u)+\frac{1}{4}m'(\sinh (u)\cosh (u) -u)\operatorname{sech}^2 (u).</math>
* cn(''u''):
::<math>\operatorname{cn} (u,m)\approx \operatorname{sech} (u)-\frac{1}{4} m'(\sinh (u)\cosh (u) -u)\tanh u \operatorname{sech} (u).</math>
* dn(''u''):
::<math>\operatorname{dn} (u,m) \approx \operatorname{sech} (u)+\frac{1}{4} m'(\sinh (u)\cosh(u) +u)\tanh (u) \operatorname{sech} (u) .</math>
* am(''u''):
::<math>\operatorname{am} (u,m) \approx \operatorname{gd} (u)+\frac{1}{4}m'(\sinh (u)\cosh (u) -u)\operatorname{sech} (u) .</math>

== {{anchor|arcsn|arccn|arcdn}}Inverse functions ==

The inverses of the Jacobi elliptic functions can be defined similarly to the [[inverse trigonometric functions]]; if <math>x=\operatorname{sn}(\xi, k)</math>, <math>\xi=\mathrm{arcsn}(x, k)</math>. They can be represented as elliptic integrals,<ref>{{dlmf|title=§22.15 Inverse Functions|first1=W. P.|last1=Reinhardt|first2=P. L.|last2=Walker|id=22.15}}</ref><ref>{{cite web|last=Ehrhardt|first=Wolfgang|title=The AMath and DAMath Special Functions: Reference Manual and Implementation Notes|url=http://www.wolfgang-ehrhardt.de/specialfunctions.pdf|accessdate=17 July 2013|page=42}}</ref><ref>{{cite book|last1=Byrd|first1=P.F.|last2=Friedman|first2=M.D.|title=Handbook of Elliptic Integrals for Engineers and Scientists|date=1971|publisher=Springer-Verlag|location=Berlin|edition=2nd}}</ref> and power series representations have been found.<ref>{{cite journal|last=Carlson|first=B. C.|title=Power series for inverse Jacobian elliptic functions|journal=Mathematics of Computation|year=2008|volume=77|issue=263|pages=1615–1621|url=http://www.ams.org/journals/mcom/2008-77-263/S0025-5718-07-02049-2/S0025-5718-07-02049-2.pdf|accessdate=17 July 2013|doi=10.1090/s0025-5718-07-02049-2}}</ref><ref name="DLMF22"/>

*<math>\mathrm{arcsn}\,(x,k) = \int_0^x \frac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2t^2)}}</math>
*<math>\mathrm{arccn}\,(x,k) =\int_x^1 \frac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2+k^2t^2)}}</math>
*<math>\mathrm{arcdn}\,(x,k) = \int_x^1 \frac{\mathrm{d}t}{\sqrt{(1-t^2)(t^2+k^2-1)}}</math>

== 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]]
* [[theta function|Jacobi theta function]]
* [[Ramanujan theta function]]
* [[Dixon elliptic functions]]
* [[Abel elliptic functions]]
* [[Weierstrass's elliptic functions]]

==Notes==
<references/>

==References==
* {{AS ref|16|569}}
* [[Naum Akhiezer|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}}
* [[Alfred Cardew Dixon|A. C. Dixon]] [https://archive.org/details/117736039  The elementary properties of the elliptic functions, with examples] (Macmillan, 1894)
* [[Alfred George Greenhill]] [https://archive.org/details/applicationselli00greerich The applications of elliptic functions] (London, New York, Macmillan, 1892)
* H. Hancock [https://archive.org/details/lecturestheorell00hancrich Lectures on the theory of elliptic functions] (New York, J. Wiley & sons, 1910)
*{{Citation | last1=Jacobi | first1=C. G. J. | title=Fundamenta nova theoriae functionum ellipticarum | url=https://archive.org/details/fundamentanovat00jacogoog | publisher=Königsberg | language=Latin | isbn=978-1-108-05200-9 | id=Reprinted by Cambridge University Press 2012 | year=1829}}
*{{dlmf|first=William P. |last=Reinhardt|first2=Peter L. |last2=Walker|id=22|title=Jacobian Elliptic Functions}}
* {{fr icon}} [[P. Appell]] and E. Lacour [https://archive.org/details/principestheorie00apperich Principes de la théorie des fonctions elliptiques et applications] (Paris, Gauthier Villars, 1897)
* {{fr icon}} G. H. Halphen [https://archive.org/details/traitedesfonctio01halprich Traité des fonctions elliptiques et de leurs applications (vol. 1)] (Paris, Gauthier-Villars, 1886–1891)
* {{fr icon}} G. H. Halphen [https://archive.org/details/traitedesfonctio02halprich Traité des fonctions elliptiques et de leurs applications (vol. 2)] (Paris, Gauthier-Villars, 1886–1891)
* {{fr icon}} G. H. Halphen [https://archive.org/details/traitedesfonctio03halprich Traité des fonctions elliptiques et de leurs applications (vol. 3)] (Paris, Gauthier-Villars, 1886–1891)
* {{fr icon}} J. Tannery and J. Molk [http://gallica.bnf.fr/notice?N=FRBNF37258233 	Eléments de la théorie des fonctions elliptiques. Tome I, Introduction. Calcul différentiel. Ire partie] (Paris : Gauthier-Villars et fils, 1893)
* {{fr icon}} J. Tannery and J. Molk  [http://gallica.bnf.fr/notice?N=FRBNF37258241 	Eléments de la théorie des fonctions elliptiques. Tome II, Calcul différentiel. IIe partie] (Paris : Gauthier-Villars et fils, 1893)
* {{fr icon}} J. Tannery and J. Molk  [http://gallica.bnf.fr/notice?N=FRBNF37258245 	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)
* {{fr icon}} J. Tannery and J. Molk  [http://gallica.bnf.fr/notice?N=FRBNF37258246 	Eléments de la théorie des fonctions elliptiques. Tome IV, Calcul intégral. IIe partie, Applications] (Paris : Gauthier-Villars et fils, 1893)
* {{fr icon}} C. Briot and J. C. Bouquet [http://gallica.bnf.fr/notice?N=FRBNF30162167 Théorie des fonctions elliptiques] (  Paris : Gauthier-Villars, 1875)

==External links==
* {{springer|title=Jacobi elliptic functions|id=p/j054050}}
* {{MathWorld|JacobiEllipticFunctions|Jacobi Elliptic Functions}}
* {{youtube|DCXItCajCyo|title=Elliptic Functions and Elliptic Integrals}}, lecture by William A. Schwalm (4 hours)

{{DEFAULTSORT:Jacobi Elliptic Functions}}
[[Category:Jacobi elliptic functions]]