テトレーション

テトレーション（英: tetration）とは、冪乗の次となる4番目のハイパー演算である。つまり、自らの冪乗を指定された回数反復する二項演算である。超冪（ちょうべき）ともいう。テトレーションという語はルーベン・グッドスタインによって、「4」を意味する接頭辞 tetra- と「繰り返し」を意味する iteration から作り出された。

定義

任意の正の実数 a > 0 および非負整数 n ≥ 0 に対し、次のようにテトレーション ⁿa を再帰的に定める。

${}^{n}a:={\begin{cases}1&{\text{if }}n=0\\a^{\left[{}^{(n-1)}a\right]}&{\text{if }}n>0\end{cases}}$

冪乗の演算が右結合、すなわち 10¹⁰¹⁰ のように積みあがった指数の上側から計算していくように、テトレーションの計算も ⁿa に対する n の部分から計算していく。

定義から直ちに、次の等式が成り立つ。

$\log _{a}\left({}^{n+1}a\right)={}^{n}a$

a と 10 が互いに素であるとき、ⁿa の最後の d 桁がオイラーの定理から求められる。

表記

テトレーションには多数の表記が存在する。テトレーションに使われる表記の中にはペンテーションやヘキセーションなど、より高次のハイパー演算の表記にも使用できるものもいくつかある。

名称	表記	説明
ルーディ・ラッカーの表記	${}^{n}a$	マウラーとグッドスタインによって導入され、ルーディ・ラッカーの『無限と心』で広まった。
クヌースの矢印表記	$a\uparrow \uparrow n,a\uparrow ^{2}n$	矢印または添字を増やすことで拡張できる。
コンウェイのチェーン表記	$a\to n\to 2$	数字を増やす、またはチェーンを拡張することで拡張できる。
アッカーマン関数	${}^{n}2=\operatorname {Ack} (4,n-3)+3$	底が 2 のときに限り、アッカーマン関数による表記が可能。
指数関数の反復合成による表示	${}^{n}a=\exp _{a}^{n}(1)$	右辺の表記に関しては後述。
フーシュマンドの表記	$\operatorname {uxp} _{a}(n),a^{\underline {n}}$	フーシュマンドの論文では「ultra power」（超冪）と書かれている。
ハイパー演算子表記	$a[4]n,H_{4}(a,n)$ $\operatorname {hyper} (a,4,n),\operatorname {hyper4} (a,n)$	数字を増やすことで拡張でき、一連のハイパー演算子を与える。
ASCII表記	a^^n	ASCII文字で表現する際、冪乗をキャレット^で表すことから。
バウアーズの配列表記	$\lbrace a,b,2\rbrace$	拡張配列表記へと一般化でき、さらにBEAFおよびバードの配列表記へと一般化される。

反復指数関数

反復指数関数（英: iterated exponential function）、あるいは反復冪（英: iterated exponential）とは指数関数の反復合成、あるいはその類似の関数およびその値を指して呼ばれる関数である。以降で表記を簡単にするため、非負整数 n と正実数 a の2つのパラメータを持つ実関数 $\mathop {\mathrm {exp} } \nolimits _{a}^{n}(x)$ を次のように定義する：

$\exp _{a}^{n}(x)=a^{a^{\cdot ^{\cdot ^{a^{x}}}}}$ （n 個の a の上に x が乗っている）

この関数は他に次のような表記で書かれる。

名称	表記	説明
（指数の反復合成）	$\exp _{a}^{n}(x)$	指数関数の表記 $\exp _{a}(x)$ はオイラーによる。
クヌースの矢印表記	$\left(a{\uparrow }\right)^{n}(x)$	矢印の数を増やすことで拡張できる。en:Large numbersを参照。
ガリダキスの表記	${}^{n}(a,x)$	底の表記が小さくならない。
ASCII表記	exp_a^n(x)	標準的な表記をベースにASCII文字のみを使用した表記。
J言語表記	x^^:(n-1)x

例

以下の表では、大部分の値が指数表記による表記すら困難なほど巨大であるため、それらの表記には底を 10 とした反復指数関数を用いた。なお小数部を持つ値はすべて近似値である。

