User:Dmitrii Kouznetsov/Analytic Tetration

Fig.0. Graphic of analytic tetration $f=F_{2}(x)=\exp _{2}^{x}(1)$ at versus $x$.

Abstract

Analytic tetration is defined as mathematical function that coincides witht the tetration at integer values of the argument and is analytic outside the negative part of the real axis. Existence of such a function is postulated; and arguments in favor of uniqueness of such a function are considered. The algorithm of evaluation is suggested. Examples of evaluation, pictures and tables are supplied. The application and the generalization is discussed.

Preface

The colleagues indicated so many misprints in my papers about tetration, posted at my homepage , that I want to give them opportunity to correct them in real time.

Especially I invite

• Arthur Knoebel
• Henryk Trappmann
• Andrew Robbins

to edit this file.

I consider the topic very important and urgent. The analytic tetration should be investigated and discussed right now; overvice, the non-analytic extension may become an ugly standard in mathematics of computation; the implementation of hige numbers with non-analytic tetration would make difficult realization of arithmetic operations and cause a lot of incompatibilities.

This is my apology for posting this research now, while the rigorous proof of existence and uniqueness of the analytic tetration is not yet found. My believe is based on the numerical check of the hypothesis of the existence and uniqueness, on smallness of the residual at the substitution of the function to the tetration equation and beauty of the resulting pictures. I cannot imagine that the agreement with 14 decimal digits occurs just by occasion without deep mathematical meaning.

In such a way I apologize for postulating of statements which should be prooven by the rigorous mathematical deduction.

Introduction

Quick start

Roughly, super-exponential

(1) $F_{b}(z)=\exp _{b}^{z}(1)$

is combination of $x$ exponentials on base $b$. For example,

$~F_{b}(0)=\exp _{b}^{0}(1)=1$

$~F_{b}(1)=\exp _{b}(1)=b$

$~F_{b}(2)=\exp _{b}^{2}(1)=\exp _{b}(\exp _{b}(1))=b^{b}$

$~F_{b}(3)=\exp _{b}^{3}(1)=\exp _{b}{\Big (}\exp _{b}{\big (}\exp _{b}(1){\big )}{\Big )}=b^{b^{b}}$

and so on. However, such definition is good only for positive integer values of $x$. In general, the superexponential can be defined through the Abel equation

(2) $~\exp _{b}{\Big (}F_{b}(z){\Big )}=F_{b}(z+1)$

with additional condition that

(3) $~F_{b}(0)=1$

Then, at least for positive values of $~b$ and positive integer values of $~z$, such a definition can be used for the evaluation of tetration.

In this paper, the way to define tetration for non-integer argument is described. For real values of the argument, at $~b=2$, such a tetration is plotted on figure 0. In the following sections, I describe, why is it so important, how to define the tetration for non-integer values of the argument, how can it be evaluated with high precision and why it is the only correct way to define analytic tetration.

Additional argument

One can consider to add the additional argument, replacing $~\exp _{b}(1)$ to $~\exp _{b}(a)$. This may have sense, while $~z$ is allowed to have only integer values. However, at the implementation of "good" tetration, the "argument" $~a$ can be considered as inverse superexponential of some argument, $~a=F_{b}(\alpha )$; then, $~\exp _{b}^{z}(a)=F_{b}{\Big (}F_{b}(\alpha ){\Big )}$; in the way, similar to that of convential logarithms: it is sufficient to investigate properties of natural logarithm ln; then, any other can be expressed as $\log _{b}(z)=\ln(z)/\ln(b)$.

$\exp _{b}^{n}(\exp _{b}^{m}(z))=\exp _{b}^{m+n}(z)$

The exponentiation of tetration is equivalent to increment of its argument. While summaton operation forms the group, exponentiation does too.

$\exp _{b}^{m}(F_{b}(z))=F_{b}(z+m)$

Inverse function and group properties

In this section, I write $~F$ instead of $~F_{b}$ and $~\exp(z)$ instead of $~\exp _{b}(z)$; onitting indices. However, you may recover them at any moment.

(I am not sure which notation is best. D.)

The speculation of the previous subseciton can be written shorter.

Assume there exist function $~F^{-1}$ such that $~F{\Big (}F^{-1}(z){\Big )}=z$.

Let $~s=\exp ^{u}(t)=F{\Big (}u+F^{-1}(t){\big )}$ and $~q=\exp ^{v}(s)=F{\Big (}v+F^{-1}(s){\big )}$ .

Then $~q=\exp ^{v}{\Big (}\exp ^{u}(t){\Big )}=F{\Big (}v+F^{-1}{\Big (}F(u+F^{-1}(t)){\Big )}{\Big )}=F{\Big (}v+u+F^{-1}(t){\Big )}$.

