Tetration: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Dmitrii Kouznetsov
(→‎See also: similar functions)
mNo edit summary
 
(189 intermediate revisions by 8 users not shown)
Line 1: Line 1:
[[Image:TetrationReal.jpg|right|400px|thumb|Fig.1. Tetration  
{{subpages}}
<!--
{{Image|TetrationReal.jpg|right|400px|Fig.0. Tetration  
<math>\mathrm{tet}_b(x)</math> for  
<math>\mathrm{tet}_b(x)</math> for  
<math>b=\mathrm{e}</math>,
<math>b\!=\!\mathrm{e}</math>,
<math>b=2</math>,
<math>b\!=\!2</math>,
<math>b=\exp(1/\mathrm{e})</math>, and
<math>b\!=\!\exp(1/\mathrm{e})</math>, and
<math>b=\sqrt{2}</math> versus  
<math>b\!=\!\sqrt{2}~</math> versus  
<math>x</math>.]]
<math>x</math>.}}
<b>
old version !-->
Warning: This article is recently written; it may have poor wording and misprints.
'''Tetration''' is a rapidly growing [[mathematical function]], which was introduced in the 20th century and proposed for the representation of huge numbers in the [[mathematics of computation]]. For positive integer values of its argument <math>x</math>, tetration <math>\mathrm{tet}_b(x)</math> on base <math>b</math> can be defined with:
Consider to use article from wikipedia http://en.wikipedia.org/wiki/Tetration instead.
 
</b>
<math>{\ \mathrm{tet}_b(x) = \ \atop {\ }} {{\underbrace{b^{b^{\cdot^{\cdot^{b}}}}}} \atop {x}}~</math>
 
{{Image|Tetreal10bx10d.png|right|500px|Fig.1. Tetration
<math>\mathrm{tet}_b(x)</math> for  various <math>b\!>\!1</math> versus
<math>x</math>.}}
For real values of the argument <math>x</math> and various values of the base <math> b </math>, this <math>\mathrm{tet}_b(x)</math> is plotted in Fig.0.
 
{{Image|Tetration xb80.png|right|700px|Fig.1. Tetration
<math>f=\mathrm{tet}_b(x)</math> in the <math>x,b</math> plane with levels <math>g\!=\!\mathrm{const}</math>.}}
The map of <math> f=\mathrm{tet}_b(x)</math> in the <math>x,b</math> plane is shown in Fig.1 with levels  <math>f=\mathrm{const}</math>.
 
Up to year 2011, this function has not been listed among elementary functions, it is not implemented in [[programming languages]] and it is not used for the internal representation of data in computers.


'''Tetration''' is fastly growing [[mathematical function]], which was introduced in XX century and suggested for representation of huge numbers in [[mathematics of computation]]. However, up to year 2008, this function is not considered as elementary function, it is not implemented in [[programming languages]] and it is not used for the internal representation of data, at least in the commercial software.  
In this article, the generalizaiton of tetration for complex (and, in particular, real) values of its argument is described.
At base <math> b\!>\!1</math>, tetration is assumed to be a [[holomorphic function]], at least for positive values of the [[real part]] of its argument.
This tetration is used to construct the holomorphic extension of the iterated [[Exponential function|exponential]] <math>\exp^c(z)</math> for the case of non-integer values of the number <math>c</math> of iterations.


==Definiton==
==Definition==
For real <math>b>1</math>, Tetration <math>F=\mathrm{tet}_b</math> on the base <math>b</math> is function of complex variable, which is [[holomorphic]] at least in the range
For real <math>b>1</math>, Tetration <math>F=\mathrm{tet}_b</math> on the base <math>b</math> is a function of a complex variable, which is [[holomorphic]] at least in the range
<math> \{ z \in \mathbb{C} :~ \Re(z) > -2 \}</math>, bounded in the range  
<math> \{ z \in \mathbb{C} :~ \Re(z) > -2 \}</math>, bounded in the range  
<math> \{ z \in \mathbb{C} :~ |\Re(z)| \le 1 \}</math>, and satisfies conditions
<math> \{ z \in \mathbb{C} :~ |\Re(z)| \le 1 \}</math>, and satisfies conditions
: <math>  F(z+1) = \exp_b\! \big( F(z) \big) </math>  
: <math>  F(z\!+\!1) = \exp_b\! \big( F(z) \big) </math>  
: <math>  F(0) = 1 </math>  
: <math>  F(0) = 1 </math>  
: <math>  F\!\big(z^*\big) = F\big(z\big)^*</math>  
: <math>  F\!\big(z^*\big) = F\big(z\big)^*</math>  
at least within range <math> \Re(z)>-2 </math>.
at least within the range <math> ~ \Re(z)>-2 ~</math> .


According to this definition, tetration is [[superfunction]] of the [[exponential]]. This justifies the alternative name "superexponential" for this function and "superlogarithm" for the inverse function.
The definition above generalizes the definitions, recently suggested for the specific cases of base <math>b=\mathrm{e}</math>
<ref name="k">
{{cite journal
|author=D.Kouznetsov.
|title=Solutions of <math>F(z+1)=\exp(F(z))</math> in the complex <math>z</math>plane.
|journal=[[Mathematics of Computation]],
|year=2009
|volume=78
|pages=1647-1670
|url= http://www.ams.org/mcom/2009-78-267/S0025-5718-09-02188-7/home.html
|preprint: http://www.ils.uec.ac.jp/~dima/PAPERS/analuxp99.pdf
|doi=10.1090/S0025-5718-09-02188-7
}}</ref> and <math>b=2</math>
<ref name="k2">D.Kouznetsov. Ackermann functions of complex argument. Preprint of the Institute for Laser Science, UEC, 2008.
http://www.ils.uec.ac.jp/~dima/PAPERS/2008ackermann.pdf</ref>.


==Etymology==
==Etymology and place of tetration in the big picture of math==
Creation of word '''tetration''' is attributed to Englidh mathematician [[Reuben Louis Goodstein]]
Creation of word '''tetration''' is attributed to the English mathematician [[Reuben Louis Goodstein]]
<!-- This link seems to be wrong
<ref name="term">
<ref name="term">
"TETRATION, a term for repeated exponentiation, was introduced by Reuben Louis Goodstein".
"TETRATION, a term for repeated exponentiation, was introduced by Reuben Louis Goodstein".
Earliest Known Uses of Some of the Words of Mathematics,
Earliest Known Uses of Some of the Words of Mathematics,
http://members.aol.com/jeff570/t.html</ref>
http://members.aol.com/jeff570/t.html</ref>
!-->
<ref name="good">{{cite journal
<ref name="good">{{cite journal
|author= [[Reuben Louis Goodstein|R.L.Goodstein]]
|author= [[Reuben Louis Goodstein|R.L.Goodstein]]
Line 35: Line 68:
|volume=12
|volume=12
|year=1947}}
|year=1947}}
</ref>. This name indicates, that this operation is fourth (id est, [[tetra]]) in the hierarchy of operations after summation, multiplication, and exponentiation. In principle, one can define "pentation", "sexation", "septation" in the simlar manner, although tetration, perhaps, already has growth fast enough for the requests of XXI century.
</ref>.  


==Real values of the arguments==
The place of tetration in the [[mathematical analysis]] can be seen at the strong zoom-out of the [[big picture of math]].
Using [[mathematical notation]], the zoom-out of the [[mathematical analysis]] can be drawn as follows:
: <math>+\!+</math> has only one argument and means [[unitary increment]]
: <math>\mathrm{add}_b(1)= b+1=+\!+(b)</math> ; <math>\mathrm{add}_b(z\!+\!1)= b+(z\!+\!1)=+\!+\big(b\!+\!(z)\big)</math>
: <math>\mathrm{mul}_b(1)=b*1=b</math> ; <math>\mathrm{mul}_b(z\!+\!1)= b*(z\!+\!1)=b+\big(b*(z)\big)</math>
: <math> \exp_b(1)=b^1=b</math> ; <math> \exp_b(z\!+\!1)=b*(\exp_b(z))</math>
: <math>\mathrm{tet}_b(1)=b</math> ; <math> \mathrm{tet}_b(z\!+\!1)=\exp_b\big(\mathrm{tet}_b(z)\big)</math>
: <math>\mathrm{pen}_b(1)=b</math> ; <math> \mathrm{pen}_b(z\!+\!1)=\mathrm{tet}_b\big(\mathrm{pen}_b(z)\big)</math>
<!-- This table is quite similar to the [[Ackermann function]]s, but, historically, the Ackermann functions have some displacement in its argument and its value.
At least for integer parameters,
!-->
Except the zeroth row, each operation in the sequence above is just a recurrence of operations from the previous row.
Operation ++ could be called [[zeration]]  (although in [[programming language]]s it is called [[increment]]),
[[addition]] (or [[summation]]) could be called [[unation]],
[[multiplication]] (or [[product]]) could be called [[duation]] ,
[[exp|exponentiation]] could be called [[trination]].
The following operations ( tetration, [[pentation]]) have not been used so often, at least up to the year 2008. Although
tetration has been given many other names:
superexponentiation <ref name="w">
{{cite journal
|author= P.Walker
|title=Infinitely differentiable generalized logarithmic and exponential functions
|journal=[[Mathematics of computation]]
|volume=196
|year=1991
|pages=723-733
|url=http://www.jstor.org/pss/2938713
}}</ref>, ultraexponent <ref name="uxp">
{{cite journal
|url=http://www.tandfonline.com/doi/full/10.1080/10652460500422247
|author=M.H.Hooshmand.
|title=Ultra power and ultra exponential functions
|journal=[[Integral Transforms and Special Functions]]
|volume=17 (8)
|pages=549-558
|year=2006}}</ref>, generalized exponent <ref name="bromer"> N.Bromer. Superexponentiation. Mathematics Magazine, 60 No. 3 (1987), 169-174</ref>,
other names were not applied to the holomorphic extension of tetration, defined in the previous section.
Manipulation with the [[Holomorphism|holomorphic]] extensions and the inverses of [[summation]], [[multiplication]], [[exponentiation]] form the core of the [[mathematical analysis]].
The table above shows the place of tetration in the [[big picture of math]], in the penultimate row.
 
In the scheme above, each next operation appears as a [[superfunction]] with respect to the previous one.
In such a way, the name '''tetration''' indicates, that this operation is fourth (id est, [[tetra]]) in the hierarchy of operations after summation, multiplication, and exponentiation. In principle, one can define "pentation", "sextation", "septation" in a similar manner, although tetration, perhaps, already has a growth rate fast enough for the requirements of the 21st century.
 
==Real values of the arguments, general view==
Examples of behavior of this function at the real axis are shown in figure 1 for values  
Examples of behavior of this function at the real axis are shown in figure 1 for values  
<math>b=\mathrm{e}</math>,
<math>b=\mathrm{e}</math>,
<math>b=2</math>,
<math>b=2</math>,
<math>b=\exp(1/\mathrm{e})</math>, and for
<math>b=\exp(1\!/\!\mathrm{e})</math>, and for
<math>b=\sqrt{2}</math>. It has logarithmic singularity at <math>-2</math>, and it is monotonously increasing function.
<math>b=\sqrt{2}</math>. It has a logarithmic singularity at <math>-2</math>, and it is a [[Monotonic function|monotonic]] increasing function.


At <math>b\le \exp(1/\mathrm{e})</math> tetration <math>\mathrm{tet}_b(x)</math> approaches its limiting value as  
At <math>b\le \exp(1/\mathrm{e})</math> tetration <math>\mathrm{tet}_b(x)</math> approaches its limiting value as  
Line 48: Line 125:


===Fast growth and application===
===Fast growth and application===
At <math>b > \exp(1/\mathrm{e})</math> tetration <math>\mathrm{tet}_b(x)</math> grows faster than any exponential function. For this reason the tetration is suggested for the representation of huge numbers in [[mathematics of computation]].
For <math>~b\! >\! \exp(1\!/\!\mathrm{e})~</math> tetration <math>\mathrm{tet}_b(x)</math> grows faster than any exponential function. For this reason tetration has been proposed for the representation of huge numbers in the [[mathematics of computation]]<ref name="w">
A number, that cannot be stored as [[floating point]], could be represented as <math>\mathrm{tet}_b(x)</math> for some standard value of <math>b</math> (for example, <math>b=2</math> or <math>b=\mathrm{e}</math>) and relatively moderate value of <math>x</math>. The analytic properties of tetration could be used for the implementation of arithmetic operations with huge numbers without to convert them to the floating point representation.
<ref name="w">
P.Walker. Infinitely differentiable generalized logarithmic and exponential functions. Mathematics
P.Walker. Infinitely differentiable generalized logarithmic and exponential functions. Mathematics
of computation, 196 (1991), 723-733.</ref>
of computation, 196 (1991), 723-733.</ref>.
A number that cannot be stored as [[floating point]] could be represented as <math>\mathrm{tet}_b(x)</math> for some standard value of <math>b</math> (for example, <math>b=2</math> or <math>b=\mathrm{e}</math>) and relatively moderate value of <math>x</math>. The analytic properties of tetration could be used for the implementation of arithmetic operations with huge numbers without to convert them to the floating point representation.