$x$	${}^{2}x$	${}^{3}x$	${}^{4}x$
1	1 (11)	1 (11)	1 (11)
2	4 (22)	16 (24)	65,536 (216)
3	27 (33)	7,625,597,484,987 (327)	1.258015 × 103,638,334,640,024
4	256 (44)	1.34078 ×10154 (4256)	$\exp _{10}^{3}(2.18726)$ (8.1 × 10153 桁)
5	3,125 (55)	1.91101 × 102,184 (53,125)	$\exp _{10}^{3}(3.33928)$ (1.3 × 102,184 桁)
6	46,656 (66)	2.65912 × 1036,305 (646,656)	$\exp _{10}^{3}(4.55997)$ (2.1 × 1036,305 桁)
7	823,543 (77)	3.75982 × 10695,974 (7823,543)	$\exp _{10}^{3}(5.84259)$ (3.2 × 10695,974 桁)
8	16,777,216 (88)	6.01452 × 1015,151,335	$\exp _{10}^{3}(7.18045)$ (5.4 × 1015,151,335 桁)
9	387,420,489 (99)	4.28125 × 10369,693,099	$\exp _{10}^{3}(8.56784)$ (4.1 × 10369,693,099 桁)
10	10,000,000,000 (1010)	1010,000,000,000	$\exp _{10}^{4}(1)$ (1.0 × 1010,000,000,000 桁)

微積分

テトレーション ⁿa に対する a は正の実数に対して定義できるので、n を固定したときにそれぞれ微分と積分が定義できる。

高さが定数の微分

任意の正の整数 n に対し、 ⁿx の微分は次のようになる。

${\frac {\mathrm {d} }{{\mathrm {d} }x}}{}^{n}x={\frac {1}{x}}\sum _{k=1}^{n}(\ln x)^{k-1}\prod _{j=n-k}^{n}{}^{j}x$

高さが定数の積分

区間 x ∈ (0, 1] における y = x^x と y = x^−x のグラフ。

1/²x、 ²x の 0 から 1 までの定積分は二年生の夢と呼ばれる。

${\begin{aligned}\int _{0}^{1}{\frac {1}{{}^{2}x}}\mathrm {d} x&=\quad \sum _{n=1}^{\infty }{\frac {1}{{}^{2}n}}&{\scriptstyle (=1.29128599706266354040728259059560054149861936827\dots )}\\\int _{0}^{1}{{}^{2}x}\,\mathrm {d} x&=-\sum _{n=1}^{\infty }{{}^{2}(-n)}&{\scriptstyle (=0.78343051071213440705926438652697546940768199014\dots )}\end{aligned}}$

任意の正の整数 n に対し、ⁿx の不定積分は次のようになる。

$\int {}^{n}x{\mathrm {d} }x=\sum _{j=0}^{n}{\frac {(-1)^{j}(j+1)^{j-1}}{j!}}b_{j}(x)+\sum _{j=n+1}^{\infty }(-1)^{j}a_{n,j}b_{j}(x)+C$

ここで a_j,k は

$a_{j,k}={\begin{cases}1&{\text{if }}k=0\\{\frac {1}{k!}}&{\text{if }}j=0\\{\frac {1}{k}}\sum _{l=1}^{k}la_{j,k-l}a_{j-1,l-1}&{\text{otherwise}}\end{cases}}$

で与えられる有理数であり、 b_j (x) は第2種不完全ガンマ関数を用いて

$b_{j}(x)=\Gamma \left(j+1,\ln x\right)$

で与えられる。

拡張

テトレーションは、高さが正の整数以外の場合に拡張できる。

底が0

0の0乗が単純には定義できないため、 ⁿ0 は直接定義できないが、極限が

$\lim _{x\to 0+}{}^{n}x={\begin{cases}1&n\in 2\mathbb {Z} \\0&n\in 2\mathbb {Z} +1\end{cases}}$

と収束するため、

${}^{n}0={\begin{cases}1&n\in 2\mathbb {Z} \\0&n\in 2\mathbb {Z} +1\end{cases}}$

と定義する。（ここで、2Z、2Z+1 はそれぞれ偶数、奇数の集合を表す。）

なお、ここで 0⁰ が一意に決まらないにもかかわらず ²0 (= 0⁰) が定義できるのは、 a^b の a と b が等しいという条件下で極限を取ったからである。

底が複素数

収束または振動する点

発散する点

複素数の累乗が可能なことから、テトレーションは複素数の底に対しても定義できる。