You can put subscript $~_{b}$ to $~F$ and $~\exp$ in the defuction above, and it will be seen, that we have no need to deal with funciton of 2 variables, considering $\exp _{b}^{z}(x)$; it can be expressed in terms of $F_{b}(z)$. However, we need to specify, what set shold be $~t$ and $~s$ from in the deduction above: must they be positive integer, or they can be real, of they can be also complex numbers.

History of tetration and huge numbers

Perhaps, every researcher used to see diagnostivs "floating overflow" at the evlauation of an expression with huge numbers....

Asymptotic

Assume, the tetration $F$ is defined with the Abel equation

(10) $\exp _{b}(F(z))=F(z+1)$

and assume the condition

(11) $F(0)=1$

Asymptotic

The analytic extension of tetration $~F(z)$ is supposed to grow fast along the real axis of the complex $z$-plane, at least for some values of base $b$. However, it has no need to grow infinitely in the direction of imaginary axis. For mathematics of computation, it would be better, if $F(x+{\rm {i}}y)$ remains bounded at $y\rightarrow +\infty$. Consider this possibility.

Assume, there exist solution $F(z)$ of equations (10), (11), analytic in the $\Re (z)\geq 0$, with, probably, countable number of cuts and singularities at $\Re (z)<0$, with exponential asymptotic behavior

(12) $F(z)=L+\varepsilon (z)+o{\big (}\varepsilon (z)^{2}{\big )}$

within some range, where

(13) $\varepsilon (z)=\exp(Qz+r)$,

$Q$ and $r$ are fixed complex numbers, and $L$ is eigenvalue of logarithm, solution of equation

(14) $x=\log _{b}(x)$.

FIg.1. Example of graphic solution of equation $x=\log _{b}(x)$ for $b={\sqrt {2}}$ (two real solutions, $x=2$ and $x=4$), $b=\exp(1/{\rm {e)}}$ (one real solution $x={\rm {e}}$) $b=2$ (no real solutions).

Three examples of graphical solution of equation (14) are shown in figure 1 for $b={\sqrt {2}}$, $b=\exp(1/{\rm {e)}}$, and $b=2$.

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

At $b={\sqrt {2}}$, there exist 2 solutions, $x=2$ and $x=4$.

At $b=\exp(1/{\rm {e)}}$ there esist one solution $x={\rm {e}}$.

and $b=2$, there are no real solutions.

In general,

• at $b<\exp(1/{\rm {e)}}$ there are two real solutions
• at $b=\exp(1/{\rm {e)}}$, there is one soluition, and
• at $b>\exp(1/{\rm {e)}}$ there esist two solutions, but they are complex.

In particular, at $b={\sqrt {2}}$, the solutions are
$x=L_{{\sqrt {2}},1}=2$ and $x=L_{{\sqrt {2}},2}=4$
.

At $b=2$, the solutions are
$x=L_{2}\approx 0.824678546142074222314065+1.56743212384964786105857\!~{\rm {i}}$ and
$x=L_{2}^{*}\approx 0.824678546142074222314065-1.56743212384964786105857\!~{\rm {i}}$.

At $b={\rm {e}}$, the solutions are
$x=L_{\rm {e}}\approx 0.318131505204764135312654+1.33723570143068940890116\!~{\rm {i}}$ and
$x=L_{\rm {e}}^{*}\approx 0.318131505204764135312654+1.33723570143068940890116\!~{\rm {i}}$.

Few hundred straightforward iterations of equation (14) are sufficient to get the error smaller than the last decimal digit in the approximations above.

The solutions $x=L_{1}$ and $x=L_{2}$ of equation (14) are plotted in figure 2 versus $\beta =\ln(b)$ with thin black lines. Let $L_{1}<L_{2}$, and only at $\ln(b)=1/e$, the equality $L_{1}=L_{2}$ takes place.

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

The thin black solid curve at $\beta \geq 1/{\rm {e}}$ represents the real part of the solutions $L$ and $L^{*}$ of (14); the thin black dashed curve represents the two options for the imaginary part; the two solutions are complex conjugaitons of each other. Let $\Re (L)>0$.

Increment and periodicity

Parameter $Q$ in equaiton (13) has sense of asymptotic increment. Consider possible values of $Q$.

