Extension of tetration to real and complex heights

Abstract: The continuous tetrational function ${^x}r=\tau(r,x)$, the unique solution of equation $\tau(r,x)=r^{\tau(r,x-1)}$ and its differential equation $\tau'(r,x) =q \tau(r,x) \tau'(r,x-1)$, is given explicitly as ${^x}r=\exp_{r}^{\lfloor x \rfloor+1}[\{x\}]_q$, where $x$ is a real variable called height, $r$ is a real constant called base, $\{x\}=x-\lfloor x \rfloor$ is the sawtooth function, $\lfloor x \rfloor$ is the floor function of $x$, and $[\{x\}]_q=(q^{\{x\}}-1)/(q-1)$ is a q-analog of $\{x\}$ with $q=\ln r$, respectively. Though ${^x}r$ is continuous at every point in the real $r-x$ plane, extensions to complex heights and bases have limited domains. The base $r$ can be extended to the complex plane if and only if $x\in \mathbb{Z}$. On the other hand, the height $x$ can be extended to the complex plane at $\Re(x)\notin \mathbb{Z}$. Therefore $r$ and $x$ in ${^x}r$ cannot be complex values simultaneously. Tetrational laws are derived based on the explicit formula of ${^x}r$.

\section{Introduction}
Tetration, i.e., iterated exponentiation, is the fourth hyperoperation after addition, multiplication, and exponentiation \cite{Goodstein1947}. Tetration is defined as
\[^{n}r\colonequals \underset{n}{\underbrace{r^{r^{\udots^{r}}}}},\]
meaning that copies of $r$ are exponentiated \(n\) times in a right-to-left direction, or recursively defined as
\[^{n}r\colonequals
\begin{cases}
1 & (n=0)\\
r^{(^{n-1}r)} & (n \geq 1).
\end{cases}\]
Obviously \(^{n-1}r=\log_r (^{n}r)\). The parameter \(r\) and \(n\) are referred to as the base and height, respectively. In this paper, we use the notations to express tetration as
\[^{n}r=\tau(r,n)=\exp_{r}^{n}1,\]
and we write \(\tau(n)\) for simplicity if \(r\) is a constant.\vspace{0.2in}

\noindent
Tetration has been applied in fundamental physics \cite{Scot2006}, degree of connection in air transportation system \cite{Sun2017}, signal interpolation \cite{Khetkeeree2019}, data compression \cite{Furuya2019}, etc.
\vspace{0.2in}

\noindent
Tetration is known as a rapid growing function that may be useful for expressing superexponential growth or huge numbers, but it also shows saturation and oscillation.
\begin{figure}
\includegraphics[scale=0.43]{fig_discrete.pdf}
\centering
\caption{Behaviors of discrete tetration for different \(r\) (left), and cobweb plots (right).}
\label{fig:discrete tetration}
\end{figure}
As in Fig. \ref{fig:discrete tetration} the protean dynamics can be illustrated by cobweb plots. \vspace{0.2in}

\noindent
Euler \cite{Euler1783} proved that the limit of tetration \(^{n}r\), as \(n\) goes to infinity, converges if \(e^{(-e)}\leq r \leq e^{1/e}\). Eisenstein \cite{Eisenstein1844} gave the closed form of \(^{\infty}r\) for complex \(r\), and Corless et al. \cite{Corless1996} expressed it as \(^{\infty}r=W(-\ln r)/(-\ln r)\) by using Lambert's W function. \vspace{0.2in}

\noindent
Different from established extension of base \(r\) to real or complex values under integers or infinite height, extension of height \(n\)  to real or complex values is still under debate. The object of this paper is to derive rigorously the explicit formula of continuous tetration. \vspace{0.2in}

\begin{dfn}
\label{dfn_tau}
Let \(r>0\) be real constant and \(x>-2\) be real variable. Tetrational function \(\tau (x)\) is a continuous function defined as:
\[\tau(0)\colonequals 1,\hspace{15pt}\tau(x)\colonequals r^{\tau(x-1)}.\]
\end{dfn}

\begin{cly}
Independent of \(r\), every tetrational function \(\tau(x)\) goes through \(\tau(-1)=0\) and \(\tau(0)=1\).
\end{cly}
\begin{cly}
\label{cly_tau}
If a curve segment connecting between \(\tau(-1)=0\) and \(\tau(0)=1\) is defined, and is extended to other region by exponentiation or logarithm, then the whole curve satisfies Definition\ref{dfn_tau}.
\end{cly} \vspace{0.2in}

\noindent
Because of Colloray \ref{cly_tau}, it is taken for granted that \(\tau(x)\) cannot be uniquely determined without an extra requirement. It follows that there exist many unique solutions for different extra conditions. Hooshmand \cite{Hooshmand2006} applied a linear segment between \(\tau(-1)=0\) and \(\tau(0)=1\), and proved uniqueness of the function \(\tau(x)=\exp_{r}^{\lfloor x \rfloor+1}\{x\}\), where \(\lfloor x \rfloor\) is the floor function of \(x\) and \(\{x\}=x-\lfloor x \rfloor\) is the sawtooth function. He found that the first derivative exists at the connecting points of \(x=n\in\mathbb{N}\) only for the case of \(r=e\), and the second derivative does not exist at the connecting points. Kouznetsov \cite{Kouznetsov2009} gave a numerical solution for \(r>e^{1/e}\) in the complex plane under the requirement that  \(\tau(x+iy)\) is holomorphic, and Paulsen and Cowgill \cite{Paulsen2017} proved that such a holomorphic solution is unique.\vspace{0.2in}

\noindent
In this paper, we analytically derive an unique solution only from the minimum requirement of differentiability: the delay-differential equation in Lemma \ref{lem_delaydiff}.
\begin{lem}
\label{lem_delaydiff}
Let \(q=\ln{r}\) be a real constant.\\
If a tetrational function is differentiable, the following delay-differential equation holds at every point in \(x>-2\):
\[\frac{d\tau(x)}{dx}=q\tau(x)\frac{d\tau(x-1)}{dx}\]
\end{lem} \vspace{0.2in}
\begin{proof}
The derivative of \(\tau(x)\) in Definition \ref{dfn_tau}.
\end{proof}\vspace{0.2in}

\noindent
This paper consists of eight sections. In Section \ref{sec_multiple}, a new function, named multi-tetrational function, is defined to express the restriction placed in the delay-differential equation. In Section \ref{sec_Taylor}, the tetrational function and multi-tetrational function are then expressed in Taylor series with coefficients characteristic to these functions. In Section \ref{sec_trans}, based on the general forms of coefficients, the Taylor series are transformed into new expressions. In section \ref{sec_explicit} explicit formulae of tetrational functions and multi-tetrational functions, the main results of this paper, are derived. In Section \ref{sec_complex}, Analytical properties as well as the extension to the complex heights and bases are discussed. In Section \ref{sec_rule}, calculation rules of tetration are presented. Concluding remarks are given in section \ref{sec_summary}.

\section{Multi-tetrational function}
\label{sec_multiple}
In this section we define a key function, and by using this, we transform the delay-differential function into an ordinary differential equation.\vspace{0.2in}

\noindent
Let us consider a continuous function \(\mu(x)\) which goes through the following discrete values:
\[\mu(0)=\tau(0)=1\]
\[\mu(1)=\tau(1)\tau(0)=r\]
\[\mu(2)=\tau(2)\tau(1)\tau(0)=r^{r}r\]
\[\vdots\]
\[\mu(n)=\tau(n)\mu(n-1)\]
It is reasonable to define \(\mu(x)\) by replacing \(n\) with \(x\).\vspace{0.2in}

\noindent
\begin{dfn}
\label{dfn_mu}
A continuous function \(\mu(x)\) is called a multi-tetrational function if it satisfies the following relations:
\[\mu(0)\colonequals\tau(0),\hspace{15pt}\mu(x)\colonequals\tau(x)\mu(x-1)\]
\end{dfn}
\begin{lem}
\label{lem_mu_value}
The equation \(\mu(-1)=\mu(0)=1\) holds independent of \(r\).
\end{lem}
\begin{proof}
From Definition \ref{dfn_tau}, we get \(\mu(0)=\tau(0)=1\). Assigning \(x=0\) to second equation in Definition\ref{dfn_mu} gives \(\mu(-1)=\mu(0)/\tau(0)=1\).
\end{proof}\vspace{0.2in}

\noindent
The delay-differential equation in Lemma \ref{lem_delaydiff} is transformed in to an ordinary differential equation including \(\mu(x)\) as in Theorem \ref{thm_dif_tau_mu}, which places strong restriction on both \(\tau(x)\) and \(\mu(x)\) .\vspace{0.2in}

