Tetration: Difference between revisions
imported>Dmitrii Kouznetsov (→Evaluation of tetration: make-uo formulas. Misprints.) |
mNo edit summary |
||
(160 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
{{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>.}} | ||
old version !--> | |||
( | '''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: | ||
<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. | |||
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. | |||
==Definition== | |||
== | 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 | ||
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 | |||
<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> . | ||
The definition above generalizes the definitions suggested | 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. | ||
<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]], | 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. | <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>. | 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 | 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 41: | Line 68: | ||
|volume=12 | |volume=12 | ||
|year=1947}} | |year=1947}} | ||
</ref>. This | </ref>. | ||
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. | |||
==Real values of the arguments== | 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 | <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 54: | Line 125: | ||
===Fast growth and application=== | ===Fast growth and application=== | ||
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"> | |||
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 | 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=== | ||
Line 74: | 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 | 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 84: | Line 155: | ||
\right. | \right. | ||
</math> | </math> | ||
The generalization of the 4th Ackermann | 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 | 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>. | 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 | 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 | 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|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 this equation are called [[fixed point]]s of logarithm. | ||
===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 | 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 | <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 | *at <math>b\!=\!\exp(1\!/\!\rm e)</math> , there is one solution, and | ||
*at <math>b>\exp(1/\rm e)</math> there | *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> | <math>L\!=\!\!L_{\sqrt{2},1}=2</math> and | ||
<math> | <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> | <math>L\!=\!L_2 \!\approx\! 0.824678546142074222314065+1.56743212384964786105857 \!~\rm i </math> and<br> | ||
<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> | <math>L=L_{\rm e} \approx 0.318131505204764135312654+1.33723570143068940890116 \!~\rm i</math> | ||
and<br> | and<br> | ||
<math> | <math>L=L_{\rm e}^* \approx 0.318131505204764135312654-1.33723570143068940890116 \!~\rm i</math>. | ||
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|right|700px|FIg.3. parameters of asymptotic of tetration versus logarithm of the base}} | |||
The solutions | The solutions | ||
<math> | <math>L=L_1 </math> and | ||
<math> | <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 168: | 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 | <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> | 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 | 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 < | 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 | 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== | ||
As the | 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 | 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 | <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>}} | |||
For <math> 1<b<\exp(1/\mathrm{e}) </math>, the period <math>T</math> is imaginary. | |||
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>. | ||
There is 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 | 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 | 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>b\!=\!\sqrt{2}</math>, these limiting values are | ||
<math>L\!=\!4</math> and | <math>L\!=\!4</math> and | ||
Line 211: | Line 284: | ||
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 solutions of the recursive equation <math>\exp_b\!\big(F(z)\big)\!=\!F(z\!+\!1)</math>, | 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 | 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 | ===Case <math> b=\exp(1/\mathrm{e}) </math>=== | ||
[[Image:HolomorphicTetrationBaseExp1e01.jpg|right|thumb| | [[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>, the | 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}( | : <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. | ||
Line 229: | 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>, but the hump of the function at the cut reduces as | 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>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>. | <math>|z|\rightarrow \infty</math>, the function approaches its limiting value <math>L=\mathrm{e}\approx 2.71</math>. | ||
almost everywhere, although there is set of singularities at | 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 | 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|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 | ||
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 | <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. | ||
There is cut at <math>x<-2</math>, <math>y=0</math>. The jump of the function at this cut approaches | 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> | <math>L\!-\!L^*</math> | ||
almost everywhere, although there is set of singularities at | almost everywhere, although there is set of singularities at negative integer values of <math>x<-2</math>. | ||
Along the real axis, tetration for | 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=== | ===Behavior along the real axis=== | ||
The growth of tetration along the real axis is | {{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. | |||
At <math> 1\! <\!b \!< \exp(1 \! / \mathrm{e} )</math> the growth is limited by the | 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. | <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 | 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. | <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 | 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 | 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 | <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. | ||
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 | Tetration <math>\mathrm{tet}_b</math> shows explosively-fast growth along the real axis | ||
<!-- and quasi-periodic behavior in the complex plane !--> | <!-- and quasi-periodic behavior in the complex plane !--> | ||
only at values <math>b \! >\! \exp(1 \! / \! \mathrm{e} )</math>. | 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. | [[Image:Kslog00.jpg|right|thumb|350px| | ||
The inverse of tetration <math>t=\mathrm{ | <!-- | ||
{{#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 | 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| | [[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{ | :<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 | <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 296: | 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 | 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|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=== | ||
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 | 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|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 325: | 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 337: | 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 354: | Line 808: | ||
In particular, at <math>u\!=\!1</math>, <math>v\!=\!-1</math>, we have | 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>. | <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== | ||
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]]. | |||
===Non-holomorphic modification of tetration=== | |||
{{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. | |||
However the difference becomes seen at the complex values of the argument. | |||
In figure 11, function <math>J(z)=z\!+\!h(z)</math> is plotted in the complex <math>z\!=\!x\!+\!\mathrm{i} y</math> 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 <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. | |||
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. | |||
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 <math>z</math> of equation <math>z+h(z)=n</math> for integer <math>n<-1</math> are singularities (branchpoints) of the modified tetration. | |||
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>. | |||
<math> | |||
<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== | ||
{{reflist}} | |||
[[Category:Suggestion Bot Tag]] | |||
Latest revision as of 06:00, 26 October 2024
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:
For real values of the argument and various values of the base , this is plotted in Fig.0.
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
- .
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
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
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 (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\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
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
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
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
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
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
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.
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
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
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
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
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
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
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
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
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
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.
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
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.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.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 - ↑ R.L.Goodstein (1947). "Transfinite ordinals in recursive number theory". Journal of Symbolic Logic 12.
- ↑ 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 - ↑ M.H.Hooshmand. (2006). "Ultra power and ultra exponential functions". Integral Transforms and Special Functions 17 (8): 549-558.
- ↑ N.Bromer. Superexponentiation. Mathematics Magazine, 60 No. 3 (1987), 169-174
- ↑ W.Ackermann. ”Zum Hilbertschen Aufbau der reellen Zahlen”. Mathematische Annalen 99(1928), 118-133.
- ↑ 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
- ↑ 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.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)
- ↑ D.Kouznetsov. Holomorphic ackermanns. http://mizugadro.mydns.jp/PAPERS/2014acker.pdf under consideration
- ↑ 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