Substitution of expression (13) into the eq. (12) gives
$L+\varepsilon (z){\rm {e}}^{Q}+o(\varepsilon (z))=\exp \!{\Big (}\ln(b)(L+\varepsilon (z)+o(\varepsilon (z)){\Big )}$ Using equation (14) gives

$L+\varepsilon (z){\rm {e}}^{Q}+o(\varepsilon (z))=L~{\Big (}1+\ln(b)\varepsilon (z)+o(\varepsilon (z)){\Big )}$

Then

(16) $Q=\ln(\ln(b))+\ln(L)$.

At $\ln(b)<1/{\rm {e}}$, the two real increments correspond the the two real eigenvalues of logarithm. Let $Q_{1}=\ln(\ln(b))+\ln(L_{1})$ and $Q_{2}=\ln(\ln(b))+\ln(L_{2})$. Both $Q_{1}$ and $Q_{2}$ are plotted in Figure 2 with thick solid lines; $Q_{1}\leq 0$ and $Q_{2}\geq 0$, and only at $\ln(b)=1/{\rm {e}}$ the equality $Q_{1}=Q_{2}=0$ holds.

In such a way, at $\ln(b)<1/{\rm {e}}$, there exist positive increment, which corresponds to $L_{2}$ and means, that the perturbation $\epsilon (z)$ of the stationary solution grows up in the diredtion of real axis; there exist negative increment (decrement, if you like), which corresponds to $L_{1}$ and indicates, that the perturbation grows up in the opposite direction. At $0<\ln(b)<2/{\rm {e}}$, there may exist tetration $F(z)$ with two different asumptotics. At large positive values of $\Re (z)$ it decays to the $L_{1}$ and At large negative values of $\Re (z)$ it decays to the $L_{2}$.

At $\ln(b)\geq 1/{\rm {e}}$, there exist two possible increments, $Q$ and $Q^{*}$. The real part of $Q$ is plotted with thick solid line, and the imaginary part is shown with thick dashed line. The real part is positive; id est, the perturbation $\varepsilon (z)$ grows up in the direction of real axis.

For computational mathematics, it would be good ot construct the solution, that decays both in the direction of the imaginary asis and in the opposite direction. Such a solution should have asymptotics (12),(13) at large positive values of the imaginary part of the argument and its conjugation at large negative imaginary part of the argument.

The exponential asymptotics implies the asumptotic perioditity. The period

$T={\frac {2\pi {\rm {i}}}{Q}}$

is shown in FIgure 2 with dotted lines; for the case when period is negative, the minus period is plotted. It is amaising, that at $\ln(b)>1/{\rm {e}}$, the imajinary part is slightly larger than unity and remains almost constant.

In particular, at $b={\sqrt {2}}$, the periods are

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_{\rm b,1} = \frac{2\pi \rmi}{Q_{b,1}}\approx -17.143148179354847104 \!~\rm i] }
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_{\rm b,2} = \frac{2\pi \rmi}{Q_{b,2}}\approx 19.236149042042854710 \!~\rm i] }

At $b=2$, the asymptotic period
$T=5.58414243554338946020010+1.05421836033693734654000\!~{\rm {i}}$;
and at $b={\rm {e}}$, the asymptotic period
$T=4.44695072006700782711227+1.05793999115693918376341\!~{\rm {i}}$.

Uniqueness

While one deal with solutions of system (10), (11) at the real axis, one imagine some real-analytic extension shown in Figure 0 and consider also

(100) $G(z)=F(z+k(z))$

where

(101) $k(z)=\sum _{n=-\infty }^{\infty }\alpha _{n}\exp(2\pi {\rm {i}}nz)$

Such a funciton is also soluiton of the Abel equation; at

(102) $\alpha _{-n}=(\alpha _{n})^{*}$

and

(103) $\alpha _{0}=-\sum _{n\neq 0}\alpha _{n}$

function $~G$ is real and passes through the same points as </math>~F</math> at integer values of the argument.

and at small values of coefficients $\alpha$, function $~G$ looks smooth, and it is difficult to guess, which of them is "true". For this reason, for the standard mathematical representation, the non-analytic stepwice funciton uxp was suggested .

However, the difference between functions $~F$ and $~G$ becomes seen, if one of them, for example, $~F$, is analytic and regular in some wide range, function $~G(z)$ will be analytic only within the strip $~|z|<y_{0}$; order of magnitude of $~y_{0}$ can be estimated with

$y_{0}\approx {\rm {MIN}}_{n}~{\frac {\ln(2/|\alpha _{n}|)}{2\pi n}}$

At larger values of the imaginary part of the argument, the periodic function $k$ takes huge values, including various negative integers. Namely at these values, function $~F$ has singulatities.

Various functions may satisfy equations (2),(3). For an abstract excersise in comlex functional analysis, all animals are created equal. However, for the applicaitons in the computational mathematics, some of them are more equal than other . As such a more equal animal we should choose the solution with simplest behavior, with minimum of singulatities, and easiest for the evaluation.