例えばテトレーション ⁿi は対数関数の主枝を用いて定められる。このときオイラーの公式から次の式が得られる。

$i^{a+bi}=e^{{\frac {1}{2}}\pi i(a+bi)}=e^{-{\frac {1}{2}}\pi b}\left(\cos {\frac {\pi a}{2}}+i\sin {\frac {\pi a}{2}}\right)$

従って任意の ⁿi = a + bi に対して ^{n + 1}i = a′ + b′i が次のように再帰的に定義できる。

${\begin{aligned}a'&=e^{-{\frac {1}{2}}{\pi b}}\cos {\frac {\pi a}{2}}\\b'&=e^{-{\frac {1}{2}}{\pi b}}\sin {\frac {\pi a}{2}}\end{aligned}}$

ここから以下の近似値が導かれる。

$^{n}i$	近似値
${}^{1}i=i$	$i$
${}^{2}i=i^{\left({}^{1}i\right)}$	$0.2079$
${}^{3}i=i^{\left({}^{2}i\right)}$	$0.9472+0.3208i$
${}^{4}i=i^{\left({}^{3}i\right)}$	$0.0501+0.6021i$
${}^{5}i=i^{\left({}^{4}i\right)}$	$0.3872+0.0305i$
${}^{6}i=i^{\left({}^{5}i\right)}$	$0.7823+0.5446i$
${}^{7}i=i^{\left({}^{6}i\right)}$	$0.1426+0.4005i$
${}^{8}i=i^{\left({}^{7}i\right)}$	$0.5198+0.1184i$
${}^{9}i=i^{\left({}^{8}i\right)}$	$0.5686+0.6051i$

同様に値を逆向きに求めていくことで ⁰i = 1、 ⁻¹i = 0 が得られる。 ⁿi の値を複素平面上にプロットすると、点列は渦巻状に極限値 0.4383 + 0.3606i へと近づく。この値は n→∞ のときと解釈できる。

このようなテトレーションの列に関する研究はオイラーの時代から続けられてきているものの、そのカオス的な振る舞いのために不明な所が多い。これまでに発表された研究のほとんどは、無限反復指数関数の収束について焦点を当てたものである。現在の研究は高性能のコンピュータを用いたフラクタルと数式処理システムの出現に大きな恩恵を受けている。テトレーションについて分かっていることの多くは、複素力学系の一般的な知識と、指数写像の専門的な研究によるものである。

高さが無限大

$e^{-e}\leq x\leq e^{1/e}$ における $\textstyle \lim _{n\to \infty }{}^{n}x$ のグラフ

関数 $\left|{\frac {\mathrm {W} (-\ln {z})}{-\ln {z}}}\right|$ の複素平面上のグラフ。実数値無限反復指数関数を黒い曲線で示した。

ある範囲の底 a に対して limn→∞ ⁿa は有限の値に収束するので、この範囲においてテトレーションは高さ無限大の場合へ拡張できる。例えば limn→∞ ⁿ(√2) は収束して、その値は 2 であるから、 ^∞(√2) = 2 であると言える。

一般に無限反復指数関数 $x^{x^{\cdot ^{\cdot ^{\cdot }}}}$ は、 n が無限大に向かうときの ⁿx の極限として定義される。これが e^−e ≤ x ≤ e^1/e （およそ 0.066 から 1.44）の範囲で収束することはオイラーによって示された。

極限値 ^∞x が存在するとき、これは方程式

${}^{\infty }x=x^{{}^{\infty }x}$

(1)

を満たす正の実数に等しい。

式(1)より x = (^∞x)^1/∞x であり、このとき右辺の最大値が e^1/e であることから、 x > e^1/e については極限値が存在しないことがわかる。また $x\mapsto {}^{\infty }x$ は $y\mapsto y^{1/y}$ の逆関数（の分岐の下部）であることがわかる。

式(1)から、極限値 ^∞z をランベルトのW関数 W を用いて次のように定義することで複素数の底 z に対しても拡張される。

${}^{\infty }z=z^{z^{\cdot ^{\cdot ^{\cdot }}}}={\frac {\mathrm {W} (-\ln {z})}{-\ln {z}}}$

(2)

高さが非正

定義より

${}^{n}a=\log _{a}\left\{{}^{(n+1)}a\right\}$