\begin{thm}
\label{thm_dif_tau_mu}
Let \(q=\ln{r}\) be a real constant.\\
Let \(\omega\in\mathbb{C}\) be q-dependent constant. \\
The delay-defferentail equation \ref{lem_delaydiff} is equivalent to the differential equation:
\[\frac{d\tau(x)}{dx}=\omega q^{x}\mu(x)\]
\end{thm}
\begin{proof}
By assigning integer heights \(x=n\in\mathbb{N}\) to the equation \ref{lem_delaydiff}, we get:
\[\tau'(0)=q\tau(0)\tau'(-1)\]
\[\tau'(1)=q\tau(1)\tau'(0)=q^{2}\tau(1)\tau(0)\tau'(-1)\]
\[\tau'(2)=q\tau(2)\tau'(1)=q^{3}\tau(2)\tau(1)\tau(0)\tau'(-1)\]
\[\vdots\]
\[\tau'(n)=q^{n+1}\mu(n)\tau'(-1).\]
By defining \(\omega\colonequals{q}\tau'(-1)\), we get
\[\tau'(n)={\omega}q^{n}\mu(n).\]
By definitions \ref{dfn_tau} and \ref{dfn_mu}, the continuous functions \(\tau(x)\), \(\mu(x)\) and \(q^{x}\) go through all of the discrete points \(\tau(n)\), \(\mu(n)\) and \(q^{n}\), respectively. Therefore equation above is extended to real values:
\[\tau'(x)={\omega}q^{x}\mu(x).\]
Conversely, if \(\tau'(x)={\omega}q^{x}\mu(x)\) and \(\mu(x)=\tau(x)\mu(x-1)\) are given, then \[\frac{\tau'(x)}{\tau'(x-1)}=\frac{{\omega}q^{x}\mu(x)}{{\omega}q^{x-1}\mu(x-1)}=\frac{q\mu(x)}{\mu(x-1)}=q\tau(x).\]
Therefore the delay-defferential equation \(\tau'(x)=q\tau(x)\tau'(x-1)\) is reproduced.\\
\end{proof}

\section{Taylor series}
\label{sec_Taylor}
In this section, we will find general forms of the coefficients in Taylor series of the tetrational function and multi-tetrational function. \vspace{0.2in}

\noindent
First we express \(\tau(x)\) and \(\mu(x)\) as Taylor series.

\begin{dfn}
\label{dfn_taylor}
Let \(a_n, b_n, c_n\) and \(d_n\) be real Taylor coefficients.
Taylor series of tetrational and multi-tetrational functions  near the points \(x=0\) or \(x=-1\) are expressed as:
\[\tau(x)\colonequals\sum_{n=0}^{\infty}\frac{1}{n!}{a_n}{x^{n}},\hspace{15pt}\tau(x-1)\colonequals\sum_{n=0}^{\infty}\frac{1}{n!}{b_n}{x^{n}},\]
\[\mu(x)\colonequals\sum_{n=0}^{\infty}\frac{1}{n!}{c_n}{x^{n}},\hspace{15pt}\mu(x-1)\colonequals\sum_{n=0}^{\infty}\frac{1}{n!}{d_n}{x^{n}},\]
\end{dfn}

\begin{cly}
The following coefficients are given as follows independent of \(r\):
\[a_0=1,\hspace{15pt}b_0=0,\hspace{15pt}c_0=1,\hspace{15pt}d_0=1\]
\end{cly}
\begin{proof}
Corollary \ref{cly_tau} gives \(a_0=\tau(0)=1\) and \(b_0=\tau(-1)=0\). Lemma \ref{lem_mu_value} gives \(c_0=\mu(0)=1\), and  \(d_0=\mu(-1)=1\).
\end{proof}\vspace{0.2in}

\noindent
Next we find the correlations of coefficients \(a_n, b_n, c_n\) and \( d_n\).

\begin{lem}
\label{lem_coefficient_abcd}
Let \(s=\ln{q}\) be the real or complex constant dependent on \(q\), where \(s\) is the principal value for \(q<0\).\\
Let \(\binom{n}{k}\) be the binomial coefficient.\\
Taylor coefficients \(a_n, b_n, c_n\),and \( d_n\) for \(n\geq1\) are related to each other in the following ways:
\begin{eqnarray}
a_n&=&\omega\sum_{k=1}^{n}\binom{n-1}{k-1}c_{k-1}s^{n-k},\\
qb_n&=&\omega\sum_{k=1}^{n}\binom{n-1}{k-1}d_{k-1}s^{n-k},\\
c_n&=&\sum_{k=0}^{n}\binom{n}{k}a_{n-k}d_{k}.
\end{eqnarray}
\end{lem}
\begin{proof}
Repeated differentiation of the equation in Theorem \ref {thm_dif_tau_mu} and assigning \(x=0\) and  \(x=-1\) to it give relations (1) and (2). The relation (3) is given by assigning \(x=0\) to the equation in Definition \ref{dfn_mu} after repeated differentiation of it.
\end{proof}
\noindent
From Lemma \ref{lem_delaydiff}, we can derive another relationship between \(a_{n}\) and \(b_{n}\) similar to Lemma \ref{lem_coefficient_abcd}. The relation, however, can be derived from equations in Lemma \ref{lem_coefficient_abcd}.\vspace{0.2in}

\noindent
Finally we show that general formulae of \(a_n, b_n, c_n\) and \( d_n\) have the common structure with the Striling numbers of the second kind \(\stirling{m}{n}\) and constants \(A_{n}\) or \(B_{n}\) e.g.,
\[a_{0}=A_{0},\]
\[a_{1}=A_{0}{\omega},\]
\[a_{2}=A_{0}s{\omega}+A_{1}{\omega}^{2},\]
\[a_{3}=A_{0}s^{2}{\omega}+3A_{1}s{\omega}^{2}+A_{2}{\omega}^{3},\]
\[a_{4}=A_{0}s^{3}{\omega}+7A_{1}s^{2}{\omega}^{2}+6A_{2}s{\omega}^{3}+A_{3}{\omega}^{4}.\]
Generally,
\[a_{0}=A_{0},\hspace{10pt}a_{n}=\sum_{k=1}^{n}\stirling{n}{k}A_{k}s^{n-k}{\omega}^{k}\hspace{10pt}(n\geq{1}).\]
The proof is given in Theorem \ref{thm_coefficient_AB}.

\begin{thm}
\label{thm_coefficient_AB}
Let \(A_{n}\) and \(B_{n}\) be a real constants having the relation:
\[A_{0}\colonequals 1,\hspace{10pt}A_{1}\colonequals 1,\hspace{10pt}B_{0}\colonequals 0,\hspace{10pt}B_{1}\colonequals 1,\]
\[A_{n}\colonequals\sum_{k=1}^{n}\binom{n}{k-1}A_{n-k}B_{k}.\hspace{15pt}(n\geq 2)\]
Let \(\stirling{n}{m}\) be the Stirling numbers of the second kind.\\
Then the coefficients \(a_n, b_n, c_n\),and \( d_n\) for \(n\geq 1\) with the following formulae satisfy the relations given in Lemma \ref{lem_coefficient_abcd}:
\[a_{n}=\sum_{k=1}^{n}\stirling{n}{k}A_{k}s^{n-k}\omega^{k},\hspace{15pt}qb_{n}=\sum_{k=1}^{n}\stirling{n}{k}B_{k}s^{n-k}\omega^{k},\]
\[c_{n}=\sum_{k=1}^{n}\stirling{n}{k}A_{k+1}s^{n-k}\omega^{k},\hspace{10pt}d_{n}=\sum_{k=1}^{n}\stirling{n}{k}B_{k+1}s^{n-k}\omega^{k}.\]
\end{thm}
\begin{proof}
We shall prove this by induction.\\
For \(n=1\), relations in Lemma \ref{lem_coefficient_abcd} (1), (2) and (3) are \(a_{1}=\omega\), \(qb_{1}=\omega\) and \(c_{1}=a_{1}d_{0}+a_{0}d_{1}\), respectively.
Clearly, the formulae of \(a_{1}=A_{1}\omega=\omega\) and \(qb_{1}=B_{1}\omega=\omega\) satisfy the first and second relations above, respectively. The relation \ref{lem_coefficient_abcd} (3) also holds since the left hand side is \(c_{1}=A_{2}\omega\), and the right hand side is \(a_{1}d_{0}+a_{0}d_{1}=(A_{1}B_{1}+A_{0}B_{2})\omega=A_{2}\omega\).\vspace{0.1in}

\noindent
Then suppose the given formulae hold for \(n=k\). \\
By using \(c_{n}\) \((n\leq{k})\) and the following relationship \cite{Abramowitz1964}\cite{Graham1994}: \[\stirling{n+1}{k+1}=\sum_{j=k}^{n}\binom{n}{j}\stirling{j}{k},\] the relation of \ref{lem_coefficient_abcd} (1) for \(n=k+1\) is expressed as
\[a_{k+1}=\omega\sum_{j=1}^{k+1}\binom{k}{j-1}c_{j-1}s^{k+1-j}\]
\[={\omega}\left[\binom{k}{0}c_{0}s^{k}+\binom{k}{1}c_{1}s^{k-1}+\cdots+\binom{k}{k}c_{k}\right]\]
\[=\binom{k}{0}\stirling{0}{0}A_{1}s^{k}{\omega}+\binom{k}{1}\stirling{1}{1}A_{2}{\omega}^{2}s^{k-1}+\cdots+\binom{k}{k}\left(\sum_{j=1}^{k}\stirling{k}{j}A_{j+1}s^{k-j}\omega^{k+1}\right)\]
\[=\binom{k}{0}\stirling{0}{0}A_{1}s^{k}{\omega}+\left[\binom{k}{1}\stirling{1}{1}+\binom{k}{2}\stirling{2}{1}\right]A_{2}s^{k-1}{\omega}^{2}+\cdots+\binom{k}{k}\stirling{k}{k}A_{k+1}{\omega}^{k+1}\]
\[=\stirling{k+1}{1}A_{1}s^{k}{\omega}+\stirling{k+1}{2}A_{2}s^{k-1}{\omega}^2+\cdots+\stirling{k+1}{k+1}A_{k+1}{\omega}^{k+1}\]
Therefore, if the relation holds for \(n=k\), it holds for \(n=k+1\).
It follows that the first relation is generally satisfied by the given formulae. \\
Similarly, proof for the relation of \ref{lem_coefficient_abcd} (2) is given by replacing \(a_{n}\) and \(c_{n}\) with \(qb_{n}\) and \(d_{n}\), respectively. Again, if the relation holds for \(n=k\), it holds for \(n=k+1\). \\
Therefore relation \ref{lem_coefficient_abcd} (1) and (2) are generally satisfied by given expressions.
Next we check the given formulae of \(a_{n}\), \(c_{n}\) and \(d_{n}\) generally satisfy the relation of \ref{lem_coefficient_abcd} (3) at any \(k\).
\[c_{k}=\sum_{j=0}^{k}\binom{k}{j}a_{k-j}d_{j}\]
\[=\binom{k}{0}a_{k}d_{0}+\binom{k}{1}a_{k-1}d_{1}+\cdots+\binom{k}{k}a_{0}d_{k}\]
By replacing each \(a_{n}\) and \(d_{n}\) by the given formulae, we have
\[c_{k}=\binom{k}{0}\left[\stirling{k}{1}A_{1}s^{k-1}{\omega}+\cdots+\stirling{k}{k}A_{k}s^{0}{\omega}^{k}\right]B_{1}\]
\[+\binom{k}{1}\left[\stirling{k-1}{1}A_{1}s^{k-2}{\omega}+\cdots+\stirling{k-1}{k-1}A_{k-1}s^{0}{\omega}^{k-1}\right]B_{2}\omega+\cdots\]
\[+\binom{k}{k}A_{0}\left[\stirling{k}{1}B_{2}s^{k-1}{\omega}+\cdots+\stirling{k}{k}B_{k+1}s^{0}{\omega}^{k}\right].\]
By rearrangement, we have
\[c_{k}=\left[\binom{k}{0}\stirling{k}{1}A_{1}B_{1}+\binom{k}{k}\stirling{k}{1}A_{0}B_{2}\right]s^{k-1}\omega\]
\[+\left[\binom{k}{0}\stirling{k}{2}A_{2}B_{1}+\sum_{h=0}^{k-2}\binom{k}{j+h}\stirling{k-j-h}{2-j}\stirling{j+h}{j}A_{1}B_{2}+\binom{k}{k}\stirling{k}{2}A_{0}B_{3}\right]s^{k-1}{\omega}^{2}\]
\[+\cdots\]
\[+\left[\binom{k}{0}\stirling{k}{k}A_{k}B_{1}{\omega}^{k}+\binom{k}{1}\stirling{k-1}{k-1}A_{k-1}B_{2}+\cdots+\binom{k}{k}\stirling{k}{k}A_{0}B_{k+1}\right]{\omega}^{k}.\]
By using the following relation for \(k<m\) \cite{Abramowitz1964}\cite{Graham1994}:
\[\stirling{n}{m}\binom{m}{j}=\sum_{k=0}^{n-m}\binom{n}{j+k}\stirling{n-j-k}{m-j}\stirling{j+k}{j},\]
We have
\[c_{k}=\stirling{k}{1}\left[\binom{1}{0}A_{1}B_{1}+\binom{1}{1}A_{0}B_{2}\right]s^{k}\omega\]
\[+\stirling{k}{2}\left[\binom{2}{0}A_{2}B_{1}+\binom{2}{1}A_{1}B_{2}+\binom{2}{2}A_{0}B_{3}\right]s^{k-1}{\omega}^{2}\]
\[+\cdots\]
\[+\stirling{k}{k}\left[\binom{k-1}{0}A_{k-1}B_{1}+\binom{k-1}{1}A_{k-2}B_{2}+\cdots+\binom{k-1}{k-1}A_{0}B_{k-1}\right]\]
By using the relationship \(A_{n}=\sum_{k=1}^{n}\binom{n}{k-1}A_{n-k}B_{k}\), finally we get
\[c_{k}=\stirling{k}{1}A_{2}s^{k}{\omega}+\stirling{k}{2}A_{3}s^{k-1}{\omega}^{2}+\cdots+\stirling{k}{k-1}A_{k}s{\omega}^{k-1}+\stirling{k}{k}A_{k+1}{\omega}^{k},\]
which is consistent with the expression of \(c_{k}\) given above.\\
Therefore the relation of \ref{lem_coefficient_abcd} (3) is generally satisfied by the given formulae.\\
\end{proof}\vspace{0.2in}

\noindent
Now we have the general form of Taylor series of the tetrational and muti-tetrational functions near the points \(x=0\) or \(x=-1\):
\[\tau(x)=A_{0}+A_{1}{\omega}x+\frac{1}{2}(A_{1}s{\omega}+A_{2}{\omega}^{2})x^{2}+\frac{1}{3!}(A_{1}s^{2}{\omega}+3A_{2}s{\omega}^2+A_{3}{\omega}^{3})x^{3}+\cdots,\]
\[\mu(x)=A_{1}+A_{2}{\omega}x+\frac{1}{2}(A_{2}s{\omega}+A_{3}{\omega}^{2})x^{2}+\frac{1}{3!}(A_{2}s^{2}{\omega}+3A_{3}s{\omega}^2+A_{4}{\omega}^{3})x^{3}+\cdots,\]
\[\tau(x-1)=\frac{B_{0}}{q}+\frac{B_{1}{\omega}}{q}x+\frac{1}{2}\frac{B_{1}s{\omega}+B_{2}{\omega}^{2}}{q}x^{2}+\frac{1}{3!}\frac{B_{1}s^{2}{\omega}+3B_{2}s{\omega}^2+B_{3}{\omega}^{3}}{q}x^{3}+\cdots,\]
\[\mu(x-1)=B_{1}+B_{2}{\omega}x+\frac{1}{2}(B_{2}s{\omega}+B_{3}{\omega}^{2})x^{2}+\frac{1}{3!}(B_{2}s^{2}{\omega}+3B_{3}s{\omega}^2+B_{4}{\omega}^{3})x^{3}+\cdots.\]

\section{Tetrational series}
\label{sec_trans}
In this section, the Taylor series of the tetrational functions and multi-tetrational functions are transformed into new expressions based on the property of Stirling numbers of the second kind.\vspace{0.2in}

\noindent
First we define a key function to express new series.
\begin{dfn}
\label{dfn_phi}
Let \(\varphi(x)\) be a function defined as
\[\varphi(x)\colonequals\frac{{\omega}(q^{x}-1)}{s}.\]
\end{dfn}
\begin{cly}
\label{cly_diff_Phi}
The following equation holds:
\[\frac{d\varphi(x)}{dx}={\omega}q^{x}.\]
\end{cly}\vspace{0.2in}

\noindent
By using \(\varphi(x)\), the Taylor series are then simplified as in Theorem \ref{thm_tetrational_series}.

\begin{thm}
\label{thm_tetrational_series}
Let Taylor series of tetrational functions and mutiple tetrational functions near the points \(x=0\) or \(x=-1\) be absolutely convergent, then they are expressed in power series of \(\varphi(x)\), named tetrational series:
\[\tau(x)=\sum_{n=0}^{\infty}\frac{1}{n!}A_{n}\varphi^{n}(x),\hspace{15pt}\tau(x-1)=\frac{1}{q}\sum_{n=0}^{\infty}\frac{1}{n!}B_{n}\varphi^{n}(x),\]
\[\mu(x)=\sum_{n=0}^{\infty}\frac{1}{n!}A_{n+1}\varphi^{n}(x),\hspace{15pt}\mu(x-1)=\sum_{n=0}^{\infty}\frac{1}{n!}B_{n+1}\varphi^{n}(x).\]
\end{thm}
\begin{proof}
Since the Taylor series of tetrational function is absolutely convergent, it can be rearranged as follows:
\[\tau(x)=A_{0}+A_{1}{\omega}x+\frac{1}{2}(A_{1}s{\omega}+A_{2}{\omega}^{2})x^{2}+\frac{1}{3!}(A_{1}s^{2}{\omega}+3A_{2}s{\omega}^2+A_{3}{\omega}^{3})x^{3}+\cdots\]
\[=A_{0}+A_{1}{\omega}\left(x+\frac{1}{2}sx^{2}+\frac{1}{3!}s^{2}x^{3}+\cdots\right)+A_{2}{\omega}^{2}\left(\frac{1}{2}x^2+\frac{1}{3!}3sx^{3}+\cdots\right)+\cdots\]
\[=A_{0}+A_{1}\frac{\omega}{s}\left(sx+\frac{1}{2}s^{2}x^{2}+\frac{1}{3!}s^{2}x^{3}+\cdots\right)+A_{2}\frac{\omega^{2}}{s^{2}}\left(\frac{1}{2}s^{2}x^2+\frac{1}{3!}3s^{3}x^{3}+\cdots\right)+\cdots.\]
By using the exponential generating function of the Stiring numbers of the second kind \cite{Abramowitz1964}\cite{Graham1994}:
\[\frac{1}{k!}(e^{x}-1)^{k}=\sum_{n=k}^{\infty}\stirling{n}{k}\frac{x^{n}}{n!},\]
we have
\[\tau(x)=A_{0}+A_{1}\left[\frac{{\omega}(e^{sx}-1)}{s}\right]+\frac{1}{2}A_{2}\left[\frac{{\omega}(e^{sx}-1)}{s}\right]^{2}+\frac{1}{3!}A_{3}\left[\frac{{\omega}(e^{sx}-1)}{s}\right]^{3}+\cdots.\]
Since \(s=\ln{q}\), we can replace \(e^{sx}\) with \(q^{x}\). The Taylor series
are expressed as the power series of \(\varphi(x).\)
Similarly, we can transform \(\tau(x-1)\), \(\mu(x)\) and \(\mu(x-1)\) into the tetrational series:
\begin{align*}
\tau(x)&=A_{0}+A_{1}\varphi(x)+\frac{1}{2}A_{2}\varphi^{2}(x)+\frac{1}{3!}A_{3}\varphi^{3}(x)+\cdots,\\
\mu(x)&=A_{1}+A_{2}\varphi(x)+\frac{1}{2}A_{3}\varphi^{2}(x)+\frac{1}{3!}A_{4}\varphi^{3}(x)+\cdots,\\
\tau(x-1)&=\frac{B_{0}}{q}+\frac{B_{1}}{q}\varphi(x)+\frac{1}{2}\frac{B_{2}}{q}\varphi^{2}(x)+\frac{1}{3!}\frac{B_{3}}{q}\varphi^{3}(x)+\cdots,\\
\mu(x-1)&=B_{1}+B_{2}\varphi(x)+\frac{1}{2}B_{3}\varphi^{2}(x)+\frac{1}{3!}B_{4}\varphi^{3}(x)+\cdots.
\end{align*}
\end{proof}\vspace{0.2in}

\begin{cly}
\label{cly_dtau_dphi_mu}
The following equations hold near \(x=0\):
\[\frac{\partial\tau(x)}{\partial\varphi}=\mu(x),\hspace{10pt}\frac{\partial\tau(x-1)}{\partial\varphi}=q^{-1}\mu(x-1).\]
\end{cly}\vspace{0.2in}

\noindent
We check the consistency of the formulae as shown in Lemma \ref{lem_consistency}.
\begin{lem}
\label{lem_consistency}
Formulae \ref{thm_tetrational_series} satisfy the differential equation \ref{thm_dif_tau_mu}:
\[\frac{d\tau(x)}{dx}={\omega}q^{x}\mu(x).\]
\end{lem}
\begin{proof}
Since \ref{cly_diff_Phi} and \ref{cly_dtau_dphi_mu}, with chain rule,
\[\frac{d\tau(x)}{dx}=\frac{\partial\tau(x)}{\partial\varphi}\frac{\partial\varphi(x)}{\partial{x}}={\omega}q^{x}\mu(x),\]
\[\frac{d\tau(x-1)}{dx}=\frac{\partial\tau(x-1)}{\partial\varphi}\frac{\partial\varphi(x)}{\partial{x}}={\omega}q^{x-1}\mu(x-1).\]
\end{proof}\vspace{0.2in}

\noindent
We further describe the general case of tetrational series with an integer shift \(n\) as in Lemma \ref{lem_integer_shift}. The formulae above are special case of \(n=0\).

\begin{lem}
\label{lem_integer_shift}
Let \(0\leq{x}<1\) be the real heights.\\
Let \(n\) and \(m\) be integers and natural numbers, respectively.\\
Let \(A_{m}^{[n]}\) and \(B_{m}^{[n]}\) be n-dependent real constants.\\
The tetrational and multi-tetrational functions, having horizontal translation by integers \(n\) or \(n-1\), are expressed as:
\begin{align*}
\tau(x+n)&=q^{n}\sum_{k=0}^{\infty}\frac{x^{k}}{k!}\sum_{j=1}^k\stirling{k}{j}A_{j}^{[n]}s^{k-j}{\omega}^{j}=q^{n}\sum_{k=0}^{\infty}\frac{1}{k!}A_{k}^{[n]}{\varphi}^{k}(x),\\
\mu(x+n)&=\sum_{k=0}^{\infty}\frac{x^{k}}{k!}\sum_{j=1}^k\stirling{k}{j}A_{j+1}^{[n]}s^{k-j}{\omega}^{j}=\sum_{k=0}^{\infty}\frac{1}{k!}A_{k+1}^{[n]}{\varphi}^{k}(x),\\
\tau(x+n-1)&=q^{n-1}\sum_{k=0}^{\infty}\frac{x^{k}}{k!}\sum_{j=1}^k\stirling{k}{j}B_{j}^{[n]}s^{k-j}{\omega}^{j}=q^{n-1}\sum_{k=0}^{\infty}\frac{1}{k!}B_{k}^{[n]}{\varphi}^{k}(x),\\
\mu(x+n-1)&=\sum_{k=0}^{\infty}\frac{x^{k}}{k!}\sum_{j=1}^k\stirling{k}{j}B_{j+1}^{[n]}s^{k-j}{\omega}^{j}=\sum_{k=0}^{\infty}\frac{1}{k!}B_{k+1}^{[n]}{\varphi}^{k}(x).
\end{align*}
Coefficients have the relationships of:
\[A_{0}^{[n]}=\tau(n),\hspace{15pt}B_{0}^{[n]}=\tau(n-1),\]
\[A_m^{[n]}=q^{n}\sum_{k=1}^{m}\binom{m-1}{k-1}A_{m-k}^{[n]}B_{k}^{[n]}\hspace{10pt}(m\geq{1}).\]
Generally, \(\tau(x)\) and \(\mu(x)\) having horizontal translation by \(n\) have the relationship of:
\[\frac{\partial \tau(x+n)}{\partial\varphi}=q^{n}\mu(x+n).\]
\end{lem}
\begin{proof}
Obviously the tetrational series satisfy \(\partial\tau(x+n)/\partial\varphi=q^{n}\mu(x+n)\). It follows that these series satisfy the differential equation \(\tau'(x+n)={\omega}q^{x+n}\mu(x+n)\), as similarly in \ref{lem_consistency}. Tetrational series and the Taylor series are converted into each other with the exponential generating function of the Stirling numbers of the second kind as in Theorem \ref{thm_tetrational_series}.\\
\noindent
The relationships of coefficients are given by assigning \(x=0\) to the repeated partial derivatives by \(\varphi\) of \(\mu(x+n)=\tau(x+n)\mu(x+n-1)\).
\end{proof}\vspace{0.2in}

\section{Explicit formulae}
\label{sec_explicit}
In this section, we derive the explicit formulae of tetrational and muti-tetrational function by determining the coefficients \(A_{n}\) and \(B_{n}\).\vspace{0.2in}

\begin{lem}
\label{lem_mu_piecewise}
Let \(0\leq{x}<1\).
The following equations hold:
\[\mu(x-1)=1,\]
\[\mu(x)=\tau(x),\]
\[\mu(x+n)=\prod_{k=0}^{n}\tau(x+k).\]
\end{lem}
\begin{proof}
We defined \(\mu(x)\) in Definition \ref{dfn_mu} as a continuous function satisfying \[\mu(x)=\tau(x)\mu(x-1),\]
which goes through the following discrete points:
\[\mu(0)=\tau(0)\]
\[\mu(1)=\tau(1)\tau(0)\]
\[\mu(2)=\tau(2)\tau(1)\tau(0)\]
\[\vdots\]
\[\mu(n)=\prod_{k=0}^{n}\tau(k)=\tau(n)\mu(n-1).\]
Similarly, \(\mu(x)\) also goes through the discrete points with changes of \(0\leq{x}<1\) as follows, since continuous functions \(\mu(x)\) and \(\tau(x)\) go through \(\mu(n)\) and \(\tau(n)\) respectively:
\[\mu(x)=\tau(x)\]
\[\mu(x+1)=\tau(x+1)\tau(x)\]
\[\mu(x+2)=\tau(x+2)\tau(x+1)\tau(x)\]
\[\vdots\]
\[\mu(x+n)=\prod_{k=0}^{n}\tau(x+k)=\tau(x+n)\mu(x+n-1)\]
The relationship \(\mu(x+n)=\tau(x+n)\mu(x+n-1)\) is consistent with the definition \ref{dfn_mu}.\\
\noindent
Since from the definition \(\mu(x-1)=\mu(x)/\tau(x)\) and \(\mu(x)=\tau(x)\) as given above, we have
\[\mu(x-1)=\frac{\tau(x)}{\tau(x)}=1.\]
\end{proof}

\begin{cly}
Tetrational function \(\tau(x)\) and \(\mu(x)\) are piecewise connected functions.
\end{cly}
\begin{proof}
Obviously from Lemma \ref{lem_mu_piecewise}, multi-tetrational function \(\mu(x)\) is a piecewise connected function. Since tetrational function have the relationship of \(\tau'(x)=\omega{q}^x\mu(x)\), tetrational function \(\tau(x)\) is also a piecewise connected function.
\end{proof}

\begin{lem}
\label{lem_AB_value}
The coefficients \(A_{n}\) and \(B_{n}\) are determined as
\[A_{n}=1,\]
\[B_{n\neq{1}}=0,\hspace{15pt}B_{1}=1.\]
\end{lem}
\begin{proof}
From Theorem \ref{thm_tetrational_series},
\[\tau(x)=\sum_{n=0}^{\infty}\frac{1}{n!}A_{n}\varphi^{n}(x),\hspace{10pt}\mu(x)=\sum_{n=0}^{\infty}\frac{1}{n!}A_{n+1}\varphi^{n}(x).\]
Since \(\partial{\tau(x)}/\partial{\varphi}=\mu(x)\) from \ref{cly_dtau_dphi_mu} and \(\tau(x)=\mu(x)\) from \ref{lem_mu_piecewise}, we have
\[1=A_{0}=A_{1}=A_{2}=\cdots.\]
From the relationship of coefficients \(A_{n}=\sum_{k=1}^{n}\binom{n}{k-1}A_{n-k}B_{k}\), we get \(B_{n\neq{1}}=0\) and \(B_{1}=1\).\\
\end{proof}
\noindent
Similarly, we can get the values of \(A_{m}^{[n]}\) and \(B_{m}^{[n]}\) given in Lemma \ref{lem_integer_shift}.\vspace{0.2in}

\begin{lem}
\label{lem_omega_value}
The constant \(\omega\) is determined as
\[\omega=\frac{sq}{q-1}.\]
\end{lem}
\begin{proof}
Since we determined the values of \(B_{n}\) as in Lemma \ref{lem_AB_value}, and with the expression of \(\varphi(x)\) in Definition \ref{dfn_phi}, we can express \(\tau(x-1)\) for \(0\leq{x}<1\) as
\[\tau(x-1)=\frac{1}{q}\varphi(x)=\frac{\omega(q^{x}-1)}{sq}\]
Then by assigning \(x=1\), we get the equation \(\omega(q-1)/sq=1\) to derive \(\omega\).
\end{proof}\vspace{0.2in}

\noindent
Now in Theorem \ref{thm_explicit} we can derive the explicit solution for the delay-differential equation \ref{lem_delaydiff}. The domains of the tetrational function as well as the multi-tetrational function are naturally extended by applying the sawtooth function.

\begin{thm}
\label{thm_explicit}
Let \(r\) be real constant.\\
Let \(x\) be real variable heights in \(x\in\mathbb{R}\setminus\{n:n\in\mathbb{Z}, n\leq{-2}\}.\)\\
Let \(\lfloor{x}\rfloor\) and \(\{x\}\colonequals{x}-\lfloor{x}\rfloor\) be the floor function and the sawtooth function of \(x\), respectively.\\
Let \([x]_q\) be the q-analog of \(x\) defined as
\[[x]_{q}\colonequals\frac{q^{x}-1}{q-1},\]
where \(q=\ln{r}\) is the principal branch if \(r<0\).\\
Let \(\exp_{r}^{n}=\log_{r}^{-n}\) be the extended iterated exponential operator defined as
\[\exp_{r}^{n}{f}\colonequals
\begin{cases}
\underset{n}{\underbrace{\exp_r\exp_r\cdots\exp_r}} {f} & \hspace{15pt}(n\geq{0})\\
\hspace{3pt}\underset{n}{\underbrace{\log_r\log_r\cdots\log_r}} {f} & \hspace{15pt}(n<{0}),
\end{cases}
\]
where \(f\) is the operand.\\
Let \(\prod_{k=m}^{n}\) be the extended product operator defined as
\[\prod_{k=m}^{n}f_{k}\colonequals
\begin{dcases}
f_{m}f_{m+1}{\cdots}f_{n-1}f_{n} & (0\leq{m}\leq{n})\\
\frac{1}{f_{m+1}f_{m+2}{\cdots}f_{n+1}} & (n\leq{m}<0),
\end{dcases}
\]
where \(f_{k}\) are the operands.\\
Then the tetrational function and the multi-tetratinal functions are expressed as
\[{^x}r=\tau(x)=\exp_{r}^{\lfloor{x}\rfloor+1}[\{x\}]_q\]
\[\mu(x)=\prod_{k=0}^{\lfloor{x}\rfloor}\exp_{r}^{k+1}[\{x\}]_q=q^{-\lfloor{x}\rfloor-1}\frac{\partial\tau(x)}{\partial[\{x\}]_q}\]
\end{thm}
\begin{proof}
Let us consider the case of \(0\leq x<1\).\\
From Definition \ref{dfn_phi} and Lemma \ref{lem_omega_value}, we have
\[\varphi(x)=\frac{q(q^{x}-1)}{q-1}.\]
By using the notation of q-analog of the variable, we can express it as
\[\varphi(x)=q[x]_q.\]
Then the tetrational and multi-tetrational functions for \(n\geq{0}\) are given as follows since \(\tau(x)=\sum_{n=0}^{\infty}{(q[x]_q)}^{n}/n!=\exp_{r}[x]_{q}\) in Theorem \ref{thm_tetrational_series} and \(\tau(x+n+1)=\exp_{r}{\tau(x+n)}\) as well as \(\mu(x+n)=\tau(x+n)\mu(x+n-1)\):
\begin{align*}
\tau(x)=\exp_{r}[x]_q, & \hspace{15pt}\mu(x)=\exp_{r}[x]_q\\
\tau(x+1)=\exp_{r}^{2}[x]_q, & \hspace{15pt}\mu(x+1)=\exp_{r}^{2}[x]_q\exp_{r}[x]_q\\
\tau(x+2)=\exp_{r}^{3}[x]_q, & \hspace{15pt}\mu(x+2)=\exp_{r}^{3}[x]_q\exp_{r}^{2}[x]_q\exp_{r}[x]_q\\
& \hspace{5pt}\vdots\\
\tau(x+n)=\exp_{r}^{n+1}[x]_q, & \hspace{15pt}\mu(x+n)=\prod_{k=0}^{n}\exp_{r}^{k+1}[x]_q
\end{align*}
In the case for \(n<0\). The relation \(\tau(x-1)=[x]_q\) and \(\tau(x-n-1)=\log_{r}\tau(x-n)\), as well as \(\mu(x-1)=1\) and \(\mu(x-n)=\mu(x-n+1)/\tau(x-n+1)\), give the following series:
\begin{align*}
\tau(x-1)=[x]_q, & \hspace{15pt}\mu(x-1)=\frac{\tau(x)}{\tau(x)}=1,\\
\tau(x-2)=\log_{r}[x]_q, & \hspace{15pt}\mu(x-2)=\frac{\tau(x)}{\tau(x)\tau(x-1)}=\frac{1}{[x]_q},\\
\tau(x-3)=\log_{r}^{2}[x]_q, & \hspace{15pt}\mu(x-3)=\frac{\tau(x)}{\tau(x)\tau(x-1)\tau(x-2)}=\frac{1}{[x]_q\log_{r}[x]_q},\\
& \hspace{5pt}\vdots\\
\tau(x-n)=\log_{r}^{n+1}[x]_q, & \hspace{15pt}\mu(x-n)=\frac{\prod_{k=0}^{0}\tau(x+k)}{\prod_{k=1}^{n}\tau(x-k+1)}=\prod_{k=0}^{-n}\tau(x+k).
\end{align*}
In the last line, we applied the extended operator.\\
Therefore we can also give the same expression for \(n<0\) as
\[\tau(x+n)=\exp_{r}^{n+1}[x]_q,\hspace{15pt}\mu(x+n)=\prod_{k=0}^{n}\exp_{r}^{k+1}[x]_q.\]
Obviously, the following equation holds for both cases of \(n\geq{0}\) and \(n<0\):
\[\mu(x+n)=q^{-n-1}\frac{\partial\tau(x+n)}{\partial[x]_q}.\]
By replacing \(x\) and \(n\) above with \(\{x\}\) and \(\lfloor{x}\rfloor\), respectively, we have the expressions of tetrational and multi-tetrational functions for \(x\in(-\infty, \infty)\):
\[\tau(x)=\exp_{r}^{\lfloor{x}\rfloor+1}[\{x\}]_q,\hspace{10pt}\mu(x)=\prod_{k=0}^{\lfloor{x}\rfloor}\exp_{r}^{k+1}[\{x\}]_q=q^{-\lfloor{x}\rfloor-1}\frac{\partial\tau(x)}{\partial[\{x\}]_q}.\]
\end{proof}

\begin{cly}
\label{cly_core}
The following equation holds:
\[
^{\{x\}}r=r^{[\{x\}]_q}.
\]
\end{cly}

\begin{proof}
Since \(0\leq{\{x\}}<1\) and \(\lfloor\{x\}\rfloor=0\), the equation is given from Theorem \ref{thm_explicit}.
\end{proof}

\noindent
Let us consider the special case of \(r=e\). If \(r\rightarrow{e}\), hence \(q=\ln{r}\rightarrow{1}\), then q-analog becomes \([\{x\}]_q\rightarrow{\{x\}}\). So we have \(\tau(x)=\exp_{r}^{\lfloor{x}\rfloor+1}\{x\}.\)
This formula is the same as Hooshmand's ultra exponential function for  \(r=e\) \cite{Hooshmand2006}.\vspace{0.2in}

\noindent
The curves of continuous tetrational functions \({^x}r=\tau(x)\) for different bases \(r\) are shown in Fig. \ref{fig:continuous tetration}. It is obvious that the tetrational functions for \(r<1\) have complex values.
\begin{figure}
\includegraphics[scale=0.4]{fig_continuous.pdf}
\centering
\caption{Behaviours of the continuous tetrational functions \({^x}r\) for different bases \(r\), where only the real component are expressed for \(r=1/e\) and \(r=e^{-e}\) (left). The tetrational functions of \(r=1/e\) and \(r=e^{-e}\) are expressed in the complex plane (right).}
\label{fig:continuous tetration}
\end{figure}
\vspace{0.2in}

\noindent
Figure \ref{fig:multiple} shows the behaviours of the mutiple tetrational functions \(\mu(x)\) for different bases \(r\). In general the functions are not differentiable at \(x=n\in\mathbb{Z}\)
\begin{figure}
\includegraphics[scale=0.4]{fig_multiple.pdf}
\centering
\caption{Behaviors of the multi-tetrational functions \(\mu(x)\) for different bases \(r\), where only the real component are expressed for \(r=1/e\) and \(r=e^{-e}\) (left). The multi-tetrational functions of \(r=1/e\) and \(r=e^{-e}\) are expressed in the complex plane (right).}
\label{fig:multiple}
\end{figure}
\vspace{0.2in}

\section{Analytical properties}
\label{sec_complex}
First, let us study the properties of \({^y}x=\tau(x,y)\) as the multivariable function of base \(x\) and height \(y\). From this point, we shall freely choose the letters to express the bases and integers.

\begin{thm}
\label{thm_total_differential}
Let \(x>0\) be real bases and let \(q=\ln{x}\).\\
Let \(y\) be real heights in \(y\in\mathbb{R}\setminus\{n:n\in\mathbb{Z}, n\leq{-2}\}.\)\\
Then the tetrational function \({^{y}}x=\tau(x,y)=\exp_{x}^{\lfloor{y}\rfloor+1}[\{y\}]_{q}\) is totally differentiable:
\[d\tau(x,y)=\frac{\partial\tau(x,y)}{\partial{x}}dx+\frac{\partial\tau(x,y)}{\partial{y}}dy.\]
\end{thm}
\begin{proof}
We shall prove this by induction.
Let \(-1\leq{y}<0\), then the tetrational function is
\[\tau(x,y)=\frac{q^{y+1}-1}{q-1}=\frac{e^{(\ln{\ln{x}})(y+1)}-1}{\ln{x}-1}.\]
Then the partial derivatives for \(-1\leq{y}<0\) are:
\[\frac{\partial\tau(x,y)}{\partial{x}}=\frac{ye^{(\ln{\ln{x}})(y+1)}}{x\ln{x}(\ln{x}-1)}-\frac{e^{(\ln{\ln{x}})(y+1)}-1}{x(\ln{x}-1)^{2}}\]
\[\frac{\partial\tau(x,y)}{\partial{y}}=\frac{(\ln{\ln{x}})e^{(\ln{\ln{x}})(y+1)}}{\ln{x}-1}\]
The following equation holds for \(-1\leq{y}<0\).
\[\frac{\partial^{2}\tau(x,y)}{\partial{x}\partial{y}}=\frac{\partial^{2}\tau(x,y)}{\partial{y}\partial{x}}=\frac{(1+(y+1)\ln{\ln{x}})e^{(\ln{\ln{x}})(y+1)}}{x\ln{x}(\ln{x}-1)}-\frac{(\ln{\ln{x}})e^{(\ln{\ln{x}})(y+1)}}{x(\ln{x}-1)^{2}}\]
Therefore
\[d\tau(x,y)=\frac{\partial\tau(x,y)}{\partial{x}}dx+\frac{\partial\tau(x,y)}{\partial{y}}dy.\]
Next, let us consider the general case.\\
From the Definition \ref{dfn_tau}:
\[{^{y+1}}x=\tau(x,y+1)=e^{(\ln{x})\tau(x,y)},\]
the partial derivatives of this general equation are
\[\frac{\partial\tau(x,y+1)}{\partial{x}}=\tau(x,y+1)\left[\frac{\tau(x,y)}{x}+(\ln{x})\frac{\partial\tau(x,y)}{\partial{x}}\right],\]
\[\frac{\partial\tau(x,y+1)}{\partial{y}}=(\ln{x})\tau(x,y+1)\frac{\partial\tau(x,y)}{\partial{y}}.\]
\noindent
Let \(\tau(x,y)\) be defined in \(n\leq{y}<n+1\), and \(\tau(x,y+1)\) be defined in \(n+1\leq{y+1}<n+2\).\\
\noindent
Suppose the following equation holds for \(n\leq{y}<n+1\):
\[\frac{\partial^{2}\tau(x,y)}{\partial{x}\partial{y}}=\frac{\partial^{2}\tau(x,y)}{\partial{y}\partial{x}}\]
Since \(\tau(x,y)\) in Theorem \ref{thm_explicit} is determined so as to be continuous at integer heights \(y\), and differentiation along base \(x\) is not affected by the existence of boundary,
obviously the connection at \(y=n\) is continuous as follows:
\[\underset{y\rightarrow{n+1}}\lim\frac{\partial\tau(x,y)}{\partial{x}}=\underset{y\rightarrow{n}}\lim\frac{\partial\tau(x,y+1)}{\partial{x}}=\frac{\partial\tau(x,n+1)}{\partial{x}}\]
\[\underset{y\rightarrow{n+1}}\lim\frac{\partial\tau(x,y)}{\partial{y}}=\underset{y\rightarrow{n}}\lim\frac{\partial\tau(x,y+1)}{\partial{y}}=\frac{\partial\tau(x,n+1)}{\partial{y}}\]
\noindent
The following equation holds:
\[\frac{\partial^{2}\tau(x,y+1)}{\partial{x}\partial{y}}=\frac{\partial^{2}\tau(x,y+1)}{\partial{y}\partial{x}}\]
\[=(\ln{x})\tau(x,y+1)\frac{\partial\tau(x,y)}{\partial{y}}\left[\frac{\tau(x,y)}{x}+(\ln{x})\frac{\partial\tau(x,y)}{\partial{x}}\right]\]
\[+\tau(x,y+1)\left[\frac{1}{x}\frac{\partial\tau(x,y)}{\partial{y}}+(\ln{x})\frac{\partial^{2}\tau(x,y)}{\partial{y}\partial{x}}\right]\]
Hence:
\[d\tau(x,y+1)=\frac{\partial\tau(x,y+1)}{\partial{x}}dx+\frac{\partial\tau(x,y+1)}{\partial{y}}dy.\]
In general, therefore, tetrational function is totally differentiable.\\
\end{proof}\vspace{0.2in}

\noindent
The property above allow us to define the inverse operations of tetration, super-logarithm or super-root, in monotonic changing range: base \(r\geq{1}\) and height \(h>-2\).\vspace{0.2in}

\noindent
Next we study the behaviours of the tetrational function \({^{h}x}=\tau(x,h)\) with a variable base \(x\) and a constant height \(h\).\vspace{0.2in}

\begin{lem}
\label{lem_(^n)x:x=0}
let \(n\) be the integer constant.\\
Let \(x>0\) be the real variable.
\[
\underset{x\rightarrow{0}}\lim(^{n}x)=
\begin{cases}
0 & (n: \mathrm{even})\\
1 & (n: \mathrm{odd})
\end{cases}
\]
\end{lem}
\begin{proof}
We shall prove this by induction.
From Definition \ref{dfn_tau}, the relations \(^{-1}x=0\) and \(^{0}x=1\) hold for \(x\rightarrow{0}\).\\
Obviously we have \(\underset{x\rightarrow{0}}{\lim}(\ln{x})=-\infty\).\\
If \(\underset{x\rightarrow{0}}{\lim}({^{2n}}x)=1\), then
\[\underset{x\rightarrow{0}}{\lim}[(^{2n}x)\ln{x}]=-\infty.\]
It follows
\[\underset{x\rightarrow{0}}{\lim}({^{2n+1}}x)=\underset{x\rightarrow{0}}{\lim}\left[e^{(^{2n}x)\ln{x}}\right]=0,\]
If \(\underset{x\rightarrow{0}}{\lim}({^{2n+1}}x)=0\), then
\[\underset{x\rightarrow{0}}{\lim}[(^{2n+1}x)\ln{x}]=\underset{x\rightarrow{0}}{\lim}[\ln[x^{(^{2n+1}x)}]]=\underset{x\rightarrow{0}}{\lim}[\ln{(^{2n}x)}]=0.\]
It follows
\[\underset{x\rightarrow{0}}{\lim}({^{2n+2}}x)=\underset{x\rightarrow{0}}{\lim}\left[e^{(^{2n+1}x)\ln{x}}\right]=1.\]
Therefore the relation generally holds.\\
\end{proof}\vspace{0.2in}
\noindent
The relation in Lemma \ref{lem_(^n)x:x=0} can be confirmed by the curves shown in Figure \ref{fig:(^n)x_real} (left). The tetrational functions \({^n}x\) for integers \(n\) have real values and go to 0 or 1 for even \(n\) or odd \(n\) respectively under \(x\rightarrow{0}\). As shown in the right graph of Figure \ref{fig:(^n)x_real}, in the general case including non-integers \(h\), the functions \({^h}x\) have complex values and non-monotonic behaviours for \(x<1\) while the functions have real values and monotonic changes for \(x\geq{1}\).
\begin{figure} [ht]
\includegraphics[scale=0.35]{fig_n_x_real.pdf}
\centering
\caption{Behaviors of \({^n}x\) for integer heights \(-1\leq{n}\leq{4}\) (left). The tetrational function \({^h}x\) for various heights \(-1\leq{h}\leq{3}\) with an interval of 0.2 (right). In general, the functions have complex values for non-integer heights if \(x<1\), and only the real components are shown in the right figures.}
\label{fig:(^n)x_real}
\end{figure}

\begin{figure}[ht]
\includegraphics[scale=0.5]{fig_n_x_complex.pdf}
\centering
\caption{The trajectories of tetrational functions \({^h}x\) in complex plane for non-integer heights \(h\) with an interval of 0.1. (a) \(-1<h<0\), (b) \(0<h<1\), (c) \(1<h<2\), (d) \(2<h<3\). }
\label{fig:(^n)x_complex}
\end{figure}

\noindent
In the right graph of Figure \ref{fig:(^n)x_real}, only the real components are shown, while trajectories for \(x<1\) in complex plane are given in Figures \ref{fig:(^n)x_complex}. It is obvious that \({^h}x\) do not always go to 0 or 1 for \(x\rightarrow{0}\), e.g., the trajectory for \(h=0.5\) circulates around 0 and hence never reach 0 or 1.\vspace{0.2in}

\noindent
It is also shown in Figure \ref{fig:(^n)x_complex} (a) that the trajectories for \(h=-0.5\) is a semicircle and the distorted semicircles for \(h\) and \(-(1+h)\) are axially symmetric each other, e.g., -0.2 and -0.8. The relation between \({^h}x\) and \(^{-(1+h)}y\) is given as Lemma \ref{lem_(^n)x_symmetry}, which explains not only the symmetry but also existence of a pair function convertible each other if \(-1\leq{h}\leq{0}\).

\begin{lem}
\label{lem_(^n)x_symmetry}
Let \(-1\leq{h}\leq{0}\).\\
Let \(x>0\) and \(y>0\).\\
Let \(({^h}x)^{*}\) be the complex conjecture of \({^h}x\).\\
Then the following equation holds:
\[({^h}x)^{*}=1-({^{-(1+h)}}y),\]
or with the q-analog notation:
\[[h+1]_{\ln{x}}=1-[-h]_{\ln{y}},\]
where
\[\ln{y}=\frac{1}{\ln{x}}.\]
\end{lem}
\begin{proof}
Let \(q=\ln{x}\) and \(k>0\).\\
For \(0<x<1\), obviously \(q<0\), and \(q\) can be written as \(q=-k=ke^{i\pi}\). Hence we have
\[q^{h+1}=k^{h+1}[\cos\pi(h+1)+i\sin\pi(h+1)]\]
The tetrational function \({^h}x\) for \(-1\leq{h}\leq{0}\) is then expressed as
\[
{^h}x=[h+1]_q=\frac{q^{h+1}-1}{q-1}=\frac{h^{k+1}\cos\pi(h+1)}{-k-1}+i\frac{k^{h+1}\sin\pi(h+1)}{-k-1}.
\]
Therefore
\[
({^h}x)^{*}=-\frac{h^{k+1}\cos\pi(h+1)}{k+1}+i\frac{k^{h+1}\sin\pi(h+1)}{k+1}.
\]
On the other hand,
\[\ln{y}=\frac{1}{q}=k^{-1}e^{i\pi}.\]
Hence we get
\[
\left(\frac{1}{q}\right)^{-(h+1)}=k^{h+1}[\cos\pi(h+1)-i\sin\pi(h+1)].
\]
Then
\[
1-(^{-(1+h)}y)=1-[-h]_{1/q}=1-\frac{\left(\frac{1}{q}\right)^{-h}-1}{\frac{1}{q}-1}=1-\frac{\left(\frac{1}{q}\right)^{-(h+1)}-q}{1-q}
\]
\[
=1-\frac{k^{h+1}\cos\pi(h+1)+k}{1+k}+i\frac{k^{h+1}\sin\pi(h+1)}{1+k}
\]
\[
=-\frac{h^{k+1}\cos\pi(h+1)}{k+1}+i\frac{k^{h+1}\sin\pi(h+1)}{k+1}.
\]
Therefore \(({^h}x)^{*}=1-({^{-(1+h)}}y)\) holds.\vspace{0.2in}\\
\noindent
For \(x\geq{1}\), obviously \(q\geq{0}\), and we have
\[({^h}x)^{*}={^h}x=\frac{q^{h+1}-1}{q-1}.\]
On the other hand, we have
\[1-(^{-(1+h)}y)=1-\frac{\left(\frac{1}{q}\right)^{-h}-1}{\frac{1}{q}-1}=\frac{q^{h+1}-1}{q-1}.\]
Therefore generally \(({^h}x)^{*}=1-({^{-(1+h)}}y)\) holds.\\
\end{proof}\vspace{0.2in}

\noindent
Now let us consider the extension of the tetration to complex bases and heights.

\begin{thm}
\label{thm_complex_domains}
Let \(r\) be real constant.\\
Let \(z\) be a complex variable. \\
Then the following statements hold:
\begin{itemize}
\item \({^z}r\) is holomorphic if and only if \(\Re(z)\notin\mathbb{Z}\).
\item \({^r}z\) is holomorphic if and only if \(r\in\mathbb{Z}\).
\end{itemize}
\end{thm}
\begin{proof}
Let \(z=x+iy\).\\
(1) Proof of the first statement.\\
We shall prove this by induction.\\
Let \(s=\ln{q}=\ln{r}\).\\
Then \(q^{z+1}=e^{s(x+1+iy)}=e^{s(x+1)}[\cos(sy)+i\sin(sy)]\).\\
For \(-1\leq{x}<{0}\) we have
\[
{^x}r=\frac{q^{x+1+iy}-1}{q-1}=\frac{e^{s(x+1)}\cos(sy)-1}{q-1}+i\frac{e^{s(x+1)}\sin(sy)}{q-1}
=u_{0}+iv_{0},\]
where \(u_{0}\) and \(v_{0}\) are real and imaginary component of \({^x}r\), respectively.\\
Obviously the Cauchy-Riemann equations holds:
\[
\frac{\partial{u_{0}}}{\partial{x}}
=\frac{\partial{v_{0}}}{\partial{y}}
=\frac{se^{s(x+1)}\cos(sy)}{q-1},
\]
\[
\frac{\partial{u_{0}}}{\partial{y}}
=-\frac{\partial{v_{0}}}{\partial{v}}
=-\frac{se^{s(x+1)}\sin(sy)}{q-1}.
\]
Therefore \({^x}r\) is holomorphic in \(-1\leq{x}<{0}\).\\
Let \(n\in\mathbb{Z}\). Suppose \(^{(x+iy)}r=u_{n}+iv_{n}\) for \((n-1)\leq{x}<{n}\) is holomorphic.\\
Then
\[
^{(x+1+iy)}r=r^{(u_{n}+iv_{n})}=e^{qu_{n}}\cos(qv_{n})+ie^{qu_{n}}\sin(qv_{n})=u_{n+1}+iv_{n+1}.
\]
By using the chain rule, e.g.,
\[
\frac{\partial{u_{n+1}}}{\partial{x}}
=\frac{\partial{u_{n+1}}}{\partial{u_{n}}}
\frac{\partial{u_{n}}}{\partial{x}}
+\frac{\partial{u_{n+1}}}{\partial{v_{n}}}
\frac{\partial{v_{n}}}{\partial{x}},
\]
we can prove that Cauchy-Riemann equations holds:
\[
\frac{\partial{u_{n+1}}}{\partial{x}}
=\frac{\partial{v_{n+1}}}{\partial{y}}
=qe^{qu_{n}}\cos(qv_{n})\frac{\partial{u_{n}}}{\partial{x}}-qe^{qu_{n}}\sin(qv_{n})\frac{\partial{v_{n}}}{\partial{x}},
\]
\[
\frac{\partial{v_{n+1}}}{\partial{x}}
=-\frac{\partial{u_{n+1}}}{\partial{y}}
=qe^{qu_{n}}\sin(qv_{n})\frac{\partial{u_{n}}}{\partial{x}}-qe^{qu_{n}}\cos(qv_{n})\frac{\partial{v_{n}}}{\partial{x}}.
\]
Therefore \({^z}r\) is generally holomorphic in inside of each segment \((n-1)<x<n\).\\ However the Cauchy-Riemann equations do not hold at boundaries of segments \(x=n\) as follows:\\
The real and imaginary components of \({^z}r\) in segment \(-1\leq{x}<{0}\) approach to the following values as \(x\rightarrow{0}\):
\[
\underset{x\rightarrow{0}}{\lim}{u_{0}}=\frac{q\cos(sy)-1}{q-1},\hspace{20pt}\underset{x\rightarrow{0}}{\lim}{v_{0}}=\frac{q\sin(sy)}{q-1}.
\]
On the other hand, \({^z}r\) approaches the following values as \(x\rightarrow{0}\) since \(u_{1}+iv_{1}=e^{qu_{0}}\cos(qv_{0})+ie^{qu_{0}}\sin(qv_{0})\):
\[
\underset{x\rightarrow{0}}{\lim}{u_{1}}=\exp\left[q\frac{q\cos(sy)-1}{q-1}\right]\cos\left(q\frac{q\sin(sy)}{q-1}\right)
\]
\[
\underset{x\rightarrow{0}}{\lim}{v_{1}}=\exp\left[q\frac{q\cos(sy)-1}{q-1}\right]\sin\left(q\frac{q\sin(sy)}{q-1}\right)
\]
Obviously, in general \({^r}z\) is not continuous at the segment boundary. It is continuous if and only if the heights are real, \(y=0\) hence \(v_{1}=v_{0}=0\).\\
Similarly,
\[
\underset{x\rightarrow{0}}{\lim}{u_{n+1}}=e^{qu_{n}}\cos(qv_{n}),\hspace{20pt}\underset{x\rightarrow{0}}{\lim}{v_{n+1}}=e^{qu_{n}}\sin(qv_{n}).
\]
Then \(\underset{x\rightarrow{n}}{\lim}(u_{n+1}+iv_{n+1})=\underset{x\rightarrow{n}}{\lim}(u_{n}+iv_{n})\) holds if and only if \(v_{n+1}=v_{n}=0\).\\
Therefore, the first statement holds.\vspace{0.2in}

\noindent
(2) Proof of the second statement\\
Let \(h,k\in\mathbb{R}\) and let \(\ln{z}=e^{h+ik}=e^{h}\cos{k}+ie^{h}\sin{k}\).\\
Since \(z=x+iy=e^{e^{h}\cos{k}}\cos{(e^{h}\cos{k})}+ie^{e^{h}\cos{k}}\sin{(e^{h}\cos{k})}\), we can confirm that the Cauchy-Riemann equations hold and expressed as follows:
\[
\frac{\partial{h}}{\partial{x}}=\frac{\partial{k}}{\partial{y}}=\frac{1}{e^{(e^{h}\cos{k}+h)}\cos{(e^{h}\sin{k}+k)}},
\]
\[
\frac{\partial{k}}{\partial{x}}=-\frac{\partial{h}}{\partial{y}}=-\frac{1}{e^{(e^{h}\cos{k}+h)}\sin{(e^{h}\sin{k}+k)}}.
\]
Let \(-1\leq{r}<0\). Then we have
\[
{^r}z=\frac{(\ln{z})^{r+1}-1}{\ln{z}-1}=\frac{e^{(r+1)h}\cos{k(r+1)}+ie^{(r+1)h}\sin{k(r+1)}-1}{e^{h}\cos{k}+ie^{h}\sin{k}-1}
\]
\[
=\frac{[e^{(r+1)h}\cos{k(r+1)}+ie^{(r+1)h}\sin{k(r+1)}-1][(e^{h}\cos{k}-1)-ie^{h}\sin{k}]}{(e^{h}\cos{k}-1)^{2}+e^{2h}\sin^{2}{k}}
\]
\[
=\frac{e^{(r+2)h}\cos{rk}-e^{(r+1)h}\cos{k(r+1)}-e^{h}\cos{k}+1}{e^{2h}-2e^{h}\cos{k}+1}
\]
\[
+i\frac{e^{(r+2)h}\sin{rk}-e^{(r+1)h}\sin{k(r+1)}+e^{h}\sin{k}}{e^{2h}-2e^{h}\cos{k}+1}
\]
\[=u+iv,\]
where \(u\) and \(v\) are real and imaginary component of \({^r}z\), respectively.\\
The partial derivatives of \(u\) and \(v\) by \(h\) or \(k\) are:
\begin{align*}
\frac{\partial{u}}{\partial{h}}&=\frac{(r+2)e^{(r+2)h}\cos{kr}-(r+1)e^{(r+1)h}\cos{k(r+1)}-e^{h}\cos{k}}{e^{2h}-2e^{h}\cos{k}+1}\\
&-\frac{\left(e^{(r+2)h}\cos{rk}-e^{(r+1)h}\cos{k(r+1)}-e^{h}\cos{k}+1\right)(2e^{2h}-2e^{h}\cos{k})}{(e^{2h}-2e^{h}\cos{k}+1)^{2}},\\
\frac{\partial{u}}{\partial{k}}&=\frac{-re^{(r+2)h}\sin{kr}+(r+1)^{2}e^{(r+1)h}\sin{k(r+1)}+e^{h}\sin{k}}{e^{2h}-2e^{h}\cos{k}+1}\\
&-\frac{\left(e^{(r+2)h}\cos{rk}-e^{(r+1)h}\cos{k(r+1)}-e^{h}\cos{k}+1\right)(2e^{h}\sin{k})}{(e^{2h}-2e^{h}\cos{k}+1)^{2}},\\
\end{align*}
\begin{align*}
\frac{\partial{v}}{\partial{h}}&=\frac{(r+2)e^{(r+2)h}\sin{kr}-(r+1)e^{(r+1)h}\sin{k(r+1)}+e^{h}\sin{k}}{e^{2h}-2e^{h}\cos{k}+1}\\
&-\frac{\left(e^{(r+2)h}\sin{rk}-e^{(r+1)h}\sin{k(r+1)}+e^{h}\sin{k}\right)(2e^{2h}-2e^{h}\cos{k})}{(e^{2h}-2e^{h}\cos{k}+1)^{2}},\\
\frac{\partial{v}}{\partial{k}}&=\frac{re^{(r+2)h}\cos{kr}-(r+1)^{2}e^{(r+1)h}\cos{k(r+1)}+e^{h}\cos{k}}{e^{2h}-2e^{h}\cos{k}+1}\\
&-\frac{\left(e^{(r+2)h}\sin{rk}-e^{(r+1)h}\sin{k(r+1)}+e^{h}\sin{k}\right)(2e^{h}\sin{k})}{(e^{2h}-2e^{h}\cos{k}+1)^{2}}.
\end{align*}
It is obvious that the Cauchy-Riemann equations do not hold.
\[
\frac{\partial{u}}{\partial{x}}=\frac{\partial{u}}{\partial{h}}\frac{\partial{h}}{\partial{x}}+\frac{\partial{u}}{\partial{k}}\frac{\partial{k}}{\partial{x}}
\neq
\frac{\partial{v}}{\partial{h}}\frac{\partial{h}}{\partial{y}}+\frac{\partial{v}}{\partial{k}}\frac{\partial{k}}{\partial{y}}=\frac{\partial{v}}{\partial{y}},
\]
\[
\frac{\partial{u}}{\partial{y}}=\frac{\partial{u}}{\partial{h}}\frac{\partial{h}}{\partial{y}}+\frac{\partial{u}}{\partial{k}}\frac{\partial{k}}{\partial{y}}
\neq
-\frac{\partial{v}}{\partial{h}}\frac{\partial{h}}{\partial{x}}-\frac{\partial{v}}{\partial{k}}\frac{\partial{k}}{\partial{x}}=-\frac{\partial{v}}{\partial{x}}.
\]
Similarly the tetrational function in other segments \({^{(r+n)}}z=\exp_{z}^{n}({^r}z)\) are generally not holomorphic inside the segment \((n-1)<x<n\) \\
However, at the boundaries of the segment \(r=n\), the function is expressed as \({^n}z=\exp_{z}^{n}[1]\) and is holomorphic proved by induction as follows .\\
Let \(z=e^{x+iy}\) and \({^{(n+1)}}z=x_{n}+iy_{n}\).\\
For \(n=1\),
\[
{^1}z=e^{x+iy}=e^{x}\cos{y}+ie^{x}\sin{y}=u_{1}+iv_{1}
\]
The Cauchy-Riemann equations hold:
\[\frac{\partial{u}}{\partial{x}}=\frac{\partial{v}}{\partial{y}}=e^{x}\cos{y},\hspace{20pt}\frac{\partial{u}}{\partial{y}}=-\frac{\partial{v}}{\partial{x}}=-e^{x}\sin{y}.\]
Hence \({^1}z\) is holomorphic.\\
Suppose \({^n}z=u_{n-1}+iv_{n-1}\) is holomorphic:
\[\frac{\partial{u_{n-1}}}{\partial{x}}=\frac{\partial{v_{n-1}}}{\partial{y}},\hspace{20pt}\frac{\partial{u_{n-1}}}{\partial{y}}=-\frac{\partial{v_{n-1}}}{\partial{x}}.\]
Then \(^{(n+1)}z\) is expressed as:
\[^{(n+1)}z=z^{(^{n}z)}=e^{(x+iy)(u_{n-1}+iv_{n-1})}=e^{(xu_{n-1}-yv_{n-1})+i(xv_{n-1}+yu_{n-1})}\]
\[=e^{(xu_{n-1}-yv_{n-1})}\cos{(xv_{n-1}+yu_{n-1})}+ie^{(xu_{n-1}-yv_{n-1})}\sin{(xv_{n-1}+yu_{n-1})}\]
\[=u_{n}+v_{n}.\]
By using the relations above, it is shown that the Cauchy-Riemann equations hold:
\[
\frac{\partial{u_{n}}}{\partial{x}}=\frac{\partial{v_{n}}}{\partial{y}}
\]
\[
=\left(u_{n-1}+x\frac{\partial{u_{n-1}}}{\partial{x}}-y\frac{\partial{v_{n-1}}}{\partial{x}}\right)e^{(xu_{n-1}-yv_{n-1})}\cos{(xv_{n-1}+yu_{n-1})}
\]
\[
-\left(v_{n-1}+x\frac{\partial{v_{n-1}}}{\partial{x}}+y\frac{\partial{u_{n-1}}}{\partial{x}}\right)e^{(xu_{n-1}-yv_{n-1})}\sin{(xv_{n-1}+yu_{n-1})},
\]
\[
\frac{\partial{u_{n}}}{\partial{y}}=-\frac{\partial{v_{n}}}{\partial{x}}
\]
\[
=-\left(v_{n-1}+x\frac{\partial{v_{n-1}}}{\partial{x}}+y\frac{\partial{u_{n-1}}}{\partial{x}}\right)e^{(xu_{n-1}-yv_{n-1})}\cos{(xv_{n-1}+yu_{n-1})}
\]
\[
-\left(u_{n-1}+x\frac{\partial{u_{n-1}}}{\partial{x}}-y\frac{\partial{v_{n-1}}}{\partial{x}}\right)e^{(xu_{n-1}-yv_{n-1})}\sin{(xv_{n-1}+yu_{n-1})}.
\]
Hence \(^{(n+1)}z\) is holomorphic.\\
Therefore the second statement holds.
\end{proof}
\begin{cly}
If both the base and height are complex values without integer components, \(x,y\in\mathbb{C}\) and \(x,y\notin\mathbb{Z}\), then the tetrational function \({^y}x\) is not holomorphic.
\end{cly}

\section{Calculation rules of tetration}
\label{sec_rule}
In this section, we study calculation rules of tetration based on the explicit formulae and properties of the functions.\vspace{0.2in}

\noindent
First let us extend the exponentiation operator.

\begin{dfn}
\label{dfn_extended_exp}
Let \(r,x,y\in\mathbb{R}\).\\
Extended exponentiation operator \(\exp_{r}^{x}\) is defined as:
\[\exp_{r}^{x}1\colonequals\exp_{r}^{\lfloor{x}\rfloor}\left(r^{[\{x\}]_q}\right)={^x}r,\]
where the notations in Theorem \ref{thm_explicit} are used.
\end{dfn}

\begin{lem}
\label{lem_e^(x-a)(^a)x}
Let \(r,x,y\in\mathbb{R}\).
The following equation holds:
\[{^x}r=\exp_r^{x-y}\left({{^y}r}\right)\]
\end{lem}
\begin{proof}
Since \(x=\lfloor{x}\rfloor+\{x\}\), by using the relation \ref{cly_core}, we have the following relations:
\[
{^x}r=\exp_{r}^{x}1=\exp_{r}^{\lfloor{x}\rfloor+\{x\}}1=\exp_{r}^{\lfloor{x}\rfloor}\left({^{\{x\}}r}\right)
.\]
Then for any \(w\in\mathbb{R}\),
\[
\exp_{r}^{x}1=\exp_{r}^{\lfloor{x}\rfloor-w+\{x\}+w}1=\exp_{r}^{\lfloor{x}\rfloor-w}\left({^{\{x\}+w}r}\right)
.\]
By using \(y=\{x\}+w\) and hence \(\lfloor{x}\rfloor-w=\lfloor{x}\rfloor+\{x\}-y=x-y\), we get the relation.
\end{proof}

\begin{lem}
\label{lem_exp_summation}
Let \(r,x,y\in\mathbb{R}\).
The following equation holds:
\[\exp_{r}^{x}\exp_{r}^{y}1=\exp_{r}^{x+y}1={^{x+y}}r.\]
\end{lem}
\begin{proof}
From Lemma\ref{lem_e^(x-a)(^a)x},
\[{^x}r=\exp_{r}^{x-a}\left({^a}r\right)=\exp_{r}^{x-a}\exp_{r}^{a}1\]
By rewriting \(y=x-a\), we get the relation.
\end{proof}\vspace{0.2in}

\noindent
Next, we distinguish the exponentiation (right associative) from the power (left associative) for clarity.\vspace{0.2in}

\begin{dfn}
\label{dfn_vee_wedge}
Let \(\vee\) and \(\wedge\) be the symbols of the exponentiation (right associative) and the power (left associative), respectively.
\begin{align*}
{^x}r\vee{f} & \colonequals{\exp_r^{x}{f}},\\
f\wedge{{^x}r} & \colonequals{f}^{({^x}r)},
\end{align*}
where \(f\) is the operand.\\
\noindent
The left-hand side of \(\vee\) is not a number but the part of the exponentiation operator. The iterated tetration is given as:
\[{^x}({^y}r)\vee{f}\colonequals(\exp_{r}^{y})^{x}f=\exp_{r}^{xy}f.\]
\end{dfn}\vspace{0.2in}

\noindent
The notation \({^x}r\vee{^y}r\) refers to \(\exp_{r}^x({^y}x)\), where the operand \({^y}r\) is a number while \({^x}r\vee\) is not a number but the extended exponent operator, not referring to \(\exp_{({^x}r)}\). On the other hand, the notation \({^x}r\wedge{^y}r\) refers to \(({^x}r)^{({^y}r)}=\exp_{({^x}r)}({^y}r)\), where both \({^x}r\) and \({^y}r\) are numbers. Similarly, \({^x}({^y}r)\vee{1}\) is not \({^x}t\) with \(t={^y}r\vee{1}\).\vspace{0.2in}

\noindent
Some calculation rules are given as in Theorem \ref{thm_tet_rules}.

\begin{thm}
\label{thm_tet_rules}
Let \(r,t,x,y,z\in\mathbb{R}\).
The following relations hold:
\begin{eqnarray}
\setcounter{equation}{1}
{^x}r\vee{1}&=&\hspace{3pt}{^x}r.\\
-{^x}r&=&{^{-1}}r\vee\left(\frac{1}{^{x+1}r}\right)\\
\frac{1}{{^x}r}&=&r\vee({-^{x-1}r})\\
{^x}r+{^y}r&=&{^{-1}}r\vee({^{x+1}}r\cdot{^{y+1}}r).\\
{^x}r-{^y}r&=&{^{-1}}r\vee\left(\frac{{^{x+1}}r}{{^{y+1}}r}\right).\\
{^{x}}r\wedge{^y}r&=&{^{y+1}}r\wedge{^{x-1}r}.\\
{^x}r\cdot{^y}r&=&{^{-1}r\vee\left({^x}r\wedge{^y}r\right)}=r\vee\left({^{x-1}}r+{^{y-1}r}\right)\\
\frac{{^x}r}{{^y}r}&=&{^{-1}r}\vee({^{x+1}}r\wedge\frac{1}{{^y}r})=r\vee({^{x-1}}r-{^{y-1}}r).
\end{eqnarray}
\begin{eqnarray}
{^x}r\vee{^y}r&=&{^y}r\vee{^x}r\hspace{3pt}=\hspace{3pt}{^{x+y}}r.\\
{^y}({^x}r)&=&{^x}({^y}r)\hspace{3pt}=\hspace{3pt}{^{xy}}r.\\
{^z}({^x}r\vee{^y}r)&=&{^{z(z+y)}}r.\\
{^{\{x\}}}r&=&r\vee\frac{(\ln{r})^{\{x\}}-1}{\ln{r-1}}\\
{^x}r&=&{^{\lfloor{x}\rfloor}}r\vee{^{\{x\}}}r\\
{^{\{x\}+\{y\}-1}}r&=&{^{\{x\}-1}}r+(\ln{r})^{\{x\}}\cdot\left({^{\{y\}-1}}r\right),\hspace{10pt}\{x\}+\{y\}<1.\\
{^x}r&=&\left({^{\lfloor{x}\rfloor+1}}r\right)\vee\left(1-{^{-\{x\}}}t\right)^{*},\hspace{10pt}\ln{r}\cdot\ln{t}=1.\\
\frac{\partial({^x}r)}{\partial({^{\{x\}-1}r})}&=&(\ln{r})^{\lfloor{x}\rfloor+1}\prod_{k=0}^{\lfloor{x}\rfloor}({^k}}r\vee{^{\{x\}}r).
\end{eqnarray}

\end{thm}
\begin{proof}
The proofs are given as fellows respectively:\\
(1) Directly given from Definition \ref{dfn_extended_exp}. This is the relation between operator \({^x}r\vee\) and number \({^x}r\):
\[{^x}r=\exp_{r}^{x}1={^x}r\vee{1}.\]
(2) By Definition \ref{dfn_vee_wedge}, we have \(\log_{r}f={^{-1}}r\vee{f}\).
\[-{^x}r=-\log_{r}[{^{x+1}}r]=\log_{r}\left[\frac{1}{{^{x+1}}r}\right]={^{-1}}r\vee\left[\frac{1}{^{x+1}r}\right].\]
(3) By the exponent rule, we have
\[\frac{1}{{^x}r}=({^x}r)^{-1}=(r^{^{x-1}r})^{(-1)}=r^{-(^{x-1}r)}=r\vee(-^{x-1}r).\]
(4) By the logarithm rule,
\[{^x}r+{^y}r=\log_{r}({^{x+1}}r)+\log_{r}({^{y+1}}r)=\log_{r}({^{x+1}}r\cdot{^{y+1}}r)={^{-1}}r\vee({^{x+1}}r\cdot{^{y+1}}r).\]
(5) By the logarithm rule,
\[{^x}r-{^y}r=\log_{r}({^{x+1}}r)-\log_{r}({^{y+1}}r)=\log_{r}\left(\frac{{^{x+1}}r}{{^{y+1}}r}\right)={^{-1}}r\vee\left(\frac{{^{x+1}}r}{{^{y+1}}r}\right).\]
(6) By the exponent rule,
\[{^x}r\wedge{^y}r=({^x}r)^{({^y}r)}=r^{(^{x-1}r)({^y}r)}=r^{({^y}r)(^{x-1}r)}={^{y+1}}r\wedge{^{x-1}r}.\]
(7) By the logarithm rule,
\[{^x}r\cdot{^y}r={^x}r\cdot\log_{r}({^{y+1}}r)=\log_{r}\left[({^{y+1}}r)^{({^x}r)}\right]={^{-1}r}\vee({^{y+1}r}\wedge{^{x}r})={^{-1}r}\vee({^{x}r}\wedge{^{y}r}),\]
where the (6) is used at the last step. By the exponent rule,
\[{^x}r\cdot{^y}r=r^{({^{x-1}}r)}\cdot{r}^{({^{y-1}}r)}=r^{({^{x-1}}r+{^{y-1}}r)}=r\vee({^{x-1}}r+{^{y-1}}r).\]
(8) By the logarithm rule,
\[\frac{{^x}r}{{^y}r}=\frac{1}{{^y}r}\log_{r}({^{x+1}}r)=\log_{r}\left[({^{x+1}}r)^{\frac{1}{({^y}r)}}\right]={^{-1}r}\vee({^{x+1}}r\wedge\frac{1}{{^y}r}).\]
By the exponent rule,
\[\frac{{^x}r}{{^y}r}=r^{({^{x-1}}r)}\cdot{r}^{(-{^{y-1}}r)}=r^{({^{x-1}}r-{^{y-1}}r)}=r\vee({^{x-1}}r-{^{y-1}}r).\]
(9) Directly give from Lemma \ref{lem_e^(x-a)(^a)x} and \ref{lem_exp_summation}.\\
(10) From definition \ref{dfn_vee_wedge}, we have the relation by assining \(f=1\).\\
(11) \({^z}({^x}r\vee{^y}r)={^z}({^x}r\vee{^y}r)\vee{1}=(\exp_{r}^{x+y})^{z}1=\exp_{r}^{z(x+y)}1={^{z(x+y)}}r\vee{1}={^{z(x+y)}}r\).\\
(12) Corollary \ref{cly_core} directly give the relation.\\
(13) We have the relation from Definition \ref{dfn_extended_exp} and (12).\\
(14) This relation is the property of q-numbers:
\[
{^{\{x\}+\{y\}-1}}r=\frac{(\ln{r})^{\{x\}+\{y\}}-1}{\ln{r}-1}=\frac{(\ln{r})^{\{x\}}-1}{\ln{r}-1}+(\ln{x})^{\{x\}}\frac{(\ln{r})^{\{y\}}-1}{\ln{r}-1}
\]
\[={^{\{x\}-1}}r+(\ln{r})^{\{x\}}\cdot\left({^{\{y\}-1}}r\right).\]
(15) From Lemma \ref{lem_(^n)x_symmetry}, we have the following relation.
\[\left({^{\{x\}-1}}r\right)^{*}=1-{^{-\{x\}}}t\]
Then \({^{\{x\}-1}}r=\left(1-{^{-\{x\}}}t\right)^{*}\).\\
By applying the operation of \(\left({^{\lfloor{x}\rfloor+1}}r\vee\right)\), we have
\[{^{\lfloor{x}\rfloor}}r\vee{^{\{x\}}}r=\left({^{\lfloor{x}\rfloor+1}}r\right)\vee\left(1-{^{-\{x\}}}t\right)^{*}.\]
Then we use (13) and get the relation. \\
(16) We get the relation from Theorem \ref{thm_explicit} and Corollary \ref{cly_core}.\\
\end{proof}

\section{Concluding remarks}
\label{sec_summary}
We derived the simple explicit formulae of the continuous tetrational function as well as multi-tetrational funcion based only on the delay-differential equation. The definition of the multi-tetrational function and finding the Stirling numbers of the second kind in Taylor coefficients are the key steps of our approach. The solutions is piecewise connected function and is class \(C^{1}\) at connecting point for real heights, and the height is extended to complex heights except the segment boundaries. Our solution has advantages in simple expression as well as in continuity in real r-x plane. The series for each segment converges absolutely and has infinite radius of convergence . The solution should be compared with those of numerical approaches of class \(C^{\infty}\) \cite{Kouznetsov2009}\cite{Paulsen2017} with limited domain of base \(r>e^{1/e}\) and finite radius of convergence for Taylor expansion. The calculation rules of tetration are given by distinguishing exponentiation from power operation.

\begin{bibdiv}
 \begin{biblist}
 	\bib{Goodstein1947}{article}{
	title={Transfinite ordinals in recursive  number theory},
	volume={12},
	DOI={10.2307/2266486},
	journal={J. Symb. Log.},
	author={Goodstein, R. L.},
	year={1947},
	pages={123--129}
	}

   \bib{Scot2006}{article}{
	title={General relativity and quantum mechanics: toward a generalization of the Lambert W function},
	volume={17},
	DOI={10.1007/s00200-006-0196-1},
	journal={Appl. Algebra Eng. Commun. Comput.},
	author={Scot, T. C.},
	author={Mann, R.},
	author={Martinez II, R. E.},
	year={2006},
	pages={41--47}
	}

	\bib{Sun2017}{article}{
	title={On the topology of air navigation route systems},
	volume={170},
	DOI={10.1680/jtran.15.00106},
	journal={P I CIVIL ENG-ENERGY},
	author={Sun, X.},
	author={Wandelt, S.},
	author={Linke, F.},
	year={2017},
	pages={46--59}
	}

	\bib{Khetkeeree2019}{article}{
	title={Signal Reconstruction using Second Order Tetration Polynomial},
	DOI={10.1109/ITC-CSCC.2019.8793435},
	journal={2019 34th International Technical Conference on Circuits/Systems, Computers and Communications (ITC-CSCC), JeJu, Korea (South)},
	author={Khetkeeree, S.},
	author={Chansamorn, C.},
	year={2019},
	pages={1--4}
	}

	\bib{Furuya2019}{article}{
	title={Compaction of Church Numerals},
	VOLUME = {12},
	DOI={10.3390/a12080159},
	journal={Algorithms (Basel)},
	author={Furuya, I.},
	author={Kida, T.},
	year={2019},
	pages={159}
	}

    \bib{Euler1783}{article}{
	title={De serie Lambertina Plurimisque eius insignibus proprietatibus,},
	volume={2},
	journal={Acta Acad. Scient. Petropol.},
	author={Euler, L.},
	year={1783},
	pages={29--51}
	}

    \bib{Eisenstein1844}{article}{
	title={Entwicklung von \(\alpha^{\alpha^{\alpha^{\udots}}}\),},
	volume={28},
	doi = {10.1515/crll.1844.28.49},
	journal={J. Reine angew. Math.},
	author={Eisenstein, G.},
	year={1844},
	pages={49--52}
	}

    \bib{Corless1996}{article}{
	title={On the Lambert W function},
	volume={5},
	doi = {10.1007/BF02124750},
	journal={Adv. Comput. Math.},
	author={Corless, R. M.},
	author={Gonnet, G. H.},
	author={Hare, D. E. G.},
	author={Jeffrey, D. J.},
	author={Knuth, D. E.},
	year={1996},
	pages={329--359}
	}

    \bib{Hooshmand2006}{article}{
	title={Ultra power and ultra exponential functions},
	volume={17},
	doi = {10.1080/10652460500422247},
	journal={Integral Transforms Spec. Funct.},
	author={Hooshmand, M. H.},
	year={2006},
	pages={549--558}
	}

    \bib{Kouznetsov2009}{article}{
	title={Solution of \(F(z+1) = \exp(F(z))\) in Complex \(z\)-Plane},
	volume={78},
	doi = {10.1090/S0025-5718-09-02188-7},
	journal={Math. Comp.},
	author={Kouznetsov, D.},
	year={2009},
	pages={1647--1670}
	}

    \bib{Paulsen2017}{article}{
	title={Solving \(F(z+1)=b^{F(z)}\)  in the complex plane,},
	volume={43},
	doi = {10.1007/s10444-017-9524-1},
	journal={Adv. Comput. Math.},
	author={Paulsen, W.},
	author={Cowgill, S.},
	year={2009},
	pages={1261--1282}
	}

    \bib{Abramowitz1964}{book}{
	title={Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables},
	publisher={NationalBureau of Standards},
	isbn = {0-486-61272-4},
	author={Abramowitz, M.},
	author={Stegun, A.},
	year={1964},
	}

    \bib{Graham1994}{book}{
	title={Concrete Mathematics: A Foundation for Computer Science},
	publisher={Addison-Wesley Longman Publishing Co., Inc.},
	isbn={0201558025},
	edition={2nd},
	author={Graham, Ronald L.},
	author={Knuth, Donald E.},
	year={1994},
	}

 \end{biblist}
\end{bibdiv}

\end{document}