At small base, both values of quasiperiod are pure imaginary. This periodic or quasi-periodic behavior is similar to that of the conventional exponential, but function $~F$ does not grow to infinity in the direction of the real axis, aproaching eigenvalue $L$ of logarithm. One example of such a behavior $~b={\sqrt {2}}$ is shown in Figure 3 for $~b={\sqrt {2}}$. The following section describes, how it was ploted.

At $~b>\exp(1/{\rm {e}})$ values of period $T$ are complex; and $\Re (T)>0$. This corresponds to the exponential growth of the asymptotic solution in the direction of the real axis; at large positive values of $\Re (Qz)$ the asymptotic does not approximate the function which grows faster than any exponential.

Small base. Base $b={\sqrt {2}}$

Fig.3. Tetration $f=F(z)$ at the complex $z$-plane. Levels of integer $\Re (f)$ and levels of integer $\Im (f)$ are shown with thick lines. Grid covers with unit step range ($-10\leq \Re (z)\leq 10$, $-10\leq \Im (z)\leq 10$).

For evaluation of Tetration we need to assume that it exists, and has asymptotic (30).

Consider the case of $1<b<\exp(1/{\rm {e)}}$. at fixed values of real part of the argument, the function has periodic behavior in the direction of imaginary axis, as it is shown in Fig. 3. In this figure, the example with $b={\sqrt {2}}$ is used. Function is periodic; period $T\approx 17.1431~{\rm {i}}$.

The limiting values:

$\lim _{x\rightarrow \infty }F_{\sqrt {2}}(x+{\rm {i}}y)=2$

$\lim _{x\rightarrow -\infty }F_{\sqrt {2}}(x+{\rm {i}}y)=4$

At non-zero values of the imaginary part of the argument, the function decays to these asymptotiv values.

At the real axis, the $~F(z)$ has cut at $z\leq -2$.

Due to the periodicity, the function has cuts also at

$Z+Tn\leq -2$ for integer $~n$.

As all other tetrations, it has singularity at $z=-2$, $~F(-1)=0$ and $~F(0)=1$.

Other, more singular tetrations can be obtained by light periodic deformation of the argument.

Base 2. Ackermann function

FIg.4. Ackermann function $f=A(4,z)=F_{2}(z+3)-3$ at the complex $z$ plane. Grid covers the range ( $-10\leq \Re (z)\leq 10~,~-4\leq \Im (z)\leq 4$ with unit step. Levels of integer values of the real part and the same for imaginary part are shown with thick lines.

Tertration at base 2 can be expressed through the 4th Ackermann function $F_{2}(z)=A(4,z-3)+3$ This Ackermann function is plotted in fig.4. As base $b>\exp(1/{\rm {e)}}$, the function is not periodic. However, there is asymptotic periods $L$ and $L^{*}$ in the upper and lower half-planes. In vicinity of positive part of the real axis, the function shows rapid growth, and it is not possible to draw the levels there. All the singularities of the Ackermann function are at the negative integer values smaller than $-4$.

Base e. Natural tetration

Fig.5.Analytic tetration at the base e , $f=F_{\rm {e}}(z)$ in the complex $z$plane. The grid covers the range ($-9\leq \Re (z)\leq 9$, $-4\leq \Im (z)\leq 4$). Levels $|f|=\exp(0),\exp(\pm 1),\exp(\pm 2),..$ and levels ${\rm {Arg}}(f)=0,\pm 1,\pm 2,\pm 3$ are shown with thick lines.

Base $b-{\rm {e}}$ is the most natural choise. In this case, the increment $Q$ is equal to the asymtotic value $L$.

At the translation for inity along the real axis,

Superlogarithm

Fig.7. Inverse of the analytic tetration $f=F_{\rm {e}}^{-1}(z)$ at the complex $z$-plane.

Inverse of the natural tetration can be considered as superlogarithm. Equilines of funciton $f=F_{\rm {e}}^{-1}(z)$ are shown in Figure 7. Levels of integer real part and those of integer imaginary part are shown with thick lines. This function has two branchpoints at eigenvalues of logarithm. In the figure, the cuts are placed along the lines $\Re (f)=-2$; then, the function is regular in vicinity of the real axis, approaching the limitig value $-2$ at $-\infty$ and slowly growing up at positive values of the argument.

Discussion

Analytic tetrations are shown for base $b={\sqrt {2}}$, $b=2$, and $b={\rm {e}}$. In the similar way, the tetration on other bases (for example, $b=10$ ) can be plotted.

Knowledge of the asymptotic behavior gives the key to the efficient evaluation.

references

See also the discussion at http://math.eretrandre.org/tetrationforum/forumdisplay.php?fid=3