が成り立つので、この関係を k ≤ 0 に対しても帰納的に拡張し

${\begin{array}{rclclcccc}{}^{1}n&=&\log _{n}\left({}^{2}n\right)&=&\log _{n}\left(n^{n}\right)&=&n\log _{n}n&=&n\\{}^{0}n&=&\log _{n}\left({}^{1}n\right)&=&\log _{n}n&&&=&1\\{}^{-1}n&=&\log _{n}\left({}^{0}n\right)&=&\log _{n}1&&&=&0\end{array}}$

と定義する。

ただし、定義できるのは n = −1 までで、log 0 が存在しないため n = −2 に対しては定義できず、従って n ≤ −2 に対して拡張できない。

高さが実数

テトレーションを高さ実数または複素数へ拡張する、という一般的な問題への広く受け入れられた解答は今のところ存在しない。いくつかのアプローチについて以下で述べる。

一般にこの問題は任意の実数 a > 0 に対し、実数 x > −2 で定義され、次の条件を満たす超指数関数 f (x) = ^xa を求めるものである。

• f (0) = 1、 f (−1) = 0
• 任意の実数 x > −1 に対し f (x) = a↑↑(f (x−1))
• 三つ目の条件は通常次の中のどれかである。
  • 連続性（通常は x > 0 における a、 x についての連続性）
  • 微分可能性（x についての 1、 2、 k 回または無限回微分可能性）
  • 正則性（x についての 2 階微分可能性を含む）

  任意の実数 x > 0 に対して

  $\displaystyle {{\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}f(x)>0}$

  が成立する。

三つ目の条件は著者およびアプローチによって異なる。実数高さへの拡張には二つの主要なアプローチが存在し、一つは正則性、もう一つは微分可能性に基づいたものである。これらの二つのアプローチは、相反する結果を導くことから互いに大きく異なるとされ、調和は難しいと考えられている。

長さ 1 の区間で ^xa が定義されれば、任意の x >−2 に対し容易に拡張される。

一次近似

^xa の一次近似によるグラフ（a = 4, e, 2, 1.5, 0.5）。漸近線はx = -2

一次近似（連続性のもと微分可能性を近似）は次のように与えられる。

${}^{x}a\approx {\begin{cases}1+x&-1<x\leq 0\\a^{\left({}^{x-1}a\right)}&0<x\\\log _{a}({}^{x+1}a)&x\leq -1\end{cases}}$

ゆえに

近似	定義域
${}^{x}a\approx x+1$	$-1<x<0$
${}^{x}a\approx a^{x}$	$0<x<1$
${}^{x}a\approx a^{a^{x-1}}$	$1<x<2$
$\vdots$	$\vdots$

と以下続く。但しこの微分可能性はあくまで区分的なものである。整数 x を境に微分係数が ln a 倍されるため、（a = e の場合を除いて）整数 x において微分不可能となる。

以下は値の計算例である。

${\begin{aligned}{}^{\pi /2}e&\approx 5.868...\\{}^{-4.3}{0.5}&\approx 4.03335...\\{}^{-5.264}{0.6}&\approx -5.35997...\\{}^{3.1}{0.7}&\approx 0.7580...\end{aligned}}$

Ultra exponential function

フーシュマンドは ultra exponential function という関数を導入した。 これはテトレーションの一次近似を表し、 uxp_a (x) と表記される。 uxp_a は、次の定理によって一意に定められる。

定理 ―  $0<a\neq 1$ に対し、連続な関数 $f\colon (-2,+\infty )\to {\mathbb {R} }$ が条件

• 任意の実数 $x>-1$ に対し $f(x)=a^{f(x-1)}$
• $f(0)=1$
• $f$ は $(-1,0)$ で微分可能
• $f'$ は $(-1,0)$ で広義単調増加か広義単調減少
• $\lim _{x\to 0+}f'(x)=(\ln a)\lim _{x\to 0-}f'(x)$

を満たすとき、$f$ は次の方程式によって一意に定まる。

$f(x)=\exp _{a}^{[x]}\left(a^{\lbrace {x}\rbrace }\right)=\exp _{a}^{[x+1]}\left(\lbrace {x}\rbrace \right)$

但し $[x]$ はガウス記号であり、また $\lbrace {x}\rbrace =x-[x]$ である。