===Integer values of the argument===
===Integer values of the argument===
For integer  <math>z</math>, tetration <math>{\rm tet}_b</math> can be interpreted as iterated exponential:
For integer  <math>z</math>, tetration <math>{\rm tet}_b(z)</math> can be interpreted as iterated exponential:
: <math>\mathrm{tet}_b(1)=b</math>
: <math>\mathrm{tet}_b(1)=b</math>
: <math>\mathrm{tet}_b(2)=b^b</math>
: <math>\mathrm{tet}_b(2)=b^b</math>
Line 69: Line 145:
99(1928), 118-133.</ref>
99(1928), 118-133.</ref>
:<math>\mathrm{tet}_2(n+3)=A(4,n)+3</math>
:<math>\mathrm{tet}_2(n+3)=A(4,n)+3</math>
where Ackermann function <math>A</math> is defined for the non-negative integer values of its arguments with equations
where the Ackermann function <math>A</math> is defined, for non-negative integer values of its arguments, by the following equations
:<math>
:<math>
A(m, n) =\left\{
A(m, n) =\left\{
Line 79: Line 155:
  \right.
  \right.
</math>
</math>
The generalization of the 4th Ackermann function <math>A(4,z)</math> for the complex values of <math>z</math> is described in the preprint <ref name="k2">
http://www.ils.uec.ac.jp/~dima/PAPERS/2008ackermann.pdf</ref>  . Construction of such holomorphic extension is equivalent to construction
of tetration <math>\mathrm{tet}_b(z)</math> for the base <math>b=2</math>.


==Asymptotic behavior and properties of tetration==
==Asymptotic behavior and properties of tetration==
The analytic extension of tetration <math>~F(z)</math> grows fast along the real axis of the complex <math>z</math>-plane, at least for some values of base <math>b</math>. However, it does not grow infinitely in the direction of imaginary axis.
The analytic extension of tetration <math>~F(z)</math> grows rapidly along the real axis of the complex <math>z</math>-plane, at least for some values of base <math>b</math>. However, it does not grow infinitely in the direction of imaginary axis.
The asymptotic behavior determines the basic properties of tetration.
The asymptotic behavior determines the basic properties of tetration.


The exponential convergence of discrete interation of logarithm corresponds to the exponential asymptotic behavior
The exponential convergence of discrete iteration of logarithm corresponds to the exponential asymptotic behavior
:<math> F(z)=L+\varepsilon(z) + o \big(\varepsilon(z)^{2}\big) </math>
:<math> F(z)=L+\varepsilon(z) + o \big(\varepsilon(z)^{2}\big) </math>
<!-- between cuts, $\Im(z)>0$, !-->
<!-- between cuts, $\Im(z)>0$, !-->
where
where
:<math> \varepsilon(z)=\exp(Qz+r) </math>,
:<math> \varepsilon(z)=\exp(Qz\!+\!r) </math>,
<math>Q</math> and  <math>r</math> are fixed complex numbers,
<math>Q</math> and  <math>r</math> are fixed complex numbers,
and <math>L</math> is eigenvalue of logarithm, solution of equation
and <math>L</math> is eigenvalue of logarithm, solution of equation
:<math> x=\log_{b}(x)</math>.
:<math> x=\log_{b}(x)</math>.


[[Image:ExampleEquationLog01.png|300px|right|thumb|FIg.2. Graphic solution of equation  
{{Image|ExampleEquationLog01.png|right|330px|FIg.2. Graphic solution of equation  
<math>x=\log_b(x)</math> for
<math>x\!=\!\log_b(x)</math> for
<math>b=\sqrt{2}</math> (two real solutions, <math>x=2</math> and <math>x=4</math>),
<math>b\!=\!\sqrt{2}</math> (two real solutions, <math>x\!=\!2</math> and <math>x\!=\!4</math>),
<math>b=\exp(1/\rm e)</math> (one real solution <math>x=\rm e</math>)
<math>b\!=\!\exp(1/\rm e)</math> (one real solution <math>x\!=\!\rm e</math>), and
<math>b=2</math> (no real solutions).]]
<math>b\!=\!2</math> (no real solutions).}}
 
Solutions of this equation are called [[fixed point]]s of logarithm.


Solutions of equation (14) are called [[fixed point]]s of logarithm.
===Fixed points of logarithm===
===Fixed points of logarithm===
Three examples of graphical solution of equation (14) are shown in figure 2 for  
Three examples of graphical solution of equation for fixed points of logarithm are shown in figure 2 for  
<math>b=\sqrt{2}</math>,  
<math>b\!=\!\sqrt{2}</math>,  
<math>b=\exp(1/\rm e)</math>,  and
<math>b\!=\!\exp(1\!/\!\rm e)</math>,  and
<math>b=2</math>.
<math>b\!=\!2</math>.


The black line shows function <math> y=x</math> in the <math>x,y</math> plane.
The black line shows function <math> y=x</math> in the <math>x,y</math> plane.
The colored curves show function <math> y=\log_b(x)</math> for cases  
The colored curves show function <math> y=\log_b(x)</math> for cases  
<math>b=\sqrt{2}</math>  (red),
<math>b\!=\!\sqrt{2}</math>  (red),
<math>b=\exp(1/\rm e)</math> (green),  and
<math>b\!=\!\exp(1/\rm e)</math> (green),  and
<math>b=2</math> (blue).
<math>b\!=\!2</math> (blue).


At  
At  
<math>b=\sqrt{2}</math>, there exist 2 solutions,  
<math>b\!=\!\sqrt{2}</math>, there exist 2 solutions,  
<math>x=2</math> and  
<math>x\!=\!2</math> and  
<math>x=4</math>.
<math>x\!=\!4</math>.


At  
At  
<math>b=\exp(1/\rm e)</math> there exist one solution  
<math>b\!=\!\exp(1\!/\!\rm e)</math> there exists one solution  
<math>x=\rm e</math>.
<math>x\!=\!\rm e</math>.


and <math>b=2</math>, there are no real solutions.<br>
and <math>b\!=\!2</math>, there are no real solutions.<br>


In general,
In general,
*at <math>b<\exp(1/\rm e)</math> there are two real solutions
*at <math>b\!<\!\exp(1\!/\!\rm e)</math> , there are two real solutions;
*at <math>b=\exp(1/\rm e)</math>, there is one soluition, and
*at <math>b\!=\!\exp(1\!/\!\rm e)</math> , there is one solution, and
*at <math>b>\exp(1/\rm e)</math> there esist two solutions, but they are complex.
*at <math>b\!>\!\exp(1\!/\!\rm e)</math> , there exist two solutions, but they are complex.


In particular,  
In particular,  
at  
at  
<math> b=\sqrt{2}</math>, the solutions are <br>
<math>b\!=\!\sqrt{2}</math>, the solutions are <br>
<math>x=L_{\sqrt{2},1}=2</math> and
<math>L\!=\!\!L_{\sqrt{2},1}=2</math> and
<math>x=L_{\sqrt{2},2}=4 </math> <br>.
<math>L\!=\!L_{\sqrt{2},2}=4 </math> <br>.


At  
At  
<math> b=2</math>, the solutions are <br>
<math>b\!=\!2</math>, the solutions are <br>
<math>x=L_2 \approx 0.824678546142074222314065+1.56743212384964786105857 \!~\rm i </math> and<br>
<math>L\!=\!L_2 \!\approx\! 0.824678546142074222314065+1.56743212384964786105857 \!~\rm i </math> and<br>
<math>x=L_2^*\approx 0.824678546142074222314065-1.56743212384964786105857 \!~\rm i </math>.
<math>L\!=\!L_2^*\!\approx\! 0.824678546142074222314065-1.56743212384964786105857 \!~\rm i </math>.


At <math> b=\rm e</math>, the solutions are<br>
At <math>b\!=\!\rm e</math>, the solutions are<br>
<math>x=L_{\rm e} \approx 0.318131505204764135312654+1.33723570143068940890116 \!~\rm i</math>
<math>L=L_{\rm e} \approx 0.318131505204764135312654+1.33723570143068940890116 \!~\rm i</math>
and<br>
and<br>
<math>x=L_{\rm e}^* \approx 0.318131505204764135312654-1.33723570143068940890116 \!~\rm i</math>.
<math>L=L_{\rm e}^* \approx 0.318131505204764135312654-1.33723570143068940890116 \!~\rm i</math>.


Few hundred straightforward iterations of equation (14) are sufficient to get the error smaller than the last decimal digit in the approximations above.
A few hundred straightforward iterations of equation <math>L=\log_b(L)</math> are sufficient to get the error smaller than the last decimal digit in the approximations above.


===Basic properties===
===Basic properties===
[[Image:TetrationAsymptoticParameters01.jpg|700px|right|thumb|FIg.3. parameters of asymptotic of tetration versus logarithm of the base]]
{{Image|TetrationAsymptoticParameters01.jpg|right|700px|FIg.3. parameters of asymptotic of tetration versus logarithm of the base}}
The solutions
The solutions
<math>x=L_1          </math> and  
<math>L=L_1          </math> and  
<math>x=L_2          </math>  of equation (14) are plotted in figure 3 versus  
<math>L=L_2          </math>  of equation (14) are plotted in figure 3 versus  
<math>\beta=\ln(b)</math> with thin black lines. Let  
<math>\beta=\ln(b)</math> with thin black lines.
At <math>\beta<1/\mathrm e</math>, both <math>L_1</math> and <math>L_2</math> are real and positive. Let  
<math>L_1<L_2      </math>, and only at
<math>L_1<L_2      </math>, and only at
<math>\ln(b)=1/e  </math>, the equality  
<math>\ln(b)=1/e  </math>, the equality  
Line 160: Line 241:
<math> \beta \ge 1/\rm e</math> represents the real part of the solutions  
<math> \beta \ge 1/\rm e</math> represents the real part of the solutions  
<math>L                          </math> and  
<math>L                          </math> and  
<math>L^*                      </math> of (14); the thin black dashed curve represents the two options for the imaginary part; the two solutions are complex conjugaitons of each other.  
<math>L^*                      </math> of (14); the thin black dashed curve represents the two options for the imaginary part; the two solutions are complex conjugations of each other.  
<!-- Let <math>\Re(L)>0</math>. !-->
<!-- Let <math>\Re(L)>0</math>. !-->
Requirement of definition of tetration determine the asymptotic of the solution. Parameter <math>Q</math> tetermines periodicity of quasi-periodicity of tetration. The two solutions for <\math>Q</math>are shown in figure 3 with green lines.
Requirement of definition of tetration determine the asymptotic of the solution. Parameter <math>Q</math> determines periodicity of quasi-periodicity of tetration. The two solutions for <math>Q</math> are shown in figure 3 with green lines.


At <math>b<\exp(1/ \mathrm{e} )</math> both solutons for <\math>Q</math> are real. The negative <math>Q</math>
At <math>b<\exp(1/ \mathrm{e} )</math> both solutions for <math>Q</math> are real. The negative <math>Q</math>
corresponds to tetration, decaying to the asymptotic value <math>L</math> in the direction of real axis; positive </math>Q</math> corresponds to the solution growing along the real axis. At the real axix, such a solution remains larger than unity; this does not allow to satisfy confition <math>F(0)=1</math>. Therefore, only one negative <math>Q</math> corresponds to the asymptotic behavior of tetration.
corresponds to tetration, decaying to the asymptotic value <math>L</math> in the direction of real axis; positive <math>Q</math> corresponds to the solution growing along the real axis. At the real axis, such a solution remains larger than unity; this does not allow to satisfy condition <math>F(0)=1</math>. Therefore, only one negative <math>Q</math> corresponds to the asymptotic behavior of tetration.


At <math>b>\exp(1/ \mathrm{e} )</math>, both options for <math>Q</math> are mutually complex conjugate. The real part is shown thif thick green line; one option of the imaginary part is shown with dashed line.
At <math>b\!>\!\exp(1/ \mathrm{e} )</math>, both options for <math>Q</math> are mutually complex conjugate. The real part is shown thif thick green line; one option of the imaginary part is shown with dashed line.


Possibilities for the period (or quasi-period) <math>T=\frac{2\pi}{Q}</math> are shown in Figure 3 with fotted lines. At  <math>b>\exp(1/ \mathrm{e} )</math>, only "negative" period forresponds to tetration. At <math>b>\exp(1/ \mathrm{ e} )</math>, the periodicity can be achieved only asymptotically; and <math>T</math> is quasi-period. The real part of quasi-period is markes with black dotted line; one of two options tor the imaginary part is marked with pink dotted line.
Possibilities for the period (or quasi-period) <math>T=\frac{2\pi}{Q}</math> are shown in Figure 3 with dotted lines. At  <math>b>\exp(1/ \mathrm{e} )</math>, only "negative" period corresponds to tetration. At <math>b>\exp(1/ \mathrm{ e} )</math>, the periodicity can be achieved only asymptotically; and <math>T</math> is quasi-period. The real part of quasi-period is marked with black dotted line; one of two options tor the imaginary part is marked with pink dotted line.


Generally, at <math>1<b<\exp(1/\mathrm{e})</math>, tetration is periodic; the period is pure imaginary.
Generally, at <math>1\!<\!b\!<\!\exp(1\!/\!\mathrm{e})</math>, tetration is periodic; the period is pure imaginary.


At <math>b=\exp(1/\mathrm{e})</math>, tetration is not periodic, and no exponential asymptotic exist.
At <math>b\!=\!\exp(1\!/\!\mathrm{e})</math>, tetration is not periodic, and no exponential asymptotic exist.


<math>b>\exp(1/\mathrm{e})</math>, tetration is quasi–periodic, the quasi-period in the upper complex half-plane is conjugate to that in the lower complex half-plane. The larger is base <math>b</math>, the shorter is quasi-period. As the quasi-periods are complex conjugated, the quasi-periodicity takes place away from the real axis.
<math>b\!>\!\exp(1\!/\!\mathrm{e})</math>, tetration is quasi–periodic, the quasi-period in the upper complex half-plane is conjugate to that in the lower complex half-plane. The larger is base <math>b</math>, the shorter is quasi-period. As the quasi-periods are complex conjugated, the quasi-periodicity takes place away from the real axis.


==Evaluation of tetration==
==Evaluation of tetration==
[[Image:AnalyticTetrationBaseSqrt2v00.jpg|220px|right|thumb|Fig.4. Tetration <math>f=F(z)</math> at the complex <math>z</math>-plane.
<!--
Levels of integer <math>\Re(f)</math> and levels of integer <math>\Im(f)</math> are shown with thick lines!-->.
<!--Grid covers with unit step range
(<math>-10\le \Re(z) \le 10 </math>,
<math>-10\le \Im(z) \le 10 </math>).!-->]]


As the asymptoric of tetration is crutually depend of base <math>b</math> in ficinity of  
As the asymptotics of tetration are critically dependent on the base <math>b</math> in the vicinity of the
value <math> b=\exp(1/\mathrm{e}) </math>, the evaluation procidure is different foe the cases
value <math> b=\exp(1/\mathrm{e}) </math>, the evaluation procedure is different for the cases
<math> 1<b<\exp(1/\mathrm{e}) </math>, <math> b=\exp(1/\mathrm{e}) </math>, and <math> b>\exp(1/\mathrm{e}) </math>, and should be considered intependently.
<math> 1\!<\!b\!<\!\exp(1\!/\!\mathrm{e}) </math>, <math> b\!=\!\exp(1/\mathrm{e}) </math>, and <math> b\!>\!\exp(1\!/\!\mathrm{e}) </math>, and these should be considered separately.


===Case <math> 1<b<\exp(1/\mathrm{e}) </math>===
===Case <math> 1<b<\exp(1/\mathrm{e}) </math>===
{{Image|HolomorphicTetrationBaseSqrt2v01.jpg|right|400px|Fig.4. Tetration <math>f=\mathrm{tet}_b(x\!+\!\mathrm{i}y)</math> for <math>b\!=\!\sqrt{2}</math>}}


At <math> 1<b<\exp(1/\mathrm{e}) </math>, the period <math>T</math> is imaginary.
For <math> 1<b<\exp(1/\mathrm{e}) </math>, the period <math>T</math> is imaginary.
Period with smallest modulus corresponds to the solution that is unity at the origen of coordinates. For <math> b=\sqrt{2}</math>, function <math>f=\mathrm{tet}_b(z)</math> is shown in figure 4 with levels of constant real part and levels of constant imaginary part. Levels
The period with smallest modulus corresponds to the solution that is unity at the origin of coordinates. For <math> b=\sqrt{2}</math>, the function <math>f=\mathrm{tet}_b(z)</math> is shown in figure 4 with levels of constant real part and levels of constant imaginary part. Levels
<math>\Re(f)=2,3,4</math> and levels  
<math>\Re(f)=2,3,4</math> and levels  
<math>\Im(f)=0,\pm 1,\pm2\pm 3, \pm 4</math> aew shown with thick lines. Intermediate levels are shown with thin lines. There are branch points at   
<math>\Im(f)=0,\pm 1,\pm2\pm 3, \pm 4</math> aew shown with thick lines. Intermediate levels are shown with thin lines. There are branch points at   
<math>\Re(z)=2~,~\Im(z)=2\pi |T| m~ \forall m \in \mathbb{N}</math>; the cut lines are
<math>\Re(z)=2~,~\Im(z)=2\pi |T| m~ \forall m \in \mathbb{N}</math>; the cut lines are
<math>\Re(z)<2~,~\Im(z)=2\pi |T| m~ \forall m \in \mathbb{N}</math>.  For this value of base, the period
<math>\Re(z)<2~,~\Im(z)=2\pi |T| m~ \forall m \in \mathbb{N}</math>.  For this value of the base, the period
:<math>T=\frac{2\pi}{\ln^2(2)}\approx  −17.1431481793548471041794 ~\mathrm{i}</math>.
:<math>T=\frac{2\pi}{\ln^2(2)}\approx  −17.1431481793548471041794 ~\mathrm{i}</math>.


The solution follows asymptotic at large values of real part of the argument, exponentially approaching the limiting value. In particular, for <math>b=\sqrt{2}</math>, this limiting value is 2.
There is a cut at  <math>x\!<\!-2</math>, <math>y\!=\!0</math>; although the jump at this cut reduces at the increase of <math>-x</math>. In such a way, the function approaches its limiting value <math>L\!=\!4</math> almost everywhere, although there is set of singularities at negative integer values of <math>x<-2</math>.
 
The solution follows asymptotic at large values of real part of the argument, exponentially approaching the limiting value. In particular, for <math>b\!=\!\sqrt{2}</math>, this maximum limiting value is in the left hand side of the figure, and to its minimum limiting value in the right hand side. For
<math>b\!=\!\sqrt{2}</math>, these limiting values are
<math>L\!=\!4</math> and
<math>L\!=\!2</math>.


The trace of the solution along the real axis corresponds to the red dotted curve in Figure 1.
The trace of the solution along the real axis corresponds to the red dotted curve in Figure 1.
Other solution of the recursive equaiton () do not satisfy criteria formulated in the definition of tetration.
Other solutions of the recursive equation <math>\exp_b\!\big(F(z)\big)\!=\!F(z\!+\!1)</math>,
that may grow up along the real axis, can be constructed in a similar way, but they do not satisfy the criteria formulated in the definition of tetration; in particular, <math>F(0)\!=\!1</math>.


===Case <math> b=\exp(1/\mathrm{e}) </math>===
===Case <math> b=\exp(1/\mathrm{e}) </math>===
[[Image:TetrationExp1e.jpg|right|thumb|220px|{{#ifexist:Template:TetrationExp1e.jpg/credit|{{TetrationExp1e.jpg/credit}}<br/>|}}Fig.5. Tetration at <math>b=\exp(1/\mathrm{e})</math>.]]
[[Image:HolomorphicTetrationBaseExp1e01.jpg|right|400px|thumb|
<!--
{{#ifexist:Template:TetrationExp1e.jpg/credit|{{TetrationExp1e.jpg/credit}}<br/>|}}
!-->
Fig.5. Tetration at <math>b=\exp(1\!/\!\mathrm{e})</math>.]]


At <math> b=\exp(1/\mathrm{e}) </math>, asymptotically
At <math> b=\exp(1\!/\!\mathrm{e}) </math>, the limiting value <math>L=\mathrm{e}</math>, and, asymptotically,
: <math> F(z)=\mathrm{e}- 2 \mathrm{e}/z + \mathcal{O}(1/z^2) </math>
: <math> F(z)=\mathrm{e}- 2 \mathrm{e}/z + \mathcal{O}\!\big(\ln(z)/z^2\big) </math>


The function <math>f=F(x+\mathrm{i}y)</math> is shown in figure 5.
The function <math>f=F(x\!+\!\mathrm{i}y)</math> is shown in figure 5.


Levels <math>\Re(f)=-2,-1,0,1,2,3,4</math> are shown with thick black lines.  
Levels <math>\Re(f)=-2,-1,0,1,2,3,4</math> are shown with thick black lines.  
Line 220: Line 306:
Intermediate levels are shown with thin lines.
Intermediate levels are shown with thin lines.


There is cut at  <math>x<-2</math>, <math>y=0</math>.
There is a cut at  <math>x<-2</math>, <math>y=0</math>, but the hump of the function at the cut reduces as
<math>y</math> reduces, id est, with increasing <math>|z|</math>. In such a way, everywhere, at
<math>|z|\rightarrow \infty</math>, the function approaches its limiting value <math>L=\mathrm{e}\approx 2.71</math>.
almost everywhere, although there is a set of singularities at negative integer values of <math>z<-1</math>.


Behavior of this function at real values of argument is shown in Figure 1 with thin solid line.
Behavior of this function at real values of argument is shown in Figure 1 with thin solid line.
Other solutions of the recursive equation <math>\exp_b\!\big(F(z)\big)\!=\!F(z\!+\!1)</math> can be constructed in the similar way; they may grow up along the real axis, but they do not satisfy criteria formulated in the definition of tetration; in particular, <math>F(0)\!=\!1</math>.


===Case <math> b>\exp(1/\mathrm{e}) </math>===
===Case <math> b>\exp(1/\mathrm{e}) </math>===
[[Image:TetrationBase2andE.jpg|right|thumb|400px|{{#ifexist:Template:TetrationBase2andE.jpg/credit|{{TetrationBase2andE.jpg/credit}}<br/>|}}Fig.6. Tetration  <math>\mathrm{tet}_b(x+\mathrm{i}y)</math>at base <math>b=2</math> and <math>b=\mathrm e</math>.]]
{{Image|TetrationBase2andE.jpg|right|400px|Fig.6. Tetration  <math>\mathrm{tet}_b(x\!+\!\mathrm{i}y)</math> at base <math>b\!=\!2</math> and <math>b\!=\!\mathrm e</math>.}}
At <math> b>\exp(1/\mathrm{e}) </math>, tetration is asymptotically periodic. It exponentially
At <math> b>\exp(1/\mathrm{e}) </math>, tetration is asymptotically periodic. It decays exponentially
decays to the fixed points <math>L</math> and <math>L^*</math> in th eupper and lower halfplane. This allows to express it through its values along the imaginary axis, using the [[Cauchi integral]]. <ref name="k">D.Kouznetsov. Solutions of <math>F(z+1)=\exp(F(z))</math> in the complex <math>z</math>plane. [[Mathematics of Computation]], 2008, in press; preprint: http://www.ils.uec.ac.jp/~dima/PAPERS/analuxp99.pdf</ref>.  
to the fixed points <math>L</math> and <math>L^*</math> in the upper and lower halfplane. This allows to express it through its values along the imaginary axis, using the [[Cauchi integral]]. <ref name="k">D.Kouznetsov. Solutions of <math>F(z+1)=\exp(F(z))</math> in the complex <math>z</math>plane. [[Mathematics of Computation]], 2008, in press; preprint: http://www.ils.uec.ac.jp/~dima/PAPERS/analuxp99.pdf</ref>.  


For the base <math>b=2</math> and <math>b=\mathrm{e}</math> tetration  
For the base <math>b=2</math> and <math>b=\mathrm{e}</math> tetration  
<math>f=\mathrm{tet}_b(z)</math> is plotted in figure 6. Levels  
<math>f=\mathrm{tet}_b(z)</math> is plotted in figure 6. Levels  
<math>\Re{f}=0, \pm 1, ... \pm 16 </math> and
<math>\Re{f}=0, \pm 1, ... \pm 16 </math> and
<math>\Im{f}=0, \pm 1, ... \pm 16 </math> are drawn with thick lines. The function has [[logarithmic singulatiry]] at point -2 and cut at real values of the argument, smaller than -2. In the right hand side, symbols <math>\approx \infty</math> mean huge values that cannot be stored in the conventional [[floating point]] representation (logarithm, [[mantissa]]). In the upper left and lower left part of eath graphic, the function approaches its asymptotic values <math>L</math> and <math>L^*</math>. Function is quasi-periodic; the same fractal structure  
<math>\Im{f}=0, \pm 1, ... \pm 16 </math> are drawn with thick lines. The function has a [[logarithmic singularity]] at point -2 and cut at real values of the argument, smaller than -2. In the right hand side, symbols <math>\approx \infty</math> mean huge values that cannot be stored in the conventional [[floating point]] representation (logarithm, [[mantissa]]). In the upper left and lower left part of each graphic, the function approaches its asymptotic values <math>L</math> and <math>L^*</math>. Function is quasi-periodic; the same fractal structure reproduces again and again at the translation of argument with quasiperiod <math>T</math> in the upper halfplane and
reproduces again and again at the translation of argument with quasiperiod <math>T</math> in the upper halfplane and
<math>T^*</math> in the lower halfplane.
<math>T^*</math> in the lower halfplane.


Along the real axis, tetration for this values of base is plotted also in figure 1 with thick solid and dashed lines.
There is cut at  <math>x<-2</math>, <math>y=0</math>. The jump of the function at this cut approaches
<math>L\!-\!L^*</math>
almost everywhere, although there is set of singularities at negative integer values of <math>x<-2</math>.
 
Along the real axis, tetration for these values of the base is plotted also in figure 1 with thick solid and dashed lines.
 
===Behavior along the real axis===
{{Image|TetrationDerivativesReal.jpg|right|280px|Derivatives of tetration at <math>b=\mathrm{e}</math>.}}
The growth of tetration along the real axis is crucially determined by its base.
The graphic of this function is shown at the top of the article for <math>b=\sqrt{2}, \exp(1/\mathrm{e}), 2, \mathrm{e}</math>.
For <math>b=\mathrm{e}</math>, the derivatives of tetration are plotted in figure at the right.
: <math> \mathrm{tet}^{\prime}</math> is plotted with red;
: <math> \mathrm{tet}^{\prime\prime}</math> is plotted with green;
: <math> \mathrm{tet}^{\prime\prime\prime}</math> is plotted with blue;
: <math> \mathrm{tet}^{\prime\prime\prime\prime}</math> is plotted with pink.
 
Tetration is strictly growing function; its first derivative is positive. For <math>b=\mathrm{e}</math> the minimum of the derivative takes place in vicinity of <math>-0.5</math> and is slightly smaller than unity.
At <math> 1\! <\!b \!< \exp(1 \! / \mathrm{e} )</math>, the growth is limited by the minimum of the limiting values
<math>L</math>. Tetration approaches this limiting value exponentially.
In particular, at <math>b=\sqrt{2}</math>, the limiting value is 2.
 
At <math> b \! =\! \exp(1 \! / \! \mathrm{e} )</math>, the growth is limited by the fixed point
<math>L\!=\!\mathrm{e}</math>. Tetration approaches this limiting value as rational function.
 
<!--
Only at base<math> b \! >\! \exp(1 \! / \! \mathrm{e} )</math>, tetration shows growth faster than any exponential along the real axis.
!-->
 
The growing-up holomorphic solution of equation <math>\exp_b(F(z))\!=\!F(z\!+\!1)</math> can be constructed in the similar way also at
<math>1\! \le \! b \! \le \! \exp(1\!/\!{\mathrm e})</math>, but at the real axis, such a solution remains larger than unity, and the condition <math>F(0)\!=\!1</math> cannot be satisfied. Therefore, this solution cannot be interpreted as extension of iterative exponentiation <math>\exp_b\Big( \exp_b\big(... 1 ..)\big) </math> for non-integer number <math>z</math> of exponentiations; in this sense, such a solution is not a tetration.
 
Tetration <math>\mathrm{tet}_b</math> shows explosively-fast growth along the real axis
<!-- and quasi-periodic behavior in the complex plane !-->
only at values <math>b \! >\! \exp(1 \! / \! \mathrm{e} )</math>.
 
==Tetration at the base <math>b=\sqrt{2} </math>==
{{Image|TetrationQ2r.jpg|left|200px|Tetration at the base <math>b=\sqrt{2} </math>}}
There is hypothesis that at base <math>b=\sqrt{2} </math>, the graphic of tetration <math>y=\mathrm{tet}_b(x)</math> is symmetric with respect to the line <math>y=-x</math>. This line, together with the graphic, is shown in Figure at let. The asymptotics <math>x=-2</math> and <math>y=2</math> are also shown. The graphic looks symmetric with respect to
<math>x \leftrightarrow y</math>, but no proof of this hypothesis was suggested.
Now this hypothesis is negated; the apparent symmetry is only approximation. However,  for
<math>b=\sqrt{2} </math>, at the segment
<math>-1\le x\le 0</math>, the relation
<math>\mathrm{tet}_b(x)\approx -\mathrm{tet}_b^{-1}(-x)</math> holds with 4 decimal digits
<ref>
D.Kouznetsov, H.Trappmann. Portrait of the four regular super-exponentials to base sqrt(2). Mathematics of Computation, 2010, v.79, p.1727-1756.
http://www.ams.org/journals/mcom/2010-79-271/S0025-5718-10-02342-2/home.html
http://mizugadro.mydns.jp/PAPERS/2010q2.pdf
</ref><ref>https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0
http://www.ils.uec.ac.jp/~dima/BOOK/202.pdf
http://mizugadro.mydns.jp/BOOK/202.pdf Д.Кузнецов. Суперфункции. Lambert Academic Press, 2014. (In Russian, 328 pages.)
</ref>.
 
==Tetration at base e==
Holomorphic tetration on the natural base <math>b\!=\!\mathrm{e}</math> is the most developed, at least up to the year 2008.  In the rest of this article, it is assumed that <math>b\!=\!\mathrm{e}\approx 2.71828</math> although the majority of results allow a straightforward extension to the case of real <math>b\!>\!\exp(1/\mathrm{e})\approx 1.44467</math>.
 
Tetration has real [[fixed point]] <math>X_0</math>, id est, solution <math>z</math> of the equation
<math>z=\mathrm{tet}(z)</math>. Its approximation
<!--
:<math>X_0\approx 1.8503545290271812</math> !-->
:<math>X_0\approx 1.85035452902718</math>
can be found, iterating equation <math>z=\mathrm{slog}(z)</math>, where slog is inverse function of tetration.
 
==Pentation==
[[File:Penmap.jpg|left|240px|thumb|[[Complex map]] of natural [[pentation]], <math>u\!+\!\mathrm i v=\mathrm{pen}(x\!+\!\mathrm i y)</math> by
<ref name="lambert">
https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0
http://mizugadro.mydns.jp/BOOK/202.pdf
Д.Кузнецов. Суперфункции. [[Lambert Academic Publishing]], 2014. (In Russian)
</ref>]]
[[Superfunction]] of [[tetration]] can be called [[pentation]]; it is solution of equations
: <math>\mathrm{pen}(z\!+\!1)=\mathrm{tet} \big( \mathrm{pen}(z)\big) </math>
: <math>\mathrm{pen}(0)=1</math>
Expansion of [[tetration]] in vicinity of its fixed point <math>X_0</math> from the interval <math>(-2,-1)</math> (for real base <math>b\!>\!1</math>)
allows to define unique (and the only "true") pentation. Other superfunctions of tetration can be obtained by the periodic modification of the argument of pentation.
 
For base <math>b=\mathrm e\approx 2.71</math>, the fixed point of tetration
<math>X_0\!\approx\! -1.850354529</math>.
and [[pentation]] to this base approaches this value at large negative values of the real part of the argument
<ref>
D.Kouznetsov. Holomorphic ackermanns. http://mizugadro.mydns.jp/PAPERS/2014acker.pdf under consideration
</ref><ref name="lambert">
https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0
http://mizugadro.mydns.jp/BOOK/202.pdf
Д.Кузнецов. Суперфункции. [[Lambert Academic Publishing]], 2014. (In Russian)
</ref>. Complex map of [[natural pentation]] is shown in figure at left;  <math>u\!+\!\mathrm i v=\mathrm{pen}(x\!+\!\mathrm i y)</math>.


==Inverse of tetration==
==Inverse of tetration==
[[Image:Kslog00.jpg|right|thumb|350px|{{#ifexist:Template:Kslog00.jpg/credit|{{Kslog00.jpg/credit}}<br/>|}}Fig.7. kslog in the complex plane.]]
[[Image:Kslog00.jpg|right|thumb|350px|
The inverse of tetration <math>t=\mathrm{kslog}(z)=\mathrm{tet}^{-1}(z)</math> can be performed using the [[Newton method]], solving equation <math> \mathrm{tet}(t)=z</math>, leading to  
<!--
{{#ifexist:Template:Kslog00.jpg/credit|{{Kslog00.jpg/credit}}<br/>|}}
!-->
Fig.7. slog in the complex plane.]]
The inverse of tetration <math>t=\mathrm{slog}(z)=\mathrm{tet}^{-1}(z)</math> can be performed using the [[Newton method]], solving equation <math> \mathrm{tet}(t)=z</math>, leading to  
:<math> \mathrm{tet}\big(\mathrm{kslog}(z)\big)=z</math>.
:<math> \mathrm{tet}\big(\mathrm{kslog}(z)\big)=z</math>.
The inverse function has branchpints <math>L</math> and <math>L^*</math>. For the kslog, at base <math>b=\mathrm{e}</math>, shown in the figure 7, the cuts are placed horisontaly, along the lines
The inverse function has branchpoints <math>L</math> and <math>L^*</math>. For the kslog, at base <math>b=\mathrm{e}</math>, shown in the figure 7, the cuts are placed horizontaly, along the lines
: <math>\Re(z)<\Re(L)</math>,  <math>\Im(z)=\pm \Im(L)</math>.
: <math>\Re(z)<\Re(L)</math>,  <math>\Im(z)=\pm \Im(L)</math>.
Due to the symmetry <math>\mathrm{kslog}\big(z^*\big)=\mathrm{kslog}\big(z\big)^*</math>, it is sufficient to plot only half of the complex plane.
Due to the symmetry <math>\mathrm{kslog}\big(z^*\big)=\mathrm{kslog}\big(z\big)^*</math>, it is sufficient to plot only half of the complex plane.


[[Image:TetrationBiholomorphic.jpg|right|thumb|200px|{{#ifexist:Template:TetrationBiholomorphic.jpg/credit|{{TetrationBiholomorphic.jpg/credit}}<br/>|}} Fig.8. Shaded region is image of the upper halfplane at the mapping with kslog function.]]
[[Image:TetrationBiholomorphic.jpg|right|200px|thumb|
<!--
{{#ifexist:Template:TetrationBiholomorphic.jpg/credit|{{TetrationBiholomorphic.jpg/credit}}<br/>|}}  
!-->
Fig.8. Shaded region is image of the upper halfplane at the mapping with slog function.]]


The mapping with function kslog is shown in figure 8. For <math>x</math> in the shaded region, the relation
The mapping with function kslog is shown in figure 8. For <math>x</math> in the shaded region, the relation
:<math> \mathrm{kslog}\big(\mathrm{tet}(z)\big)=z</math>
:<math> \mathrm{slog}\big(\mathrm{tet}(z)\big)=z</math>
takes place. The upper part of the complex plane is mapped into the upper halfplane, and the lower halfplane is mapped to the lower halfplane. The real axis is mapped into the halfline <math>(-2,\infty)</math>. The fixed points
takes place. The upper part of the complex plane is mapped into the upper halfplane, and the lower halfplane is mapped to the lower halfplane. The real axis is mapped into the halfline <math>(-2,\infty)</math>. The fixed points
<math>L</math> is mapped to imaginary infinity, following the shadower strip.
<math>L</math> is mapped to imaginary infinity, following the shaded strip.


In figure 8, the images of the grid lines <math>\Re(z)=-1,0,\Re(L),1,2,3</math> and
In figure 8, the images of the grid lines <math>\Re(z)=-1,0,\Re(L),1,2,3</math> and
Line 258: Line 439:
<math>\Im(F(z))=\rm const</math> shown in figure 6 for <math>b=\rm e</math>.
<math>\Im(F(z))=\rm const</math> shown in figure 6 for <math>b=\rm e</math>.


Together, the pair of functions tet and kslog allow to evaluate any iteration (including negative, fractional and even complex) of the exponential finction.
Together, the pair of functions tet and kslog allow to evaluate any iteration (including negative, fractional and even complex) of the exponential function.
 
===Beyond the cutlines of slog===
{{Image|SuperlogarithmCutlines.jpg|right|900px|3 options to put cutlines of the superlogarithmic function}}
The function slog has branchpoints, and there are many options to pit the cutlines. The cutlines, parallel to the abscias axis, considered above is one of many options. One can put the cutlines along the level <math>\Re(\mathrm{slog}(z))=-2</math>. This cases simplify the plotting of the slog function, while the only tetration is available for the efficient evaluation. On the other hand, there is no simple expression for the parametrization of such a cutline, and the calculation of the position of the cutlines slows down the algorithm of evaluation of slog. In addition,
there is tradition, that cutlines of the special functions are placed parallel to the real or imaginary axis. One could place the cutline from the branchpoints to the abscissa axis. In this case, the additional cutlines should go along the negative part of the real axis, which is not convenient for the applications. (For example, for the evaluation of the generalized exponential <math>\exp^c</math> at non-integer <math>c</math>, it is important, that the real axis belongs to the range of holomorphism of slog. There are 3 other options to place the cutlines parallel to the coordinate axes, keeping condition <math>\mathrm{slog}(z^*)=\mathrm{slog}(z)^*</math>. There 3 options are shown in the figure.
 
The shaded region shows the domain
:<math>G= \{ z \in \mathbb{C}~:~ \Re(z)\!>\!\Re(L)~,~ |z|<|L| \}</math>
For evaluation of any of these superlogarithms, it is sufficient to have the efficient algorithm for the evaluation in this domain.
 
The left hand side picture corresponds to the slog, sa it is defined above. Within the strip between the cutlines, the function approaches its limiting value <math>-2</math>; it becomes infinite in vicinity of the branch points and varies very slow in the rest of the complex plane.
 
The central part of the picture represents slog with vertical cutlines, let us call this function slogv. In the left hand side of the complex halfplane, the evaluation of slog begins with exponentiation of its argument. The exponential is periodic function; therefore, in this part, slog is periodic:
:<math>\mathrm{slogv}(z+2\pi\mathrm{i} m)= \mathrm{slogv}(z)~ \forall  m\in \mathbb{N} ~,~ \forall z\in \mathbb{C}~:~ \Re(z)<0  </math>
 
The third option, let us denote it slogr, is to put the cutlines in the direction of the real axis. Along the real axis, the function remains holomorphic, but there are additional branchpoints in the ranges <math>\Re(z)>1~,~ 1<|\Im(z)|<2</math>. In these regions, the function slogr show the fractal behavior.
 
The slog, defined at the beginning, seems to be simpler than slogv and slogr; therefore, namely slog with horizontal cutlines, parallel to the real axis, going in the direction, opposite to the real axis, is used below.
 
==Polynomial approximation==
The [[Taylor series]] for the tetration can be written in the usual form:
:<math>\mathrm{tet}(z)=\sum_{n=0}^{\infty} c_n (z-z_0)^n</math>
where the <math>n</math>th coefficient
:<math> c_n=\frac{1}{n!}\mathrm{tet}^{(n)}(z_0)</math>
is expressed through the <math>n</math> derivative of the function.
The coefficients of the expansion can be calculated using the straightforward differentiation of the representation through the [[Cauchy integral]].
 
===Taylor expansion at zero===
{{Image|TetrationPolynomial25power.jpg|right|600px|Fig.N. Approximation of tetration with polynomial of 25th power}}
For <math> z_0=0</math>, the calculation gives the following values
:<math>c</math><sub>0</sub> =  1
:<math>c</math><sub>1</sub> ≈  1.091767351258322138
:<math>c</math><sub>2</sub> ≈  0.271483212901696469
:<math>c</math><sub>3</sub> ≈  0.212453248176258214
:<math>c</math><sub>4</sub> ≈  0.069540376139988952
:<math>c</math><sub>5</sub> ≈  0.044291952090474256
:<math>c</math><sub>6</sub> ≈  0.014736742096390039
:<math>c</math><sub>7</sub> ≈  0.008668781817225539
:<math>c</math><sub>8</sub> ≈  0.002796479398385586
:<math>c</math><sub>9</sub> ≈  0.001610631290584341
:<math>c</math><sub>10</sub>≈  0.000489927231484419
:<math>c</math><sub>11</sub>≈  0.000288181071154065
:<math>c</math><sub>12</sub>≈  0.000080094612538551
:<math>c</math><sub>13</sub>≈  0.000050291141793809
:<math>c</math><sub>14</sub>≈  0.000012183790344901
:<math>c</math><sub>15</sub>≈  0.000008665533667382
:<math>c</math><sub>16</sub>≈  0.000001687782319318
:<math>c</math><sub>17</sub>≈  0.000001493253248573
:<math>c</math><sub>18</sub>≈  0.000000198760764204
:<math>c</math><sub>19</sub>≈  0.000000260867356004
:<math>c</math><sub>20</sub>≈  0.000000014709954143
:<math>c</math><sub>21</sub>≈  0.000000046834497327
:<math>c</math><sub>22</sub>≈-0.000000001549241666
:<math>c</math><sub>23</sub>≈  0.000000008741510781
:<math>c</math><sub>24</sub>≈-0.000000001125787310
:<math>c</math><sub>25</sub>≈  0.000000001707959267
<!--
,-0.000000000377858315
, 0.000000000349577877
,-0.000000000105377012
, 0.000000000074590971
,-0.000000000027175982
, 0.000000000016460766
,-0.000000000006741873
, 0.000000000003725329
,-0.000000000001639087
, 0.000000000000858364
!-->
 
The truncated Taylor series gives the [[polynomial approximation]].
In the upper right hand side of the Figure N, the polynomial
:<math>f=\sum_{n=0}^{25} c_n z^n</math>
is shown in the complex <math>z</math> plane.
:Levels <math>\Re(f)=-2,-1,0,1,2,3,4</math>  are shown with thick black curves.
:Levels <math>\Re(f)=-1.8,-1.6,-1.4,-1.2,-0.8,-0.6,-0.4,-0.2</math>  are shown with thin red curves.
:Levels <math>\Re(f)=0.2,0.4,0.6,0.8,1.2,1.4,1.6,1.8</math>  are shown with thin thin blue curves.
:Levels <math>\Im(f)=-2,-1</math>  are shown with thick red curves.
:Levels <math>\Im(f)=1,2</math>  are shown with thick blue curves.
:Levels <math>\Im(f)=\pm \pi,\pm 3\pi</math>  are shown with thick pink curves.
:Levels <math>\Im(f)=
\pm 0.2,
\pm 0.4,
\pm 0.6,
\pm 0.8,
\pm 1.2,
\pm 1.4,
\pm 1.6,
\pm 1.8</math> are shown with thin green curves.
 
In the upper left corner of figure N, the same is shown for function <math>\log(f)</math>
 
At the bottom left, the overlap of the upper two images is shown.
 
At the bottom right, lines of constant modulus and constant phase of holomorphic tetration in the same range.
 
The good approximation of tetration takes place in the range <math>|z|</math> of order of unity or smaller; the radius of convergence of the series is 2.
 
===Expansion at 3i ===
[[Image:TetrationTailorExpansion3ipower25.jpg|500px|thumb|Taylor expansion of tetration developed at 3i, truncated at 30th power, plotted in the complex plane.]]
Coefficients of the expansion
:<math>\mathrm{tet}(z)=\sum_{n=0}^{\infty} a_n \left(\frac{z-3 \mathrm{i}}{2}\right)^n</math>
can be evaluated in the similar way:
:<math>a_{ 0}\approx 0.3709065890322851 + 1.3368216707889140 i</math>
:<math>a_{ 1}\approx 0.0366009653759846 + 0.1392221538995050 i</math>
:<math>a_{ 2}\approx-0.1688843184064154 + 0.0971853361962927 i</math>
:<math>a_{ 3}\approx-0.1268131504868087  -0.1183162876702863 i</math>
:<math>a_{ 4}\approx 0.0423580931032393  -0.1052093008832072 i</math>
:<math>a_{ 5}\approx 0.0584830639356318  -0.0081022452449608 i</math>
:<math>a_{ 6}\approx 0.0234003166529485 + 0.0180777701182038 i</math>
:<math>a_{ 7}\approx 0.0034426098470138 + 0.0181510375563591 i</math>
:<math>a_{ 8}\approx-0.0080369581444167 + 0.0091742846703500 i</math>
:<math>a_{ 9}\approx-0.0070469552816877  -0.0009395850672747 i</math>
:<math>a_{10}\approx-0.0018461796309531  -0.0032234258318168 i</math>
:<math>a_{11}\approx 0.0005406488544310  -0.0018967206101561 i</math>
:<math>a_{12}\approx 0.0010224364808881  -0.0005596865717924 i</math>
:<math>a_{13}\approx 0.0006471439639805 + 0.0002598066193583 i</math>
:<math>a_{14}\approx 0.0001044445559337 + 0.0003719947259883 i</math>
:<math>a_{15}\approx-0.0001117853540434 + 0.0001678668755219 i</math>
:<math>a_{16}\approx-0.0001063015871081 + 0.0000207220003313 i</math>
:<math>a_{17}\approx-0.0000507809881911  -0.0000357591300574 i</math>
:<math>a_{18}\approx-0.0000031474299869  -0.0000352318593759 i</math>
:<math>a_{19}\approx 0.0000134766134413  -0.0000133303413745 i</math>
:<math>a_{20}\approx 0.0000098023908240 + 0.0000004760718415 i</math>
:<math>a_{21}\approx 0.0000035549347545 + 0.0000038981621220 i</math>
:<math>a_{22}\approx-0.0000002155265265 + 0.0000029627341324 i</math>
:<math>a_{23}\approx-0.0000013167390363 + 0.0000009738135453 i</math>
:<math>a_{24}\approx-0.0000008340196081  -0.0000001866385871 i</math>
:<math>a_{25}\approx-0.0000002286961098  -0.0000003749771677 i</math>
:<math>a_{26}\approx 0.0000000537258461  -0.0000002306013659 i</math>
:<math>a_{27}\approx 0.0000001140665665  -0.0000000656951029 i</math>
:<math>a_{28}\approx 0.0000000666359546 + 0.0000000232663057 i</math>
:<math>a_{29}\approx 0.0000000139678685 + 0.0000000331511830 i</math>
:<math>a_{30}\approx-0.0000000068489056 + 0.0000000171304198 i</math>
 
The plot of approximation of tetration with the resulting polynomial of 30th power is shown in figure. This approximation can be used for plotting of camera-ready pictures of tetration, using it and its conjugation at <math>\Re(z)\le 0.5 </math>. With 50 terms, at <math>|z-3|<2</math>, such approximation returns 14 significant figures.
 
==Asymptotic expansion at large values of the imaginary part of the argument==
{{Image|TetrationAsymptoticExpansion00.jpg|right|500px|Fig.As. Deviation of tetration from its asymptotic expansion.}}
Here only the case of base <math>b\!=\!\mathrm{e}</math> is considered; although the generalization to the case <math>b\!>\!\mathrm{e}</math> is straightforward. In this section, the behavior of tetration
<math>\mathrm{tet}(z)</math> is considered for moderate values of
<math>\Re(z)</math> and large values of <math>\Im(z)</math>
 
At <math>b\!=\!\mathrm{e}</math>, in the upper halfplane, tetrations approach the fixed point
<math>L\approx 0.3+1.3\mathrm{i}</math> of the logarithm. This approach is exponential. Using the exponential
:<math>\varepsilon=\varepsilon(z)=\exp(L z +R) </math>
as a small parameter, for some complex constant <math>R</math>, the tetration can be estimated as follows:
: <math>\mathrm{tet}(z)=L+\varepsilon+a_2 \varepsilon^2 + a_3 \varepsilon^3 + \mathcal{O}(\varepsilon^4)</math>
Substitution of this expression into the equation <math>\mathrm{tet}(z\!+\!1)=\exp\big(\mathrm{tet}(z)\big)</math> gives <br>
<math>a_2=\frac{1/2}{L-1}\approx -0.151314 - 0.2967488 \mathrm{i}~,~</math>
<math>a_3=\frac{a_2+1/6}{L^2-1}\approx -0.36976 + 0.98730 \mathrm{i}</math><br>
Value <math>R\approx 1.077961437526 -0.946540963949 \mathrm{i}</math> can be found fitting the tetration
<math>\mathrm{tet}(z)</math> with function <math> L+\epsilon+a_2 \epsilon^2 + a_2 \epsilon^3 </math> at large values of
<math>\Im(z)</math>. Then, the expansion can be written in the form
:<math>\mathrm{tet}(z)=\sum_{n=0}^N a_n \varepsilon^n + \mathcal(O)\!\big(\varepsilon^{N\!+\!1}\big)</math>
assuming that <math>a_0\!=\!L</math> and <math>a_1\!=\!1</math>. At given <math>N</math>, such a representation
indicates the possible approximations of tetration. The deviations of tetration from these approximations are shown in figure Fig.As.
 
The four plots in fig.As correspond to the four asymptotic approximations. The deviations
:<math>f=\mathrm{tet}(z)-L</math>, 
:<math>f=\mathrm{tet}(z)-\Big(L+\varepsilon(z)\Big)</math>, 
:<math>f=\mathrm{tet}(z)-\Big(L+\varepsilon+a_2 \varepsilon^2 \Big)</math>, 
:<math>f=\mathrm{tet}(z)-\Big(L+\varepsilon+a_2 \varepsilon^2 + a_3 \varepsilon^3\Big)</math>
are shown in the complex plane with lines of constant phase and constant modulus.
:Levels  <math>\arg(f)=-2,-1</math> are shown with thick red lines.
:Level  <math>\arg(f)=0</math> is shown with thick black lines.
:Levels  <math>\arg(f)=1,2</math> are shown with thick blue lines.
:Levels  <math>\arg(f)=\pm \pi</math> are shown with scratched lines. (these lines reveal the step of sampling used by the plotter).
 
:Levels  <math>|f|=\exp(-0.8),\exp(-0.6),\exp(-0.4),\exp(-0.2)</math> are shown with thin red lines.
:Levels  <math>|f|=\exp(0.2),\exp(0.4),\exp(0.6),\exp(-0.8)</math> are shown with thin blue lines.
:Levels  <math>|f|=\exp(3), \exp(2), \exp(1),\exp(0), \exp(-\!1), \exp(-\!2), \exp(-\!3), \exp(-\!4), \exp(-\!5), \exp(-\!5), \exp(-\!7), \exp(-\!8)</math> are shown with thin thick black lines.
:Level <math>|f|=\exp(-10)</math> is shown with thick red line.
:Levels <math>|f|= \exp(-12),\exp(-14), \exp(-16),\exp(-18)</math> are shown with thick black lines.
:Level <math>|f|=\exp(-20)</math> is shown with thick red line.
:Levels <math>|f|= \exp(-22),\exp(-24), \exp(-26),\exp(-28)</math> are shown with thick black lines.
:Level <math>|f|=\exp(-30)</math> is shown with thick green line.
:Level <math>|f|=\exp(-31)</math> is shown with thick black line.
The plotter tried to draw also
:Level <math>|f|=\exp(-32)</math> with thick black line and
:Level <math>|f|=\exp(-33)</math> with thin dark green line, which are seen a the upper left hand side corners of the two last pictures, but the precision of evaluation of tetration is not sufficient to plot the smooth lines; for the same reason, the curve for
:<math>|f|=\exp(-31)</math> in the last picture, in the upper right band side looks a little bit irreguler; also, the pattenn in the upper left corner of the last two pictures looks chaotic; the plotter cannot distinguish the function from its asymptotic approximation.
 
The figure indicates that, at <math>\Im(z)>4 </math>, <math>\Re(z)<4\Im(z) - 25</math>, the asymptotic approximation
:<math>\mathrm{tet}(z)\approx L+\varepsilon+a_2 \varepsilon^2 + a_3 \varepsilon^3</math>
gives at least 14 correct significant figures. At large values of the imaginary part, this approximation is more precise than the evaluation of tetration through the contour integral.
 
==Approximation of tetration with [[elementary function]]s==
{{Image|TetrationApproLP100.jpg|right|300px|<math>f(z)=\log(z+2)+\sum_{n=0}^{100}c_n z^n</math> in the complex
<math>z</math>-plane.}}
Due to recurrent relation <math>\mathrm{tet}(z)=\exp(\mathrm{tet}(x+1))</math>, it is sufficient to approximate tetration in any vertical strip of unit with in the complex plane. Some of such approximations are suggested in <ref name="k">k</ref>.
In principle, the numerical approximation of tetration with implementation of the [[Cauchi integral]] with finite sums <ref name="k">k</ref> also should be considered as approximation with elementary function. However for the practical evaluation of tetration, shorter expressions are more suitable. One of such approximation comes from the [[Taylor series|Taylor expansion]] of function <math>\mathrm{tet}(z)-\log(z+2)</math>. The subtraction of logarithm remove the closest singularity that limits the radius of convergence of the Taylor series, and makes precise the approximation with finite sum. One of such approximation with one logarithm and polynomial of 100th power
:<math>f(z)=\log(z+2)+\sum_{n=0}^{100}c_n\left(\frac{z}{2}\right)^n</math>
is shown at the figure in the complex <math>z</math>-plane.
<!-- The good approximation (say 12 significant figures) takes place at <math>|z|\le2</math>. !-->
In vicinity of the origin of coordinates,
say, <math>|z|<1 </math>, the most of terms are negligibly small, and the shortened sum still approximates the tetration.
The first coefficients in this expansion are
:<math>c_0 = 0.30685281944005469058 </math>
:<math>c_1 = 1.18353470251664338875</math>
:<math>c_2 = 1.58593285160678321155 </math>
:<math>c_3 = 1.36629265207672068172 </math>
:<math>c_4 = 1.36264601823980036066  </math>
:<math>c_5 = 1.21734246689515424045 </math>
:<math>c_6 = 1.10981816083559525765 </math>
:<math>c_7 = 0.96674692974769849130 </math>
:<math>c_8 = 0.84089872598668435888 </math>
:<math>c_9 = 0.71353210966804747617 </math>
:<math>c_{10}=0.60168548504001373445 </math>
More coefficients can be extracted from the generator of the figure.
While <math>|z| \le 2 </math>, the approximation with 101 terms returns at least 14 correct significant figures.
 
==Approximation of slog==
Function slog, which is inverse of tetration, allows the approximation with elementary functions.
 
===Taylor expansion===
{{Image|SLOGtailor16.jpg|right|400px|Approximation of slog with polynomial of 16th power from the Taylor expansion at unity.}}
 
The [[Taylor series]] for the tetration can be inverted, gaining the expansion of the superlogarithm:
: <math>\mathrm{slog}(z)=\sum_{n=1}^\infty  c_n (z-1)^n</math>.
 
Approximations for the first 16 coefficients:
:<math>c_{ 1}\approx  0.91594605649953  </math>
:<math>c_{ 2}\approx -0.20861842957759 </math>
:<math>c_{ 3}\approx -0.05450400630209 </math>
:<math>c_{ 4}\approx  0.07134941925273  </math>
:<math>c_{ 5}\approx -0.02004387374438 </math>
:<math>c_{ 6}\approx -0.01101258023037 </math>
:<math>c_{ 7}\approx  0.01207268318645  </math>
:<math>c_{ 8}\approx -0.00272922880760 </math>
:<math>c_{ 9}\approx -0.00269905319156 </math>
:<math>c_{10}\approx  0.00243941500632 </math>
:<math>c_{11}\approx -0.00036220360858</math>
:<math>c_{12}\approx -0.00070125921262</math>
:<math>c_{13}\approx  0.00052782155380 </math>
:<math>c_{14}\approx -0.00002987943551</math>
:<math>c_{15}\approx -0.00018614540434</math>
:<math>c_{16}\approx  0.00011722843042 </math>
The partial sum with 16 terms (from zero to 16) is plotted in the figure in the complex <math>z</math> plane.
Lines of constant real part and constant imaginary part are drawn.
 
The radius of convergence of this series is determined by the distance to the closest singularity; the representation of
<math>\mathrm{slog}(z)</math> with the Taylor series is valid for
: <math>|z-1|< |1-L|</math>
Obviously, it fails at <math>z=L</math>. For this case, the asymptotic expansion can be used.
 
===Expansion at fixed point <math>L</math> of logarithm===
{{Image|SlogFitFixedPoint04.jpg|right|500px|Fitting of slog with the expansion around <math>L</math>.}}
 
Superlogarithm can be approximated with expansion <ref name="kneser">
H.Kneser. “Reelle analytische L¨osungen der Gleichung <math>\varphi(\varphi(x)) = e^x </math> und verwandter Funktionalgleichungen”.
Journal fur die reine und angewandte Mathematik, 187 (1950), 56-67.</ref>
:<math>\mathrm{slog}(z)=\frac{1}{L} \left( \log(z\!-\!L)+
\sum_{m,n} r_{m,n} (z\!-\!L)^{m+2\pi \mathrm{i} n/L}
\right)</math>,
where <math>L</math> is [[fixed point]] of logarithm. This expansion indicates the ways to construct the asymptotic approximations of slog.
The coefficients can be expressed from the asymptotical analysis of equation <math>\mathrm{slog}(\exp(z))= \mathrm{slog}(z)+1</math>. Also, they can be expressed from the asymptotical estimate of tetration at large values of the imaginary part of the argument.
The evaluation of first coefficients gives
:<math>r_{0,0}\approx \!  -\!1.0779614375 +  0.9465409639 i </math>
:<math>r_{1,0}\approx  ~ 0.1513148972 +  0.2967488367 i </math>
:<math>r_{2,0}\approx  \! -\!0.0607692016 +  0.0359770148 i </math>
:<math>r_{3,0}\approx  \! -\!0.0079378875  -0.0176412865 i </math>
:<math>r_{4,0}\approx  ~ 0.0051087546  -0.0000718839 i </math>
These coefficients allow to approximate slog in vicinity of the fixed point of logarithm with function
:<math>f=\frac{1}{L} \left( \log(z\!-\!L)+
\sum_{m=0}^4 r_{m,0} (z-L)^{m}
\right)</math>.
In the figure this <math>f</math> is shown in the complex <math>z</math>-plane.
:Levels <math>\Re(f)=0,\pm 1, \pm2</math> are shown with thick black lines.
:Levels <math>\Re(f)=-1.9, .. -0.1</math> are shown with thin red lines.
:Levels <math>\Im(f)=-1.9, .. -1.9</math> are shown with thin dark green lines.
:Levels <math>\Im(f)=0</math> is shown with thick green line.
(Deviation of this line from the real axis indicates the error of the approximation.)
:Level <math>\Im(f)=-1</math> is shown with thick red line
:Levels <math>\Im(f)=1,2,3</math> are shown with thin dark blue lines.
For comparison, dashed lines show the precise evaluation for some of the levels above for the robust implementation of the slog function as inverse of tetration. While <math>|z-L|<1 </math>, the deviation of these dashed lines from the levels for function <math>f</math> is not seen even at the strong zooming-in of the central part of the figure.
===Approximation of slog with elementary functions===
{{Image|SLOGappro50.jpg|right|400px|Approximation of slog with function fslog.}}
{{Image|SLOGtest50.jpg|right|600px|Numerical test of approximation of slog with function fslog.}}
 
The precision of approximation of slog (with fixed precision of the arithmetics used for the tvaluation) can be extended, takung unto account the singularitues of slog at the fixed points. From the asumptotical representation above, one can guess the robust representation for slog:
:<math> \mathrm{slog}(z)=\frac{\log(z-L)}{L} + \frac{\log(z-L^*)}{L^*} + \sum_{n=0}^\infty u_n (z-1)^n</math>
The coefficients of this expansion are real. The first coefficients:
:<math>u_0\approx\! ~1.419225215505</math>
:<math>u_1\approx\! -0.026066290298</math>
:<math>u_2\approx~ 0.001733047818</math>
:<math>u_3\approx\! -0.000019521307</math>
:<math>u_4\approx\! -0.000063070064</math>
:<math>u_5\approx~ 0.000025678960</math>
:<math>u_6\approx\! -0.000005590100</math>
:<math>u_7\approx\! -0.000000072797</math>
:<math>u_8\approx~ 0.000000651489</math>
:<math>u_9\approx\! -0.000000276981</math>
:<math>u_{10}\!\approx~ 0.000000031118</math>
:<math>u_{11}\!\approx~ 0.000000029409</math>
:<math>u_{12}\!\approx\! -0.000000018969</math>
 
This representation allows construction of approximations, truncating the series. One of such approximations
:<math> f=\mathrm{fslog}(z)=\frac{\log(z-L)}{L} + \frac{\log(z-L^*)}{L^*} + \sum_{n=0}^{50} u_n (z-1)^n</math>
is shown in figure in the complex <math>z</math> plane. Lines of constant <math>\Re(f)</math> and those of constant <math>\Im(f)</math> are plotted.
 
The range of approximation of slog with this function is wider than that with the Taylor expansion at unity. The extended range of approximation allows its validation with the numerical test of identities
: <math>\mathrm{slog}(z)=\mathrm{slog}(\exp(z))-1</math>
: <math>\mathrm{slog}(z)=\mathrm{slog}(\log(z))+1</math>
The residuals
: <math>f=\mathrm{fslog}(z)-\mathrm{fslog}(\exp(z))+1</math> and
: <math>f= \mathrm{fslog}(z)-\mathrm{fslog}(\log(z))-1</math>
are shown in the figure with levels <math>|f|=</math>const.
In the voided regions in vicinity of <math>z=0</math> and <math>z=1</math>, the residual is at the level of <math>10^{-14}</math>. (It is difficult to make the residual smaller, using the arithmetics with double complex variables.) This test indicates, that at <math>|z-1| \le |L|</math>, the approximation of slog with two the logarithms and the polynomial of 50th power gives at least 9 correct significant figures.


==Iterated exponential and <math>\sqrt{\exp}</math>==
==Iterated exponential and <math>\sqrt{\exp}</math>==
[[Image:Sqrt(exp)(z).jpg|right|500px|thumb|Fig.9. <math>\exp^c(z)</math> in the complex <math>z</math> plane for various <math>c</math>.]]
{{Image|Sqrt(exp)(z).jpg|right|500px|Fig.9. <math>\exp^c(z)</math> in the complex <math>z</math> plane for various <math>c</math>.}}


Especially interesting is the case of iteration of natural exponential, id est, <math> b=\mathrm{e}</math>. Existence of the fractional iteration, and, in particular, existence of operation <math>\sqrt{\exp}=\exp^{1/2}</math> was demonstrated in 1950 by H.Kneser. <ref name="kneser">H. Kneser.
===Pre-historic approach to the iterated exponential===
Reelle analytische Losungen der Gleichung <math>\varphi(\varphi(x))=\exp(x)</math> und verwandter Funktionalgleichungen''. Journal fur die reine und angewandte Mathematik, v.187 (1950), 56-67.
Especially interesting is the case of iteration of natural exponential, id est, <math> b=\mathrm{e}</math>. Existence of the fractional iteration, and, in particular, existence of operation <math>\sqrt{\exp}=\exp^{1/2}</math> was demonstrated in 1950 by H.Kneser. <ref name="kneser">H. Kneser. Reelle analytische Losungen der Gleichung <math>\varphi(\varphi(x))=\exp(x)</math> und verwandter Funktionalgleichungen''. Journal fur die reine und angewandte Mathematik, v.187 (1950), 56-67.
</ref>.
</ref>.
However, that time, there was no computer facility for the evlauation of such an exotic function
However, that time, there was no computer facility for the evaluation of such an exotic function
<math>F</math> that <math>F(F(z))=\exp(z)</math>; perhaps, just  absence  of an apropriate plotter did not allow Kneser to plot the distribution of fractal exponential function  
<math>F</math> that <math>F(F(z))=\exp(z)</math>; perhaps, just  absence  of an apropriate plotter did not allow Kneser to plot the distribution of fractal exponential function  
<math>\exp^c(z)</math> in the complex <math>z</math> plane for various values of <math>c</math>, shown in Fig.9.
<math>\exp^c(z)</math> in the complex <math>z</math> plane for various values of <math>c</math>, shown in Fig.9.


[[Image:Expc.jpg|right|200px|thumb|Fig.10. <math>\exp^c(x)</math> versus <math>x</math>  
{{Image|Expc.jpg|right|300px|Fig.10. <math>\exp^c(x)</math> versus <math>x</math>  
<!-- for some real values of <math>c</math> !-->]]
<!-- for some real values of <math>c</math> !-->}}


===The Implementation through the tetration===
Holomorphic tetration allows to extend the iterated exponential <math>\exp_b^c</math>
Holomorphic tetration allows to extend the iterated exponential <math>\exp_b^c</math>


Line 287: Line 777:
<math>\exp_b^1=\exp_b</math>. If the subscript is omitted, it is assumed to be <math>\mathrm{e}</math>, id est, <math>\exp^c=\exp_\mathrm{e}^c</math>
<math>\exp_b^1=\exp_b</math>. If the subscript is omitted, it is assumed to be <math>\mathrm{e}</math>, id est, <math>\exp^c=\exp_\mathrm{e}^c</math>


===Iterated exponential in the complex plane===
Function <math>f=\exp^c(z)</math> is shown in figure 9  with levels of constant real part and levels of constant imaginary part. Levels  
Function <math>f=\exp^c(z)</math> is shown in figure 9  with levels of constant real part and levels of constant imaginary part. Levels  
<math>\Re(f)=-3,-2,-1,0,1,2,3,4,5,6,7,8,9</math> and
<math>\Re(f)=-3,-2,-1,0,1,2,3,4,5,6,7,8,9</math> and
Line 299: Line 790:
At non-integer values of <math>c</math>, <math>L</math> and <math>L^*</math> are [[branch point]]s of function <math>\exp^c</math>; in figure, the cut is placed parallel to the real axis. At <math>c<0</math> there is an additional cut which goes along the negative part of the real axis. In the figure, the cuts are marked with pink lines.
At non-integer values of <math>c</math>, <math>L</math> and <math>L^*</math> are [[branch point]]s of function <math>\exp^c</math>; in figure, the cut is placed parallel to the real axis. At <math>c<0</math> there is an additional cut which goes along the negative part of the real axis. In the figure, the cuts are marked with pink lines.


===Iterated exponential of a real argument===
For real values of the argument, function <math> y=\exp^c(x)</math> is ploted in figure 10 versus <math>x</math> for values
For real values of the argument, function <math> y=\exp^c(x)</math> is ploted in figure 10 versus <math>x</math> for values
<math>c=0,\pm 0.1, \pm 0.5, \pm 0.9, \pm 1, \pm 2</math>.<br>
<math>c=0,\pm 0.1, \pm 0.5, \pm 0.9, \pm 1, \pm 2</math>.<br>
Line 306: Line 798:
<math>\exp^{-1}(x)=\ln(x)</math>,
<math>\exp^{-1}(x)=\ln(x)</math>,
<math>\exp^{-2}(x)=\ln\big(\ln(x)\big)</math> and so on.
<math>\exp^{-2}(x)=\ln\big(\ln(x)\big)</math> and so on.
At least for the real and big enough <math>x</math>, the relation <math>\exp^u\big(\exp^v(x)\big)=\exp^{u+v}(x)</math> holds, which is qute analogous of relations
<math>ux\!+\!vx=(u\!+\!v)x</math> and <math>x^u x^v=x^{u+v}</math>.
<!-- and <math> \big(x^u\big)^v=x^(uv)</ref>.!-->
However, at negative <math>u</math> or negative <math>v</math>, value <math>x</math> should be big enough, that 
<math>\exp^u(x)</math> and
<math>\exp^v(x)</math> and
<math>\exp^{u+v}(x)</math> are defined, see figure 10. For example, at <math>\min(u,v,u\!+\!v)=-2</math>, we need <math>x\!>\!1</math>.
In particular, at <math>u\!=\!1</math>, <math>v\!=\!-1</math>, we have
<math>\exp^1\!\big(\exp^{-1}(x)\big)=\exp\!\big(\ln(x)\big)=x</math>.
===Application of iterated exponential===
The iterated exponential, that can be implemented with holomorphic tetration, may have various applications. In particular,
The <math>\exp_b^c(x)</math> at <math>0\!<\!c\!<\!1</math> could describe a process that grows faster than any polynomial, but slower than any exponential. In such a way, the iterated exponential, at the proper implementation, should greatly extend the abilities of fast and precise [[fitting of functions]]. This is just analogy of function <math>f(x)=x^a</math> which, at fractal values of <math>a</math>, may be good for description of a function that grows faster than any linear function but slower than any quadratic function.


==Similar functions==
==Similar functions==
WIthdrowal of some of requirements from the definition of tetration allows the huge variety of similar functions.
Withdrawal of some of requirements from the definition of tetration allows the huge variety of similar functions.
 
===Entire solutions of <math>F(z\!+\!1)=\exp(F(z))</math>===
Withdrawal of the requirement <math>F\big(z^*\big)=F(z)^*</math> and <math>F(0)\!=\!1</math>
allows the solution by Kneser <ref name="kneser">kneser</ref>, which is entire and also could be used to build up various powers of the exponential; in particular, <math>\sqrt{\exp}</math>. Such entire function is shown in upper part of figure 1c in <ref name="k">k</ref>, in order to reveal the asymptotic behavior of holomotphic tetration.
 
Withdrawal of condition <math>F(0)\!=\!1</math> allows to construct solutions, similar to the growing tetration, for base <math>1\!<\!b\!\le\!\exp(1/\mathrm{e})</math>. Although such solutions cannot be interpreted as generalization of exponential iterated <math>z</math> times, they can be useful for generalization of [[exponential function]].


Withdrowal of the requirement <math>F(z^*)=F(z)^*</math> and <math>F(0)=1</math>
===Non-holomorphic modification of tetration===
allows the solution by Kneser <ref name="kneser">kneser</ref>, which is entire and also could be used to buld up carious powers of the exponential.  
{{Image|ZpluxSinZ.jpg|right|300px|Fig.11. Almost identical function <math>f(z)\!=\! z \!+\! 10^{-9} \sin(2 \pi z)</math> in the complex <math>z\!=\!x\!+\!\mathrm{i} y</math> plane.}}
{{Image|TetrationModified.jpg|right|400px|Fig.12. Motified tetration at the complex plane.}}
{{Image|TetrationModifiedZoom.jpg|right|400px|Fig.13. Zoom in of fig.12}}
Withdrawal of requirement of holomorphicity from the definition of tetration allows functions, which look like the tetration, at least along the real axis. Even the reduction of the range of holomorphism in the requirement allows to consider tetration with modified argument. The modified tetration can be defined as
: <math> \mathrm{tem}(z)=\mathrm{tet}\!\big(J(z)\big)= \mathrm{tet}\!\big(z\!+\! h(z)\big)</math> ,
where <math>J(z)=z+h(z)</math>, and <math>h</math> is a 1-periodic function. The simple example of such function is
:<math>h(z)=10^{-9}\sin(2 \pi z)</math>
In this case, along the real axis, the function <math>J</math> is almost identical to its argument; and values of the modified tetration are close to values of tetration. Being plotted at figure 1 or in figure 10, the deviation of such function <math>J</math> from the identity is small, and the deviation of the modified tetration from tetration is also small. If the figures are printed in the real scale, then the deviation of the curves would be of order of atomic size.


WIthdrowal of requirement of holomorphicity and boundness allow huge variety of functions. They can be expressed in the form
However the difference becomes seen at the complex values of the argument.  
<math> \tilde F(z)= F(z+h(z))</math> , where <math>h</math> is 1-periodif function. Howevr, the entire function
In figure 11, function <math>J(z)=z\!+\!h(z)</math> is plotted in the complex <math>z\!=\!x\!+\!\mathrm{i} y</math> plane.
<math>g(z)=z+h(z)</math> should take all possible values, except, perhaps, one. Therefore, there are singularities of tetration among these values. In such a way, any other solution of recursive tetration <math>F(z+1)=\exp_b(z)</math> should have singularities inside the strip <math>|\Re(z)|< 1 </math>.  
Levels of constant real part and those of constant imaginary part are drawn. In vicinity of the real axis, these lines almost coincide with the gridline; the grid is drawn with step unity and extended one step to the right and one step to the left from the graphic. In order to show that it behaves as it if would be a continuation of the plot. At <math>|\Im(z)|\approx 3</math>, the deviation from the identical function becomes visible, and at
<math>|\Im(z)|>3</math>, the <math>J(z)</math> has many points with real values, including those with various negative integer values.
The tetration has cut at negative values of <math>z<-2</math> and singularities there. Therefore, the cuts of the modified tetration are determined by the lines <math>J(z) \le -2 </math>, and modified tetration unavoidable has singularities in points <math>z</math> such that <math>J(z) \in \{ n\in \mathbb{N} ~:~ n<-1 \} </math>.
These singularities are determined by the function <math>h</math> and do not depend on the base <math>b</math> of tetration. In figure 11, the lines <math>J(z)<-2</math> are seen not only along the real axis, but also at the top and at the bottom of the figure. In such a way, figure 11 shows the cutlines of the modified tetration. One has no need to evaluate tetration in order to
see the margin of the change of holomorphism of the modified tetration.


WIthdrowal of condition <math>F(0)=1</math> allows to construct solutions, similar to the growing tetration, for base <math>1<b<\exp(1/\mathrm{e})</math>. Although such solutions cannot be interpreted as generalization of exponential iterated </math>z</math> times, they can be useful for generalization of [[exponential function]].
In figure 12, the modified tetration is plotted in the complex <math>z</math>-plane.
The additional cuts are seen in the upper and lower parts of the figure 12. Only within the strip along the real axis, the function is holomorphic. While the amplitude of sinusoidal is of order of <math>10^{-9}</math>, the strip of holomorphism is wider than unity, although this width slightly reduces along the real axis.


==See also==
In order to see the behavior of the modified tetration in vicinity of the additional singularities, in fig.13, the zooming-in of the part of figure 12 is shown. The zoom has improved resolution, so, in its turn, it can be zoomed in to the size of the screen of a computer to see the details. In each cell of the grid, the small and deformed image of the central part of the fig.12 appears.
*[[exp]]
 
*[[holomorphic function]]
One has no need to evaluate tetration in order to reveal its singularities outside the real axis.
*[[biholomorphism]]
All the solutions <math>z</math> of equation <math>z+h(z)=n</math> for integer <math>n<-1</math> are singularities (branchpoints) of the modified tetration.
*[[Ackermann functions]]
 
*[[fixed point]]
Let <math>\mathcal{M}= \{ n\in \mathbb{N}~:~n<-1 \} </math>.
 
The following theorem is suggested:
For any entire 1-periodic function <math>h</math> such that
<math>h(0)\!=\!0</math>, <math>h</math> is not a constant,
<math>h\big(z^*\big)=h(z)^*~\forall z\in \mathbb{C} </math>,
there exist <math> z\in \mathbb{C} : |\Re(z)|\!<\!1~,~ z\!+\!h(z)~ \in \mathcal{M} </math>.
sequence
Although this theorem is not yet proved, the intents to construct at least one example of function, that would contradict it, were not successful. This theorem is somehow independent from the theory of tetration, but it indicates, that any modified tetration cannot
be holomorphic in the range <math>\{z\in \mathrm{C}~:~\Re(z)>-2 \} </math>. 
 
According to the theorem above, the modified tetration does not satisfy the condition of quasi-periodicity, and does not satisfy the criterion of holomorphism in the definition of tetration. The sequence of cutlines for the specific example of modified tetration is seen in figure 13; the modified tetration is not even continuous. An addition to function <math>h</math> some highest sinusoidals brings the discontinuities even closer to the real axis. This indicates, that if in at least one point at the real axis between <math>-\!1</math> and <math>1</math>, some solution <math>F</math> of equations <math>F(z\!+\!1)\!=\!\exp(F(z))</math>, <math>F(0)\!=\!1</math> differs from tetration tet for at least <math>10^{-9}</math>, then this solution is not holomorphic in the range <math>\{z\in \mathbb{C} : |\Re(z)|<1, |\Im(z)|<4 \}</math>.
 
Even small deformation of tetration tet breaks its continuity. Similar reasons in favor of uniqueness tetration are suggested also in <ref name="k">k</ref>. There is only one tetration, that satisfies requirements of the definition, although the rigorous mathematical proof of the uniqueness is still under development.
 
<!--
In such a way, any other solution of recursive tetration <math>F(z\!+\!1)\!=\!\exp_b(z)</math> should have singularities inside the strip <math>|\Re(z)|< 1 </math>.
For any real parameter <math>a</math> the range <math>a\!-\!1\!<\! \Re(z)\! \le \!a\! </math> can be filled with ANY function (for example, the linear function),  and extended to the whole complex plane using the recurrent equation. The resulting recurrent function may have a lot of cuts and singularities; it may be very different from the holomorphic tetration outside the real axis.!-->


==References==
==References==
<references/>
{{reflist}}
 
==Free online sources==
*http://en.wikipedia.org/wiki/Tetration
*http://reglos.de/lars/ffx.html , Lars Kindermann. References about Iterative Roots and Fractional Iteration.
*http://math.eretrandre.org/tetrationforum/index.php , Discussion about tetration
*http://tetration.itgo.com/ Tetration site by Andrew Robbins
*http://www.tetration.org/ Tetration site by Daniel Geisler


[[Category:Mathematical functions]]
[[Category:Suggestion Bot Tag]]
[[Category:Functions of complex variable]]

Latest revision as of 06:00, 26 October 2024

This article is developing and not approved.
Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
Code [?]
 
This editable Main Article is under development and subject to a disclaimer.

Tetration is a rapidly growing mathematical function, which was introduced in the 20th century and proposed for the representation of huge numbers in the mathematics of computation. For positive integer values of its argument , tetration on base can be defined with:

Fig.1. Tetration for various versus .

For real values of the argument and various values of the base , this is plotted in Fig.0.

Fig.1. Tetration in the plane with levels .

The map of in the plane is shown in Fig.1 with levels .

Up to year 2011, this function has not been listed among elementary functions, it is not implemented in programming languages and it is not used for the internal representation of data in computers.

In this article, the generalizaiton of tetration for complex (and, in particular, real) values of its argument is described. At base , tetration is assumed to be a holomorphic function, at least for positive values of the real part of its argument. This tetration is used to construct the holomorphic extension of the iterated exponential for the case of non-integer values of the number of iterations.

Definition

For real , Tetration on the base is a function of a complex variable, which is holomorphic at least in the range , bounded in the range , and satisfies conditions

at least within the range .

According to this definition, tetration is superfunction of the exponential. This justifies the alternative name "superexponential" for this function and "superlogarithm" for the inverse function. The definition above generalizes the definitions, recently suggested for the specific cases of base [1] and [2].

Etymology and place of tetration in the big picture of math

Creation of word tetration is attributed to the English mathematician Reuben Louis Goodstein [3].

The place of tetration in the mathematical analysis can be seen at the strong zoom-out of the big picture of math. Using mathematical notation, the zoom-out of the mathematical analysis can be drawn as follows:

has only one argument and means unitary increment
 ;
 ;
 ;
 ;
 ;

Except the zeroth row, each operation in the sequence above is just a recurrence of operations from the previous row. Operation ++ could be called zeration (although in programming languages it is called increment), addition (or summation) could be called unation, multiplication (or product) could be called duation , exponentiation could be called trination. The following operations ( tetration, pentation) have not been used so often, at least up to the year 2008. Although tetration has been given many other names: superexponentiation [4], ultraexponent [5], generalized exponent [6], other names were not applied to the holomorphic extension of tetration, defined in the previous section.

Manipulation with the holomorphic extensions and the inverses of summation, multiplication, exponentiation form the core of the mathematical analysis. The table above shows the place of tetration in the big picture of math, in the penultimate row.

In the scheme above, each next operation appears as a superfunction with respect to the previous one. In such a way, the name tetration indicates, that this operation is fourth (id est, tetra) in the hierarchy of operations after summation, multiplication, and exponentiation. In principle, one can define "pentation", "sextation", "septation" in a similar manner, although tetration, perhaps, already has a growth rate fast enough for the requirements of the 21st century.

Real values of the arguments, general view

Examples of behavior of this function at the real axis are shown in figure 1 for values , , , and for . It has a logarithmic singularity at , and it is a monotonic increasing function.

At tetration approaches its limiting value as , and .

Fast growth and application

For tetration grows faster than any exponential function. For this reason tetration has been proposed for the representation of huge numbers in the mathematics of computation[4]. A number that cannot be stored as floating point could be represented as for some standard value of (for example, or ) and relatively moderate value of . The analytic properties of tetration could be used for the implementation of arithmetic operations with huge numbers without to convert them to the floating point representation.

Integer values of the argument

For integer , tetration can be interpreted as iterated exponential:

and so on; then, the argument of tetration can be interpreted as number of exponentiations of unity. From definition it follows, that

and

Relation with the Ackermann function

At base , tetration is related to the Ackermann function [7]

where the Ackermann function is defined, for non-negative integer values of its arguments, by the following equations

The generalization of the 4th Ackermann function for the complex values of is described in the preprint [2] . Construction of such holomorphic extension is equivalent to construction of tetration for the base .

Asymptotic behavior and properties of tetration

The analytic extension of tetration grows rapidly along the real axis of the complex -plane, at least for some values of base . However, it does not grow infinitely in the direction of imaginary axis. The asymptotic behavior determines the basic properties of tetration.

The exponential convergence of discrete iteration of logarithm corresponds to the exponential asymptotic behavior

where

,

and are fixed complex numbers, and is eigenvalue of logarithm, solution of equation

.
© Image: Dmitrii Kouznetsov
FIg.2. Graphic solution of equation for (two real solutions, and ), (one real solution ), and (no real solutions).

Solutions of this equation are called fixed points of logarithm.

Fixed points of logarithm

Three examples of graphical solution of equation for fixed points of logarithm are shown in figure 2 for , , and .

The black line shows function in the plane. The colored curves show function for cases (red), (green), and (blue).

At , there exist 2 solutions, and .

At there exists one solution .

and , there are no real solutions.

In general,

  • at , there are two real solutions;
  • at , there is one solution, and
  • at , there exist two solutions, but they are complex.

In particular, at , the solutions are
and
.

At , the solutions are
and
.

At , the solutions are
and
.

A few hundred straightforward iterations of equation are sufficient to get the error smaller than the last decimal digit in the approximations above.

Basic properties

FIg.3. parameters of asymptotic of tetration versus logarithm of the base

The solutions and of equation (14) are plotted in figure 3 versus with thin black lines. At , both and are real and positive. Let , and only at , the equality takes place. Basic properties of tetration are determined by the base . The main parameters versus are plotted in figure 3. The thin black solid curve at represents the real part of the solutions and of (14); the thin black dashed curve represents the two options for the imaginary part; the two solutions are complex conjugations of each other. Requirement of definition of tetration determine the asymptotic of the solution. Parameter determines periodicity of quasi-periodicity of tetration. The two solutions for are shown in figure 3 with green lines.

At both solutions for are real. The negative corresponds to tetration, decaying to the asymptotic value in the direction of real axis; positive corresponds to the solution growing along the real axis. At the real axis, such a solution remains larger than unity; this does not allow to satisfy condition . Therefore, only one negative corresponds to the asymptotic behavior of tetration.

At , both options for are mutually complex conjugate. The real part is shown thif thick green line; one option of the imaginary part is shown with dashed line.

Possibilities for the period (or quasi-period) are shown in Figure 3 with dotted lines. At , only "negative" period corresponds to tetration. At , the periodicity can be achieved only asymptotically; and is quasi-period. The real part of quasi-period is marked with black dotted line; one of two options tor the imaginary part is marked with pink dotted line.

Generally, at , tetration is periodic; the period is pure imaginary.

At , tetration is not periodic, and no exponential asymptotic exist.

, tetration is quasi–periodic, the quasi-period in the upper complex half-plane is conjugate to that in the lower complex half-plane. The larger is base , the shorter is quasi-period. As the quasi-periods are complex conjugated, the quasi-periodicity takes place away from the real axis.

Evaluation of tetration

As the asymptotics of tetration are critically dependent on the base in the vicinity of the value , the evaluation procedure is different for the cases , , and , and these should be considered separately.

Case

Template:CC-Image
Fig.4. Tetration for

For , the period is imaginary. The period with smallest modulus corresponds to the solution that is unity at the origin of coordinates. For , the function is shown in figure 4 with levels of constant real part and levels of constant imaginary part. Levels and levels aew shown with thick lines. Intermediate levels are shown with thin lines. There are branch points at ; the cut lines are . For this value of the base, the period

Failed to parse (syntax error): {\displaystyle T=\frac{2\pi}{\ln^2(2)}\approx −17.1431481793548471041794 ~\mathrm{i}} .

There is a cut at , ; although the jump at this cut reduces at the increase of . In such a way, the function approaches its limiting value almost everywhere, although there is set of singularities at negative integer values of .

The solution follows asymptotic at large values of real part of the argument, exponentially approaching the limiting value. In particular, for , this maximum limiting value is in the left hand side of the figure, and to its minimum limiting value in the right hand side. For , these limiting values are and .

The trace of the solution along the real axis corresponds to the red dotted curve in Figure 1. Other solutions of the recursive equation , that may grow up along the real axis, can be constructed in a similar way, but they do not satisfy the criteria formulated in the definition of tetration; in particular, .

Case

Fig.5. Tetration at .

At , the limiting value , and, asymptotically,

The function is shown in figure 5.

Levels are shown with thick black lines.

Levels are shown with thick red lines.

Levels are shown with thick blue lines.

Intermediate levels are shown with thin lines.

There is a cut at , , but the hump of the function at the cut reduces as reduces, id est, with increasing . In such a way, everywhere, at , the function approaches its limiting value . almost everywhere, although there is a set of singularities at negative integer values of .

Behavior of this function at real values of argument is shown in Figure 1 with thin solid line. Other solutions of the recursive equation can be constructed in the similar way; they may grow up along the real axis, but they do not satisfy criteria formulated in the definition of tetration; in particular, .

Case

(CC) Image: Dmitrii Kouznetsov
Fig.6. Tetration at base and .

At , tetration is asymptotically periodic. It decays exponentially to the fixed points and in the upper and lower halfplane. This allows to express it through its values along the imaginary axis, using the Cauchi integral. [1].

For the base and tetration is plotted in figure 6. Levels and are drawn with thick lines. The function has a logarithmic singularity at point -2 and cut at real values of the argument, smaller than -2. In the right hand side, symbols mean huge values that cannot be stored in the conventional floating point representation (logarithm, mantissa). In the upper left and lower left part of each graphic, the function approaches its asymptotic values and . Function is quasi-periodic; the same fractal structure reproduces again and again at the translation of argument with quasiperiod in the upper halfplane and in the lower halfplane.

There is cut at , . The jump of the function at this cut approaches almost everywhere, although there is set of singularities at negative integer values of .

Along the real axis, tetration for these values of the base is plotted also in figure 1 with thick solid and dashed lines.

Behavior along the real axis

(CC) Image: Dmitrii Kouznetsov
Derivatives of tetration at .

The growth of tetration along the real axis is crucially determined by its base. The graphic of this function is shown at the top of the article for . For , the derivatives of tetration are plotted in figure at the right.

is plotted with red;
is plotted with green;
is plotted with blue;
is plotted with pink.

Tetration is strictly growing function; its first derivative is positive. For the minimum of the derivative takes place in vicinity of and is slightly smaller than unity. At , the growth is limited by the minimum of the limiting values . Tetration approaches this limiting value exponentially. In particular, at , the limiting value is 2.

At , the growth is limited by the fixed point . Tetration approaches this limiting value as rational function.


The growing-up holomorphic solution of equation can be constructed in the similar way also at , but at the real axis, such a solution remains larger than unity, and the condition cannot be satisfied. Therefore, this solution cannot be interpreted as extension of iterative exponentiation for non-integer number of exponentiations; in this sense, such a solution is not a tetration.

Tetration shows explosively-fast growth along the real axis only at values .

Tetration at the base

Template:CC-Image
Tetration at the base

There is hypothesis that at base , the graphic of tetration is symmetric with respect to the line . This line, together with the graphic, is shown in Figure at let. The asymptotics and are also shown. The graphic looks symmetric with respect to , but no proof of this hypothesis was suggested. Now this hypothesis is negated; the apparent symmetry is only approximation. However, for , at the segment , the relation holds with 4 decimal digits [8][9].

Tetration at base e

Holomorphic tetration on the natural base is the most developed, at least up to the year 2008. In the rest of this article, it is assumed that although the majority of results allow a straightforward extension to the case of real .

Tetration has real fixed point , id est, solution of the equation . Its approximation

can be found, iterating equation , where slog is inverse function of tetration.

Pentation

Complex map of natural pentation, by [10]

Superfunction of tetration can be called pentation; it is solution of equations

Expansion of tetration in vicinity of its fixed point from the interval (for real base ) allows to define unique (and the only "true") pentation. Other superfunctions of tetration can be obtained by the periodic modification of the argument of pentation.

For base , the fixed point of tetration . and pentation to this base approaches this value at large negative values of the real part of the argument [11][10]. Complex map of natural pentation is shown in figure at left; .

Inverse of tetration

Fig.7. slog in the complex plane.

The inverse of tetration can be performed using the Newton method, solving equation , leading to

.

The inverse function has branchpoints and . For the kslog, at base , shown in the figure 7, the cuts are placed horizontaly, along the lines

, .

Due to the symmetry , it is sufficient to plot only half of the complex plane.

Fig.8. Shaded region is image of the upper halfplane at the mapping with slog function.

The mapping with function kslog is shown in figure 8. For in the shaded region, the relation

takes place. The upper part of the complex plane is mapped into the upper halfplane, and the lower halfplane is mapped to the lower halfplane. The real axis is mapped into the halfline . The fixed points is mapped to imaginary infinity, following the shaded strip.

In figure 8, the images of the grid lines and images of the grid lines are shown. These curves reproduce levels and shown in figure 6 for .

Together, the pair of functions tet and kslog allow to evaluate any iteration (including negative, fractional and even complex) of the exponential function.

Beyond the cutlines of slog

(CC) Image: Dmitrii Kouznetsov
3 options to put cutlines of the superlogarithmic function

The function slog has branchpoints, and there are many options to pit the cutlines. The cutlines, parallel to the abscias axis, considered above is one of many options. One can put the cutlines along the level . This cases simplify the plotting of the slog function, while the only tetration is available for the efficient evaluation. On the other hand, there is no simple expression for the parametrization of such a cutline, and the calculation of the position of the cutlines slows down the algorithm of evaluation of slog. In addition, there is tradition, that cutlines of the special functions are placed parallel to the real or imaginary axis. One could place the cutline from the branchpoints to the abscissa axis. In this case, the additional cutlines should go along the negative part of the real axis, which is not convenient for the applications. (For example, for the evaluation of the generalized exponential at non-integer , it is important, that the real axis belongs to the range of holomorphism of slog. There are 3 other options to place the cutlines parallel to the coordinate axes, keeping condition . There 3 options are shown in the figure.

The shaded region shows the domain

For evaluation of any of these superlogarithms, it is sufficient to have the efficient algorithm for the evaluation in this domain.

The left hand side picture corresponds to the slog, sa it is defined above. Within the strip between the cutlines, the function approaches its limiting value ; it becomes infinite in vicinity of the branch points and varies very slow in the rest of the complex plane.

The central part of the picture represents slog with vertical cutlines, let us call this function slogv. In the left hand side of the complex halfplane, the evaluation of slog begins with exponentiation of its argument. The exponential is periodic function; therefore, in this part, slog is periodic:

The third option, let us denote it slogr, is to put the cutlines in the direction of the real axis. Along the real axis, the function remains holomorphic, but there are additional branchpoints in the ranges . In these regions, the function slogr show the fractal behavior.

The slog, defined at the beginning, seems to be simpler than slogv and slogr; therefore, namely slog with horizontal cutlines, parallel to the real axis, going in the direction, opposite to the real axis, is used below.

Polynomial approximation

The Taylor series for the tetration can be written in the usual form:

where the th coefficient

is expressed through the derivative of the function. The coefficients of the expansion can be calculated using the straightforward differentiation of the representation through the Cauchy integral.

Taylor expansion at zero

(CC) Image: Dmitrii Kouznetsov
Fig.N. Approximation of tetration with polynomial of 25th power

For , the calculation gives the following values

0 = 1
1 ≈ 1.091767351258322138
2 ≈ 0.271483212901696469
3 ≈ 0.212453248176258214
4 ≈ 0.069540376139988952
5 ≈ 0.044291952090474256
6 ≈ 0.014736742096390039
7 ≈ 0.008668781817225539
8 ≈ 0.002796479398385586
9 ≈ 0.001610631290584341
10≈ 0.000489927231484419
11≈ 0.000288181071154065
12≈ 0.000080094612538551
13≈ 0.000050291141793809
14≈ 0.000012183790344901
15≈ 0.000008665533667382
16≈ 0.000001687782319318
17≈ 0.000001493253248573
18≈ 0.000000198760764204
19≈ 0.000000260867356004
20≈ 0.000000014709954143
21≈ 0.000000046834497327
22≈-0.000000001549241666
23≈ 0.000000008741510781
24≈-0.000000001125787310
25≈ 0.000000001707959267

The truncated Taylor series gives the polynomial approximation. In the upper right hand side of the Figure N, the polynomial

is shown in the complex plane.

Levels are shown with thick black curves.
Levels are shown with thin red curves.
Levels are shown with thin thin blue curves.
Levels are shown with thick red curves.
Levels are shown with thick blue curves.
Levels are shown with thick pink curves.
Levels are shown with thin green curves.

In the upper left corner of figure N, the same is shown for function

At the bottom left, the overlap of the upper two images is shown.

At the bottom right, lines of constant modulus and constant phase of holomorphic tetration in the same range.

The good approximation of tetration takes place in the range of order of unity or smaller; the radius of convergence of the series is 2.

Expansion at 3i

Taylor expansion of tetration developed at 3i, truncated at 30th power, plotted in the complex plane.

Coefficients of the expansion

can be evaluated in the similar way:

The plot of approximation of tetration with the resulting polynomial of 30th power is shown in figure. This approximation can be used for plotting of camera-ready pictures of tetration, using it and its conjugation at . With 50 terms, at , such approximation returns 14 significant figures.

Asymptotic expansion at large values of the imaginary part of the argument

(CC) Image: Dmitrii Kouznetsov
Fig.As. Deviation of tetration from its asymptotic expansion.

Here only the case of base is considered; although the generalization to the case is straightforward. In this section, the behavior of tetration is considered for moderate values of and large values of

At , in the upper halfplane, tetrations approach the fixed point of the logarithm. This approach is exponential. Using the exponential

as a small parameter, for some complex constant , the tetration can be estimated as follows:

Substitution of this expression into the equation gives

Value can be found fitting the tetration with function at large values of . Then, the expansion can be written in the form

assuming that and . At given , such a representation indicates the possible approximations of tetration. The deviations of tetration from these approximations are shown in figure Fig.As.

The four plots in fig.As correspond to the four asymptotic approximations. The deviations

,
,
,

are shown in the complex plane with lines of constant phase and constant modulus.

Levels are shown with thick red lines.
Level is shown with thick black lines.
Levels are shown with thick blue lines.
Levels are shown with scratched lines. (these lines reveal the step of sampling used by the plotter).
Levels are shown with thin red lines.
Levels are shown with thin blue lines.
Levels are shown with thin thick black lines.
Level is shown with thick red line.
Levels are shown with thick black lines.
Level is shown with thick red line.
Levels are shown with thick black lines.
Level is shown with thick green line.
Level is shown with thick black line.

The plotter tried to draw also

Level with thick black line and
Level with thin dark green line, which are seen a the upper left hand side corners of the two last pictures, but the precision of evaluation of tetration is not sufficient to plot the smooth lines; for the same reason, the curve for
in the last picture, in the upper right band side looks a little bit irreguler; also, the pattenn in the upper left corner of the last two pictures looks chaotic; the plotter cannot distinguish the function from its asymptotic approximation.

The figure indicates that, at , , the asymptotic approximation

gives at least 14 correct significant figures. At large values of the imaginary part, this approximation is more precise than the evaluation of tetration through the contour integral.

Approximation of tetration with elementary functions

(CC) Image: Dmitrii Kouznetsov
in the complex -plane.

Due to recurrent relation , it is sufficient to approximate tetration in any vertical strip of unit with in the complex plane. Some of such approximations are suggested in [1]. In principle, the numerical approximation of tetration with implementation of the Cauchi integral with finite sums [1] also should be considered as approximation with elementary function. However for the practical evaluation of tetration, shorter expressions are more suitable. One of such approximation comes from the Taylor expansion of function . The subtraction of logarithm remove the closest singularity that limits the radius of convergence of the Taylor series, and makes precise the approximation with finite sum. One of such approximation with one logarithm and polynomial of 100th power

is shown at the figure in the complex -plane. In vicinity of the origin of coordinates, say, , the most of terms are negligibly small, and the shortened sum still approximates the tetration. The first coefficients in this expansion are

More coefficients can be extracted from the generator of the figure. While , the approximation with 101 terms returns at least 14 correct significant figures.

Approximation of slog

Function slog, which is inverse of tetration, allows the approximation with elementary functions.

Taylor expansion

(CC) Image: Dmitrii Kouznetsov
Approximation of slog with polynomial of 16th power from the Taylor expansion at unity.

The Taylor series for the tetration can be inverted, gaining the expansion of the superlogarithm:

.

Approximations for the first 16 coefficients:

The partial sum with 16 terms (from zero to 16) is plotted in the figure in the complex plane. Lines of constant real part and constant imaginary part are drawn.

The radius of convergence of this series is determined by the distance to the closest singularity; the representation of with the Taylor series is valid for

Obviously, it fails at . For this case, the asymptotic expansion can be used.

Expansion at fixed point of logarithm

(CC) Image: Dmitrii Kouznetsov
Fitting of slog with the expansion around .

Superlogarithm can be approximated with expansion [12]

,

where is fixed point of logarithm. This expansion indicates the ways to construct the asymptotic approximations of slog. The coefficients can be expressed from the asymptotical analysis of equation . Also, they can be expressed from the asymptotical estimate of tetration at large values of the imaginary part of the argument. The evaluation of first coefficients gives

These coefficients allow to approximate slog in vicinity of the fixed point of logarithm with function

.

In the figure this is shown in the complex -plane.

Levels are shown with thick black lines.
Levels are shown with thin red lines.
Levels are shown with thin dark green lines.
Levels is shown with thick green line.

(Deviation of this line from the real axis indicates the error of the approximation.)

Level is shown with thick red line
Levels are shown with thin dark blue lines.

For comparison, dashed lines show the precise evaluation for some of the levels above for the robust implementation of the slog function as inverse of tetration. While , the deviation of these dashed lines from the levels for function is not seen even at the strong zooming-in of the central part of the figure.

Approximation of slog with elementary functions

(CC) Image: Dmitrii Kouznetsov
Approximation of slog with function fslog.
(CC) Image: Dmitrii Kouznetsov
Numerical test of approximation of slog with function fslog.

The precision of approximation of slog (with fixed precision of the arithmetics used for the tvaluation) can be extended, takung unto account the singularitues of slog at the fixed points. From the asumptotical representation above, one can guess the robust representation for slog:

The coefficients of this expansion are real. The first coefficients:

This representation allows construction of approximations, truncating the series. One of such approximations

is shown in figure in the complex plane. Lines of constant and those of constant are plotted.

The range of approximation of slog with this function is wider than that with the Taylor expansion at unity. The extended range of approximation allows its validation with the numerical test of identities

The residuals

and

are shown in the figure with levels const. In the voided regions in vicinity of and , the residual is at the level of . (It is difficult to make the residual smaller, using the arithmetics with double complex variables.) This test indicates, that at , the approximation of slog with two the logarithms and the polynomial of 50th power gives at least 9 correct significant figures.

Iterated exponential and

(CC) Image: Dmitrii Kouznetsov
Fig.9. in the complex plane for various .

Pre-historic approach to the iterated exponential

Especially interesting is the case of iteration of natural exponential, id est, . Existence of the fractional iteration, and, in particular, existence of operation was demonstrated in 1950 by H.Kneser. [12]. However, that time, there was no computer facility for the evaluation of such an exotic function that ; perhaps, just absence of an apropriate plotter did not allow Kneser to plot the distribution of fractal exponential function in the complex plane for various values of , shown in Fig.9.

(CC) Image: Dmitrii Kouznetsov
Fig.10. versus

The Implementation through the tetration

Holomorphic tetration allows to extend the iterated exponential

For non-integer values of . It can be defined as

If in the notation the superscript is omitted, it is assumed to be unity; for example . If the subscript is omitted, it is assumed to be , id est,

Iterated exponential in the complex plane

Function is shown in figure 9 with levels of constant real part and levels of constant imaginary part. Levels and are drown with thick lines. Red corresponds to a negative value of the real or the imaginaryt part, black corresponds to zero, and blue corresponds to the positeive values. Levels are shown with thin red lines. Levels are shown with thin green lines. Levels and Levels are marked with thick green lines, where is fixed point of logarithm. At non-integer values of , and are branch points of function ; in figure, the cut is placed parallel to the real axis. At there is an additional cut which goes along the negative part of the real axis. In the figure, the cuts are marked with pink lines.

Iterated exponential of a real argument

For real values of the argument, function is ploted in figure 10 versus for values .
in programming languages, inverse function of exp is called log.

For logarithm on base e, notation ln is also used. In particular, , and so on.

At least for the real and big enough , the relation holds, which is qute analogous of relations and . However, at negative or negative , value should be big enough, that and and are defined, see figure 10. For example, at , we need . In particular, at , , we have .

Application of iterated exponential

The iterated exponential, that can be implemented with holomorphic tetration, may have various applications. In particular, The at could describe a process that grows faster than any polynomial, but slower than any exponential. In such a way, the iterated exponential, at the proper implementation, should greatly extend the abilities of fast and precise fitting of functions. This is just analogy of function which, at fractal values of , may be good for description of a function that grows faster than any linear function but slower than any quadratic function.

Similar functions

Withdrawal of some of requirements from the definition of tetration allows the huge variety of similar functions.

Entire solutions of

Withdrawal of the requirement and allows the solution by Kneser [12], which is entire and also could be used to build up various powers of the exponential; in particular, . Such entire function is shown in upper part of figure 1c in [1], in order to reveal the asymptotic behavior of holomotphic tetration.

Withdrawal of condition allows to construct solutions, similar to the growing tetration, for base . Although such solutions cannot be interpreted as generalization of exponential iterated times, they can be useful for generalization of exponential function.

Non-holomorphic modification of tetration

(CC) Image: Dmitrii Kouznetsov
Fig.11. Almost identical function in the complex plane.
(CC) Image: Dmitrii Kouznetsov
Fig.12. Motified tetration at the complex plane.
(CC) Image: Dmitrii Kouznetsov
Fig.13. Zoom in of fig.12

Withdrawal of requirement of holomorphicity from the definition of tetration allows functions, which look like the tetration, at least along the real axis. Even the reduction of the range of holomorphism in the requirement allows to consider tetration with modified argument. The modified tetration can be defined as

,

where , and is a 1-periodic function. The simple example of such function is

In this case, along the real axis, the function is almost identical to its argument; and values of the modified tetration are close to values of tetration. Being plotted at figure 1 or in figure 10, the deviation of such function from the identity is small, and the deviation of the modified tetration from tetration is also small. If the figures are printed in the real scale, then the deviation of the curves would be of order of atomic size.

However the difference becomes seen at the complex values of the argument. In figure 11, function is plotted in the complex plane. Levels of constant real part and those of constant imaginary part are drawn. In vicinity of the real axis, these lines almost coincide with the gridline; the grid is drawn with step unity and extended one step to the right and one step to the left from the graphic. In order to show that it behaves as it if would be a continuation of the plot. At , the deviation from the identical function becomes visible, and at , the has many points with real values, including those with various negative integer values. The tetration has cut at negative values of and singularities there. Therefore, the cuts of the modified tetration are determined by the lines , and modified tetration unavoidable has singularities in points such that . These singularities are determined by the function and do not depend on the base of tetration. In figure 11, the lines are seen not only along the real axis, but also at the top and at the bottom of the figure. In such a way, figure 11 shows the cutlines of the modified tetration. One has no need to evaluate tetration in order to see the margin of the change of holomorphism of the modified tetration.

In figure 12, the modified tetration is plotted in the complex -plane. The additional cuts are seen in the upper and lower parts of the figure 12. Only within the strip along the real axis, the function is holomorphic. While the amplitude of sinusoidal is of order of , the strip of holomorphism is wider than unity, although this width slightly reduces along the real axis.

In order to see the behavior of the modified tetration in vicinity of the additional singularities, in fig.13, the zooming-in of the part of figure 12 is shown. The zoom has improved resolution, so, in its turn, it can be zoomed in to the size of the screen of a computer to see the details. In each cell of the grid, the small and deformed image of the central part of the fig.12 appears.

One has no need to evaluate tetration in order to reveal its singularities outside the real axis. All the solutions of equation for integer are singularities (branchpoints) of the modified tetration.

Let .

The following theorem is suggested: For any entire 1-periodic function such that , is not a constant, , there exist . sequence Although this theorem is not yet proved, the intents to construct at least one example of function, that would contradict it, were not successful. This theorem is somehow independent from the theory of tetration, but it indicates, that any modified tetration cannot be holomorphic in the range .

According to the theorem above, the modified tetration does not satisfy the condition of quasi-periodicity, and does not satisfy the criterion of holomorphism in the definition of tetration. The sequence of cutlines for the specific example of modified tetration is seen in figure 13; the modified tetration is not even continuous. An addition to function some highest sinusoidals brings the discontinuities even closer to the real axis. This indicates, that if in at least one point at the real axis between and , some solution of equations , differs from tetration tet for at least , then this solution is not holomorphic in the range .

Even small deformation of tetration tet breaks its continuity. Similar reasons in favor of uniqueness tetration are suggested also in [1]. There is only one tetration, that satisfies requirements of the definition, although the rigorous mathematical proof of the uniqueness is still under development.


References

  1. 1.0 1.1 1.2 1.3 1.4 1.5 D.Kouznetsov. (2009). "Solutions of in the complex plane.". Mathematics of Computation, 78: 1647-1670. DOI:10.1090/S0025-5718-09-02188-7. Research Blogging. Cite error: Invalid <ref> tag; name "k" defined multiple times with different content Cite error: Invalid <ref> tag; name "k" defined multiple times with different content Cite error: Invalid <ref> tag; name "k" defined multiple times with different content Cite error: Invalid <ref> tag; name "k" defined multiple times with different content Cite error: Invalid <ref> tag; name "k" defined multiple times with different content
  2. 2.0 2.1 D.Kouznetsov. Ackermann functions of complex argument. Preprint of the Institute for Laser Science, UEC, 2008. http://www.ils.uec.ac.jp/~dima/PAPERS/2008ackermann.pdf Cite error: Invalid <ref> tag; name "k2" defined multiple times with different content
  3. R.L.Goodstein (1947). "Transfinite ordinals in recursive number theory". Journal of Symbolic Logic 12.
  4. 4.0 4.1 P.Walker (1991). "Infinitely differentiable generalized logarithmic and exponential functions". Mathematics of computation 196: 723-733. Cite error: Invalid <ref> tag; name "w" defined multiple times with different content
  5. M.H.Hooshmand. (2006). "Ultra power and ultra exponential functions". Integral Transforms and Special Functions 17 (8): 549-558.
  6. N.Bromer. Superexponentiation. Mathematics Magazine, 60 No. 3 (1987), 169-174
  7. W.Ackermann. ”Zum Hilbertschen Aufbau der reellen Zahlen”. Mathematische Annalen 99(1928), 118-133.
  8. D.Kouznetsov, H.Trappmann. Portrait of the four regular super-exponentials to base sqrt(2). Mathematics of Computation, 2010, v.79, p.1727-1756. http://www.ams.org/journals/mcom/2010-79-271/S0025-5718-10-02342-2/home.html http://mizugadro.mydns.jp/PAPERS/2010q2.pdf
  9. https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0 http://www.ils.uec.ac.jp/~dima/BOOK/202.pdf http://mizugadro.mydns.jp/BOOK/202.pdf Д.Кузнецов. Суперфункции. Lambert Academic Press, 2014. (In Russian, 328 pages.)
  10. 10.0 10.1 https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0 http://mizugadro.mydns.jp/BOOK/202.pdf Д.Кузнецов. Суперфункции. Lambert Academic Publishing, 2014. (In Russian)
  11. D.Kouznetsov. Holomorphic ackermanns. http://mizugadro.mydns.jp/PAPERS/2014acker.pdf under consideration
  12. 12.0 12.1 12.2 H.Kneser. “Reelle analytische L¨osungen der Gleichung und verwandter Funktionalgleichungen”. Journal fur die reine und angewandte Mathematik, 187 (1950), 56-67. Cite error: Invalid <ref> tag; name "kneser" defined multiple times with different content Cite error: Invalid <ref> tag; name "kneser" defined multiple times with different content