証明は、三番目から五番目の条件より f が [−1, 0] で線型となることから従う。

フーシュマンドはさらに次のような一意性定理を導いた。

定理 ―  連続な関数 $f\colon (-2,+\infty )\to {\mathbb {R} }$ が条件

• 任意の $x>-1$ に対し $f(x)=e^{f(x-1)}$
• $f(0)=1$
• $f$ は $(-1,0)$ で下に凸
• $\lim _{x\to 0-}f'(x)\leq \lim _{x\to 0+}f'(x)$

を満たすとき、$f$ は $\operatorname {uxp} (=\operatorname {uxp} _{e})$ に等しい。

証明は先とほぼ同様である。漸化式より $\lim _{x\to -1+0}f'(x)=\lim _{x\to +0}f'(x)$ となること、三番目の条件より f が (−1, 0) で線型となることから従う。

定理より、 x >−1 に対し f (x) = exp (f (x−1)) および f (0) = 1 であって、かつ (−1, +∞) で下に凸であるような関数 f は唯一 uxp のみである。f が十分微分可能であるためには (−1, 0) で極値を持つ必要がある。

より高次の近似

^xa の二次近似によるグラフ（a = 4, e, 2, 1.5, 0.5）

（微分可能性についての）二次近似は次のように与えられる。

${}^{x}a\approx {\begin{cases}\log _{a}({}^{x+1}a)&x\leq -1\\1+{\frac {2\ln(a)}{1+\ln(a)}}x-{\frac {1-\ln(a)}{1+\ln(a)}}x^{2}&-1<x\leq 0\\a^{\left({}^{x-1}a\right)}&0<x\end{cases}}$

これは任意の x > 0 について微分可能であるが二階微分可能でない。a = e のときこれは一次近似に等しくなる。

三次近似および高次への一般化は次のように与えられる。

高さが複素数

複素平面上にテトレーション $f=F(x+iy)$を解析接続したものを描画。 $|f|=1,e^{\pm 1},e^{\pm 2},\ldots$と$\arg(f)=0,\pm 1,\pm 2,\ldots$を太い曲線で示した。

次の条件を満たす関数 F が一意に定まる事が証明されている 。

• F (z + 1) = exp (F (x))
• F (0) = 1
• z→±i∞ のとき F (z) が対数関数の不動点（およそ 0.318 ± 1.337i ）に近づく
• 実数 z < −2 を除く複素平面全域で正則

この関数 F を右図に示す。また、底が e ではない場合についても、底が$e^{1/e}$よりも大きい場合については同様に証明されている。倍精度浮動小数点数近似はオンラインで公開されている。

一意性

テトレーションを一意に定めるためには正則性の条件が重要となる。いま、関数 F に対し関数 S を次のように構成する。

${\begin{aligned}A(z)&=\sum _{n=1}^{\infty }\sin(2\pi nz)\alpha _{n}\\B(z)&=\sum _{n=1}^{\infty }{\bigl (}1-\cos(2\pi nz){\bigr )}\beta _{n}\\S(z)&=F\left(z+A(z)+B(z)\right)\end{aligned}}$

ここで α_n、 β_n は十分速く減衰する実数列であり、少なくとも実軸の近くで A (z)、 B (z) を収束させるとする。

この関数 S は F と同様に最初の二つの条件 S (z + 1) = exp (S (z))、 S (0) = 1 を満たす。また α_n、 β_n が十分速く 0 に近づくとき、S は正の実軸近傍で解析的となる。しかし α_n、 β_n が全て 0 でない場合、S は新たに大量の特異点と不連続線を複素平面上に持つことになる。これは sin (z)、 cos (z) が虚軸に沿って指数関数的に増大するためである。これらの特異点は α_n、 β_n が小さければ小さいほど実軸から離れていくため、 S が正則であるためには全ての α_n、 β_n が 0となる、即ち S = F であればよい。

実解析上のテトレーションは一意的に定まらないので、複素平面への拡張は一意性に必要である。

未解決問題

• ⁿπ, ⁿe が整数になるような正の整数 n は存在するか。特に ⁴π(≈9.080222455390617769723931713×10^{666262452970848503}) は整数か。
• 与えられた自然数$n$と$q\in \mathbb {Q} \cap \left(\left(0,\infty \right)\setminus \mathbb {Z} \right)$に対し、$^{n}q$は整数か。 特に ⁴x = 2 の正の解 x は有理数か。

逆関係

冪は冪根と対数の二つの逆関係を持つ。これに倣って、以下テトレーションの逆関係をそれぞれ超冪根（英: super-root）と超対数（英語: super-logarithm）と呼ぶ。

超冪根

超冪根はテトレーションの底に関する逆関係である。

超平方根

y = √x_s のグラフ。

超平方根（英: super square root）は ²x の逆であり、二つの等価な表記 ssrt(x), √x_s を持つ。

この関数は次のようなランベルトのW関数による表示を持つ。

$\mathrm {ssrt} (x)=e^{W(\mathrm {ln} (x))}={\frac {\mathrm {ln} (x)}{W(\mathrm {ln} (x))}}$

またこの関数により、冪根と対数の間の鏡映的な関係が表れる。次の方程式は y = ssrt x のときに真となる。

${\sqrt[{y}]{x}}=\log _{y}x$

平方根と同様に超平方根は一つとは限らない。ただし、平方根と異なり超平方根の個数を決定するのは容易とは言えない。一般に、e^−1/e < x < 1 のとき x は二つの正の超平方根を 0 と 1 の間に持ち、1 < x のとき x は 1 より大きい一つの正の超平方根を持ち、0 < x < e^−1/e のとき x は超平方根を実数の範囲で持たない。しかし上の式より、任意の x (≠ 1) は可算無限個の超平方根を複素数の範囲で持つことが従う。

超平方根はネットワークのクラスタサイズを決定するのに使用される。

その他の超冪根

任意の整数 n > 2 に対して ⁿx は定義され、x ≥ 1 のとき増加となり ⁿ1 = 1 を満たす。従って x ≥ 1 のとき ⁿ√x_s は存在する。しかし上述した一次近似を用いた場合、−1 < y ≤ 0 のとき ^yx は x によらず y + 1 となり、従ってこの場合 x = ^y√y + 1_s は存在しない。

超平方根のほか n 次の超冪根も同様の記号を用いて ⁿ√x_s と表すことができる。

超冪根は高さが無限大の場合へと拡張することができ、これは 1/e ≤ x ≤ e の場合に限り問題なく定義される（#高さが無限大を参照）。^∞x が存在するとき ^∞x = x^∞x が成り立つことから、無限次の超冪根は初等関数によって ^∞√x_s = x^1/x と表すことができる。例えば ^∞√2_s = 2^1/2 = √2 となる。

n を任意の正の整数とするとゲルフォント＝シュナイダーの定理より、超平方根 √n_s は整数または超越数となり、超立方根 ³√n_s は整数または無理数（これが超越数かどうかは知られていない）となる。

超対数

超対数はテトレーションの高さに関する逆関係である。

テトレーション ^xa を x に関して連続的に増加するものとして定義すると、任意の実数 x に対し超対数 slog_a x (a > 1) が定義される。

この関数 slog_a x は以下の式を満たす。

${\begin{array}{lcl}\operatorname {slog} _{a}{^{x}a}&=&x\\\operatorname {slog} _{a}a^{x}&=&1+\operatorname {slog} _{a}x\\\operatorname {slog} _{a}x&=&1+\operatorname {slog} _{a}\log _{a}x\\\operatorname {slog} _{a}x&>&-2\end{array}}$

さらなる拡張

テトレーションは$a\uparrow ^{2}b\uparrow ^{2}\cdot \cdot \cdot \uparrow ^{2}z$という風に拡張できる。そしてテトレーションの回数を数え上げるペンテーションを定義することができ、$a\uparrow ^{3}b$と表せる。

同じようにヘキセーションも定義できる($a\uparrow ^{4}b$)。この拡張を一般化してクヌースの矢印表記ができる。

またテトレーションはハイパーE表記で$a\uparrow \uparrow x=E(a)a\#x$と書ける。

関連項目

• 巨大数
• アッカーマン関数
• 二重指数関数
• ハイパー演算
• ペンテーション、ヘキセーション
• 反復対数
• en:Symmetric level-index arithmetic

外部リンク

• Daniel Geisler's site on tetration
• Tetration Forum
• Tetration - TORI - Mizugadro, the research site by Dmitrii Kouznetsov
• Gottfried Helms' site on tetration