Extension of Algebraic Solutions Using the Lambert W Function Sidney Edwards

Abstract

The Lambert W function has utility for solving various exponential and logarithmic equations arranged in the form of f(x)eg(x). This document presents a variety of categorized inversion formulas and identities making use of the W function and tetration. Related techniques are then used to derive polar forms of exponential and related functions.

Contents

1 Introduction 1 1.1 Overview Of The W Function 1 1.2 Tetration 2 1.3 Identities of the W function 2 2 New Solutions to Equations 3 2.1 Examples of Applications 6 2.2 Combining Branches of the W Function 6 3 Representation in polar coordinates 7 3.1 Rotations in Cartesian coordinates 8 4 Derivations 11 5 References 28 arXiv:1902.08910v3 [math.GM] 26 Oct 2020

[email protected] Extension of Algebraic Solutions Using the Lambert W Function 1

1 Introduction

This document presents a series of algebraic applications of the Lambert W function, specifically in using it for deriving new solutions to various exponential and logarithmic equations. Section 2, beginning on page 3 summarizes these solutions with proof of derivations starting on page 11.

Although the W function has been defined since the time of Euler, it has only recently been applied to common sciences and mathematics. It is used in models for mechanics [5], electrical engineering [6], quantum physics [7] and computational science [8]. For example, consider an equation used in chemical engineering for modeling the time evolution of film thickness of particles in a reaction, with D(t) representing film thickness over time, a and b representing reaction constants and W(t) is the Lambert W function [10]: {(t, D(t)) ∈ R2}

b h 2 2 i D(t) = 1 + W (−e−1−a t/b ) . a

To understand how to arrive at an exact solution for time, one must know the W function’s definition and properties in the following sections. The exact solution is shown on page 6.

1.1 Overview Of The W Function

The Lambert W function (or “W function”) denoted as W (z) (for any z ∈ C), is defined as the inverse of zez [1] and satisfies the relationship z = W (z)eW (z). It was originally introduced by Johann Heinrich Lambert when considering functions containing combinations of polynomials and exponents. Euler later utilized it for his famous power tower function and defined its relation to tetration [9]. Although it is analytically defined in the complex plane, this document restricts the W function and consequent derivations to real-valued functions.

1 Figure 1: The W function of the real variable x is shown below as multivalued for x ∈ R over [− e , 0) and is comprised of two branches [2]. The principal branch, denoted W0(x), is the Lambert W function on 1 1 the interval [− e , ∞) with W (x) ≥ −1, while the second branch, W−1(x), lies on the interval [− e , 0) with W (x) ≤ −1 [1]. Extension of Algebraic Solutions Using the Lambert W Function 2

1.2 Tetration The Lambert W function was shown to be related to tetration by Leonard Euler in his work on iterated exponentials [3]. Tetration is a term coined for the next iterative operator beyond exponentiation or otherwise called iterated exponentiation [4]. A notation similar to an exponent can be used to represent this operation:

a... na = aa | {z } n times with a being the base and n and being the height of tetration [8]. This can more simply be suggested as “a to the height of n”.

Tetration has two different inverse operations, one being an iterated logarithm [8] and other other being what is only colloquially called the super root or tetra root. The tetra root is the inverse operator that returns the base for a given nth height. For instance, in the case of xx, or 2x, the tetra square root would be used to return x, x tr2(x ) = x. This relates to the W function which can be used to define the tetra square root, derived from finding the inverse of f(x) = xx:

y = xx ln(y)) = x ln(x) = eln(x) ln(x) W (ln(y) = ln(x) eW (ln(y) = x.

This shows that the inverse of xx is eW (ln(f(x)) explicitly notates the tetra square root. For real numbers, − 1 this function’s principal branch is defined over the interval [e e , ∞).

In working with intricate exponents, it is useful to condense the notation for certain operations. One of these is the nth tetra root, which can be denoted generally as trn(x). For clarification, t represents tetra, r as root and n the height of tetration it inverts.

1.3 Identities of the W function Starting from W (z) z = W (z)e , z ∈ C [1] the following identities can be derived:

z W Exponential Identity z = W (z)eW (z) → eW (z) = W (z)

  W Log-Difference Identity ln(z) = ln W (z)eW (z) → ln(W (z) = ln(z) − W (z)

   nzn  W Product Identity nW (z) = W nW (z)enW (z) → nW (z) = W W (z)n−1

W Sum Identity Extension of Algebraic Solutions Using the Lambert W Function 3

    1 1  W (x) + W (y) = W [W (x) + W (y)]eW (x)+W (y) → W (x) + W (y) = W xy + W (y) W (x)

A change of base formula can be derived by generalizing the base of the exponent:

x logb(a)x y = xa logb(a)y = logb(a)xb

W0,a(y) = x W0,b(logb(a)y) = logb(a)x

logb(a)x W (log (a)y) y = xb 0,b b = x logb(a)

W0,b(logb(a)y) and thus the change of base formula is W0,a(y) = . logb(a)

2 New Solutions to Equations

Below is the series of various exponential and logarithmic equations with their respective inverses alluded to earlier, arranged by category. Inverses of functions similar in form to zez harness the W function, starting with the basic example y = xeax. This function is inverted by manipulating the factor into the same form as the exponent to implore the use of the W function:

y = xeax W (ay) = ax ay = axeax W (ay) = x. a

This is also true for more complicated functions containing combinations of logarithms, exponents and specific types of polynomials (See pages 11-27 for proof of derivations).

In the following sections, generalized real constants are represented as a, b, c, d, f and g though with f(x) specifically representing a function of the real variable x in equations (4), (5) and (8). “e” in all derivations denotes Euler’s number. In application of any following solutions, note the previously mentioned domain restrictions. Only one branch of the W function will result if the argument of the inverse function is strictly greater than 0.

Products of Exponentials and Polynomials   c dx+f c d 1 db − f b y = (ax + b) e W y c e ac c − =x (1) d ac a

 1/c  c [a + c]y1/c  c W [a+c]y y = eaW (x)+bxc W e a+c ceb/c =x (2) a + c ceb/c 1  cd(y−f)d/b  d d bW ( ) y = axbecx + f bad/b =x (3)  cd 

bd/a dfx/a dfx/a ! ax + b W ( cde xe ) − W ( dfxe ) f(x) = f 0 a −1 a − bd/a =x (4) ce−dx − f d

−b −b −b −b ! W ( xe a ) − W ( xe a ) f(x) = ax coth(bx) − ax f 0 a −1 a =x. (5) 2b Extension of Algebraic Solutions Using the Lambert W Function 4

Products of Logarithms and Polynomials

 −f/b 1/c  ln(axb)c W −f[a dy] y = e a1/bbc =x (6) dxf 1− b+c  1  b b c b+c 1− b+c b hy i b+c y = ax W (x) ( ) b W =x (7) b b + c a  2 (d−cb/a)x  W ((d−cb/a)xe(d−cb/a)x)−W (− c(d−cb/a) xe ) ln(ax + b) e a − b f(x) = f =x (8) cx + d  a 

Products of Logarithms and Exponentials

x 2 y y = W (x)e ln( tr3(e )) =x (9)

Products of Logarithms

y ln(x)  1 − y  y−1 y = =x (10) W (x) W (1 − y)   1 eab 1 − eab y−c aW (bx) W ( y−c ) y = + c =x (11) x(W (bx) + 1) eb

Sums of Exponentials and Polynomials    dx 1 a dc d y−b y = ax + b + ce ln W e a =x (12) d dc a

Sums of Logarithms and Polynomials

y − b  c y = ax + b + ceW (x) W ec/a − =W (x) (13) a a ad y − ab x W ( e c c ) y = a(x + b) + c ln c =x (14) d a 1/a bW (eay/b) y = xa + ln(xb) =x (15) a

Sums of Logarithms

1 b " a+b by/a # aW ( a+b eby/a) aW ( e ) a y = ln(xa) + W (xb) a e b(a+b) =x (16) a + b y  y  y = W (ax) − W (bx) exp =x (17) abey − a2 bey − a y = W (x) + W (1/x) (18) ey/2 h i1/2 √ y2ey − 2 + p(y2ey − 2)2 − 4 2 √ √ −1/2 −y/2h 2 y 2 y 2 i ·e− 2e y e −2+ (y e −2) −4 = x

e−y/2 h i−1/2 √ y2ey − 2 − p(y2ey − 2)2 − 4 2 Extension of Algebraic Solutions Using the Lambert W Function 5

−y/2 h √ i−1/2 e √ y2ey −2− (y2ey −2)2−4 ·e 2 = x y = 2W (ax) − W ([bx]2) (19) r ! a b2 − e−y 1 − yey − 1 b2 a2

2 q 2  − a e−y 1− b yey −1 ·e b2 a2 = x

r ! a b2 e−y 1 − yey + 1 b2 a2

2 q 2  a e−y 1− b yey +1 ·e b2 a2 = x Pn   Qn ak 1/ k=1 ak  [ k=1 bk ] y n exp W Pn a X k=1 k y = x ak ln(bkx) Pn =x (20) Qn ak 1/ k=1 ak k=1 [ k=1 bk ]

Composites of Exponentials and Polynomials   c dx+f c d 1 db − f b y = (ax + b) e W y c e ac c − =x (21) d ac a 1  cd(y−f)d/b  d d bW ( ) y = axbecx + f bad/b =x (22)  cd 

d/b W ( a d ln(y))/d d e bc y = (axb)cx =x (23) a1/b

Composites of Exponentials and Logarithms

 1/c  c [a + c]y1/c  c W [a+c]y y = eaW (x)+bxc W e a+c ceb/c =x (24) a + c ceb/c √ − ln(b)± ln2(b)−4 ln(c)(ln(y)−ln(a)) y = a(bxln(c))ln(x) e 2 ln(c) =x (25)

y c y c 1 bsr (xc) √1 [W (2 ln(( ) b )) ln(( ) a ] 2 y = ax 2 e c 2 a a =x (26)

  d/b   exp W a d ln(y) /d d bc y = (axb)cx =x (27) a1/b 2   y b/c tr3 a y = a tr (bx)cx =x (28) 2 b Composites of Logarithms and Polynomials

ay 1/a bW (e b ) y = xa + ln(xb) =x (29) a " #1/b ! aW ( a+b eby/a) aW ( a+b eby/a) y = ln(xa) + W (xb) a exp a =x (30) a + b b(a + b) Extension of Algebraic Solutions Using the Lambert W Function 6

2.1 Examples of Applications Different physical and mathematical applications can be seen harnessing the W function. Starting with the diode article [6], it defines the model of the current and voltage in a solar cell it solves with the W function:

q V − IRs nK T (V −IRs) I = I0(e b − 1) + − Iph. Rsh

Here, I and V denote variable current and voltage respectively. To avoid confusion, I0 represents the reverse saturation current and Iph is the solar cell’s photocurrent. Let us denote these constants and variables more generally as this allows the equation to be more easily related to the above derivations: y − xc x + f + a = aeb[y−xc] + (31) d I[Rsh + Rs] + Rsh[Iph + I0]−  qR [I + I + I ] nKbT qI0Rsh sh ph 0 q W nK T exp = V b (nKbT )

Consider an equation for modeling film thickness in a chemical reaction [10]. The film thickness D(t) is represented as a function of time t with a characterizing the reaction constant in reaction a and b the diffusion in reaction b Solving for time using properties of the W function can be done as follows:

b h 2 2 i D(t) = 1 + W (−e−1−a t/b ) (32) a b2 a b2 a b2 a2 [ b D − 1)] + a2 ln( b D − 1) + a2 = t

The article on Cheillini integrability shows a solution to a differential model using the W function [5]. It specifies a general equation for the non-relativistic potential V (x) derived for a quadratically damped harmonic oscillator in Lagrangian mechanics: µ h i V ±(x ) = e−η± eG±(xi)G (x ) − eG±(x0)G (x ) . i 4 ± i ± 0

± indicates the sign of the quadratic damping term, xi represents the x coordinate of a turning point where the velocity of the particle is zero in phase space (indicating it models consecutive turning points), and G±(xi) is a substitution for the integral of the a polynomial term with two sub-cases. After making a substitution and simplifying, the article goes on to demonstrate the steps for explicitly calculating the recursive relation xi+1 itself using the W function:

−1 G±(xi) xi+1 = G± (W (e G±(xi))).

2.2 Combining Branches of the W Function One relatively new concept of interest is that the W function’s properties can allow both of its real branches to be used simultaneously as the identities of the W function remain invariant of the specific branch by its definition. Take for instance a function related to the generating function for Bernoulli numbers that was recently shown to be invertible using branch differences of the W function: x f(x) = [2]. e−x − 1 An algebraic method of solving this was derived by Jeffrey and Jankowski by breaking the equation into separate products using branch differences [2]. Another method of derivation can be found that exemplifies how both branches can be used from a functional Extension of Algebraic Solutions Using the Lambert W Function 7

perspective as in derivations (4), (5) and (8). One should consider however that, just as with branches of other functions such as the square root function, only one definition of the function may be assumed for a given equation.

Starting with x f(x) = , e−x − 1 we can substitute for a branch difference to observe how it transforms the function,

W0(x) − W−1(x) f(W0(x) − W−1(x)) = . e−W0(x)−W−1(x) − 1 Paying close attention to the exponent containing the branches, we can observe that it can be rewritten as

W0(x) − W−1(x) f(W0(x) − W−1(x)) = xW0(x) − 1 xW−1(x) which allows the x in the numerator and denominator to simplify. This fraction of branches in the denominator can be manipulated and added with the constant:

W0(x) − W−1(x) f(W0(x) − W−1(x)) = xW0(x) − 1 xW−1(x)

W0(x) − W−1(x) f(W0(x) − W−1(x)) = W0(x)−W−1(x) W−1(x)

f(W0(x) − W−1(x)) = W−1(x)

Finishing the last step requires having defined the W function in terms of the second branch to make use of the substitution xex,

x x x f(W0(xe ) − W−1(xe )) = W−1(xe ) x f(W0(xe ) − x) = x.

x From this process, one can conclude Jeffrey and Jankowski’s statement that the inverse function of e−x−1 is x ax −1 W0(xe ) − x. This function can be further generalized to several constants. If f(x) = e−bx−c , then f (x) is shown as ! W ( bx ebcx/a) − W ( bcx ebcx/a) f 0 a −1 a = x. (33) b

3 Representation in polar coordinates

Another instance where the W function arises is in representing Cartesian functions as polar functions. Below, one can see the conversion of several functions:

y = ln(x), y = ex, y = W (x), y = xex. Extension of Algebraic Solutions Using the Lambert W Function 8

Using conventional substitution, rewrite x and y in terms of r and θ to solve for r:

y = ln(x) y = ex r sin(θ) = ln(r cos(θ)) r sin(θ) = er cos(θ) er sin(θ) = r cos(θ) 1 = e−r cos(θ) er sin(θ) r sin(θ) = cos(θ) 1 r = re−r cos(θ) 1 sin(θ) re−r sin(θ) = cos(θ) cos(θ) − = −r cos(θ)e−r cos(θ) −1 sin(θ) −re−r sin(θ) = cos(θ) W (− cot(θ)) = −r cos(θ)

−r sin(θ) − sin(θ) W (− cot(θ)) −r sin(θ)e = = r cos(θ) − cos(θ) −r sin(θ) = W (− tan(θ) W (− tan(θ) r = . − sin(θ) y = W (x) r sin(θ) = W (r cos(θ)) r sin(θ)er sin(θ) = r cos(θ) y = xex sin(θ)er sin(θ) = cos(θ) r sin(θ) = r cos(θ)er cos(θ) r sin(θ) cos(θ) r cos(θ) e = sin(θ) = cos(θ)e sin(θ) sin(θ) = er cos(θ) r sin(θ) = ln(cot(θ)) cos(θ) ln(cot(θ)) ln(tan(θ)) r = = r sin(θ) cos(θ)

From these representations, it is easy to see one can use polar coordinates to transform a function such as ln(x) to its Cartesian inverse

π W (− tan(π − [θ + 2 ])) W (− cot(θ)) π = . − sin(π − [θ + 2 ]) − cos(θ)

3.1 Rotations in Cartesian coordinates Moving further with this notion, one can also use the polar form to solve for rotations of the curves created by these equations in Cartesian coordinates. Take for instance the equation y = ex, put into terms of r and θ as r sin(θ) = er cos(θ). Starting from this, a π rotation of the curve by 4 can be induced via the argument which then allows one to solve for the rotated curve in terms of Cartesian coordinates using the W function: π r sin(θ) = er cos(θ) → r sin(θ − ) = er cos(θ−π/4). 4 In Figure 2 on page 9 are the graphics for these implicit curves. The left image is of the unrotated exponential π function (in polar form) while the right image shows this same curve but rotated by 4 in polar coordinates. Using conventional trigonometric identities, the right-hand side can be expanded out

sin(θ) cos(θ) π r cos(θ− π ) sin(θ) cos(θ) r( √ + √ ) r sin(θ − ) = e 4 → r( √ − √ ) = e 2 2 4 2 2 Extension of Algebraic Solutions Using the Lambert W Function 9

3 3

2 2

1 1

-3 -2 -1 1 2 3 -3 -2 -1 1 2 3

-1 -1

-2 -2

-3 -3

Figure 2: The original exponential function y = ex shown in the left image is plotted in polar form in Mathematica using the preceding equations. The left image contains the same curve but rotated in polar π coordinates by 4 .

r sin(θ) r cos(θ) It is convenient here to make two substitutions of u = √ and v = √ , transforming the preceding 2 2 equation into the form of

u − v = eu+v.

We have seen a similar equation to this in (12) that is solved using the W function. A similar technique can be applied in this specific instance to solve for u and thus for its related Cartesian form y:

−v = eu+v − u −2v = eu+v − u − v 2v = −eu+v + u + v u+v e2v = e−e +u+v u+v e2v = eu+ve−e u+v −e2v = −eu+ve−e W (−e2v) = −eu+v ln(−W (−e2v)) = u + v ln(−W (−e2v)) − v = u     r cos(θ) r cos(θ) r sin(θ) ln −W − exp 2 √ − √ = √ 2 2 2 √ √ 2 ln(−W (− exp( 2x))) − x = y

The top-most graph in Figure 3.1 on page 10 shows the plots in polar form. Below shows the plots in their Cartesian form. We can see on the top of the two graphs that the original y = ex is plotted while on the π right the same curve is rotated by 4 and plotted in the form derived from above.

This rotation can be generalized to any angle ”φ” in a similar manner with the exception of angles of πn x r cos(θ) 2 . Starting with the original equation, y = e , it is represented in polar coordinates as r sin θ = e Extension of Algebraic Solutions Using the Lambert W Function 10

Figure 3: The two images show the standard exponential function in Cartesian form in the upper graph. π In the lower graph is a plot of the curve rotated by 4 expressed in Cartesian form. The left concave curve represents where the real component of both branches of the rotated function overlap.

nπ which can be rotated by an arbitrary angle φ (for φ 6= 2 ) and solved for y in Cartesian form: y = AeBx (34) → r sin(θ + φ) = AeBr cos(θ+φ) u = r cos(θ), a = cos(φ), v = r sin(θ), b = sin(φ), a2 + b2 a 1 ln( a W ( ABb exp(B u))) − u = v Bb ABb a a b csc(φ) cot(φ) B sec(φ)x B ln( AB W (AB tan(φ)e )) − cot(φ)x = y

Both branches of the W function are needed to describe the real rotation of curve excluding when the argument of W is strictly greater than or equal to zero. Almost the same exact process can be used to find Extension of Algebraic Solutions Using the Lambert W Function 11 the rotation for the natural logarithm function in Cartesian form as well: 1  x  y = ln (35) B A r cos(θ + φ) → sin(θ + φ) = 1 ln B A u = r cos(θ), a = cos(φ), v = r sin(θ), b = sin(φ),

2 2 1 b ABa Ba +Bb u b v = ln( W ( e b )) − u Ba ABa b a sec(φ) tan(φ) B csc(φ)x y = B ln( AB W (AB cot(φ)e )) − tan(φ)x

One may notice that the trigonometric functions within the rotations of both ex and ln(x) have recurring singularities for certain values of φ. The Cartesian rotation equation is invalid for interval rotations of πn πn φ = 2 . The proper procedure for φ = 2 should therefore be to derive the Cartesian rotation equation π separately. In the instance of φ = 2 , the resulting back-conversion to Cartesian form is y = − ln(x).

4 Derivations

The following pages contain the algebraic proofs of derivations in previous sections of this document. The first section of derivations corresponds to section 2 on page 3, arranged in the same order but without repeating identical derivations. Later derivations are proofs of the example applications in section 2 and polar representations in section 3.

When working with such formulas, it is important to care for the domain restrictions of the W function. x Since W (x) is not bijective, W (xe ) = x cannot always be assumed [1]. The principal branch W0(x) is real only over x ∈ [−1/e, ∞) and since W (−1e−1) = −1, it is necessary to specify the constraint that g(x) ≥ −1 g(x) g(x) for W0(g(x)e ) = g(x) to hold. Similarly, it is also necessary to that g(x) ≤ −1 for W−1(g(x)e ) = g(x) to hold.

 af  a − +b W y cd e d In derivation (1), the inverse function of y is proven as a . Since the domain of the principal a − af +b branch is real for x ∈ [−1/e, ∞), the argument y cd e d must also be greater than -1/e.

a y cd a af −b−1 cd acf−bcd−cd ≥ −1/e, y cd ≥ −e d , y ≥ (−1) a e a . af −b e d

 cd cd This reveals that one must require a ∈ Z | (∃z) a = z to avoid complex values in the argument. One acf−bcd−cb cd acf−bcd−cd must also constrain the exponentiated term a such that (∀ a ) = 2z + 1, a ≥ −1. In the cd case that a = 2z, the exponential has a positive coefficient and thus the exponentiated term can be taken cd as any real number that satisfies the first constraint. Similarly for W−1, a = 2z + 1 for the argument of W−1 to remain in [−1/e, 0). Extension of Algebraic Solutions Using the Lambert W Function 12

Derivation (1), (21) Derivation (2), (24)

y = (ax + b)cedx+f y = eaW (x)+bxc c 1/c dx + f   y = (ax + b)e c c = eaW (x)+b W (x)eW (x) 1/c −f/c dx/c y e = (ax + b)e = eaW (x)+becW (x)W (x)c d d y y1/ce−f/c = (ax + b)edx/c = eaW (x)ecW (x)W (x)c ac ac eb dx db aW (x)+cW (x) c = + edx/c = e W (x) c ac 1 y c [a+c]   c W (x) d 1/c db − f dx db dx + db b = e W (x) y e ac c = + e c ac e c ac c ac 1/c [a + c]y [a+c]   c W (x) a+c d 1/c db − f dx db = e W (x) c W y e ac c = + ceb/c ac c ac [a + c]y1/c    W = W (x) a+c d 1/c db − f db dx b/c c W y e ac c − = ce ac ac c  [a+c]y1/c    cW b c d 1/c db − f b ce c W y e ac c − = x = W (x) d ac a a + c    1/c  [a+c]y1/c [a+c]y cW   cW b b ce c ce c e a+c = x a + c

Derivation (3), (22)

d y = axbecx + f

d y − f = axbecx

y − f d = xbecx a 1/b (y − f) cxd = xe b a1/b d/b (y − f) cdxd d b d = x e a b d/b d cd(y − f) cdx cdxd b d = e ba b b cd(y − f)d/b  cdxd W d = ba b b  cd(y−f)d/b  bW d ba b = xd cd 1/d   cd(y−f)d/b  bW bad/b = x  cd  Extension of Algebraic Solutions Using the Lambert W Function 13

Derivation (4)

ax + b f(x) = . ce−dx − f W (x) − W (x) a [W (x) − W (x)] + b f 0 −1 = d 0 −1 . d ce−[W0(x)−W−1(x)] − f W (x) − W (x)  a [W (x) − W (x) − bd/a] + b f 0 −1 − bd/a = d 0 −1 . d ce−[W0(x)−W−1(x)] − f a [W (x) − W (x)] = d 0 −1 . ce−[W0(x)−W−1(x)]−bd/a − f a [W (x) − W (x)] = d 0 −1 . cebd/a W0(x) − f W−1(x) W (x) − W (fx)  a [W (x) − W (fx)] f 0 −1 − bd/a = d 0 −1 . d cfebd/a W0(x) − f W−1(fx) a [W (x) − W (fx)] = d 0 −1 f[cebd/a W0(x) − 1] W−1(fx)  bd/a  a bd/a W0(ce x) − W−1(fx) d [W0(ce x) − W−1(fx)] f − bd/a = bd/a d f[ W0(ce x) − 1] W−1(fx) a [W (cebd/ax) − W (fx)] = d 0 −1 h bd/a i f W0(ce x)−W−1(fx) W−1(fx) W (cebd/axefx) − W (fxefx)  a f 0 −1 − bd/a = x d d bd/a dfx/a dfx/a ! W ( cde xe ) − W ( dfxe ) f 0 a −1 a − bd/a = x d Extension of Algebraic Solutions Using the Lambert W Function 14

Derivation (5)

f(x) = ax coth(bx) + ax = ax (coth(bx) + 1) x ax f( ) = (coth(x/2) + 1) 2b 2b ax ex/2 + e−x/2 = ( + 1) 2b ex/2 − e−x/2 ax ex/2 + e−x/2 ex/2 − e−x/2  = + 2b ex/2 − e−x/2 ex/2 − e−x/2 ax  2ex/2  = 2b ex/2 − e−x/2 ax  ex/2  = b ex/2 − e−x/2 W (x) − W (x) a(W (x) − W (x))  e(W0(x)−W−1(x))/2  f 0 −1 = 0 −1 2b b (eW0(x)−W−1(x))/2 − e−(W0(x)+W−1(x))/2 q  W−1(x)  a(W0(x) − W−1(x)) W0(x) = q q  b W−1(x) − W0(x) W0(x) W−1(x)  q  W−1(x) a(W (x) − W (x)) W (x) = 0 −1  0  b  √W−1(x)−W0(x)  W0(x)W−1(x) q  W−1(x)  −a W (x) = 0 b  √ 1  W0(x)W−1(x) s −a W (x)W (x)2 = 0 −1 b W0 −a = W (x) b −1 −b −b ! −b a −b a   W0( a xe ) − W−1( a xe ) −a −b −b f = W xe a 2b b −1 a −a (−b) = x b a = x Extension of Algebraic Solutions Using the Lambert W Function 15

Derivation (6) Derivation (7)

ln(axb)c y = axbW (cxd)f y = dxf y = xbW (cxd)f ln(axb)c a dy = f hy i1/b x = xW (cxd)f/b b 1 ln(ax ) a (dy) c = d/b f hy i d d df/b x c = x W (cx ) b a 1 ln(ax ) (dy) c = hy id/b ln(x) f c = cxdW (cxd)df/b e c a 1 b −f ln(x) (dy) c = ln(ax )e c d W (cxd) d df = W (cx )e W (cx ) b 1 1 −f ln(x) d (dy) c = b ln(a b x)e c = eW (cx )W (cxd)df/b+1 1 d −f ln(a b ) 1 1 W (cx ) d (df+b)/b e c (dy) c 1 −f(ln(x)+ln(a b ) = e W (cx ) = ln(a b x)e c d/(df+b)  d  b hy i bW (cx ) d − f 1 1 c = exp W (cx ) (a b dy) c 1 −f ln(a b x) a df + b = ln(a b x)e c b bc hy id/(df+b) bW (cxd) bW (cxd) f 1 = exp − 1 b −f(a b dy) c 1 −f ln(a x) −f b c df + b a df + b df + b bc = c ln(a x)e  h id/(df+b) d − f 1 bc y bW (cx ) −f(a b dy) c 1 W = W ( ) = ln(a b x) df + b a df + b bc  d/(df+b) f h i − 1 df + b bc y d W ( −f(a b dy) c ) 1 W = W (cx ) e bc = a b x b df + b a

f − 1 !  d/(df+b) −f(a b dy) c df+b bc hy i W 1 b W df+b e a b bc = x a   hy id/(df+b) · exp df+b W bc = cxd b df+b a  hy id/(df+b) df+b W bc bc df+b a   hy id/(df+b) · exp df+b W bc = xd b df+b a 1   d  d hy i df+b df+b W bc bc df+b a 1    d  d hy i df+b · exp df+b W bc = x b df+b a Extension of Algebraic Solutions Using the Lambert W Function 16

Derivation (8)

ln(ax + b) f(x) = cx + d x − b ln(x) f( ) = a x−b
c a + d ln(x) = c cb a x + (d − a ) ln(x) = cb  c  (d − a ) a(d−cb/a) x + 1 ex − b x f( ) = a cb  cex  (d − a ) a(d−cb/a) + 1 eW (ax)−W (cx) − b W (ax) − W (cx) f( ) = a cb  ceW (ax)−W (cx)  (d − a ) a(d−cb/a) + 1 W (ax) − W (cx) = cb  W (cx)  (d − a ) (d−cb/a)W (ax) + 1 W ( ax )−W (cx) ! ax e d−cb/a − b W ( ) − W (cx) f = d−cb/a a cb  W (cx)  (d − ) ax + 1 a W ( d−cb/a ) W ( ax )−W (−cx) ! ax e d−cb/a − b W ( ) − W (−cx) f = d−cb/a −W (−cx)+W ( ax ) a cb  d−cb/a  (d − ) ax a W ( d−cb/a )  ax  W d−cb/a = cb d − a  2 (d−cb/a)x  W ((d−cb/a)xe(d−cb/a)x)−W (− c(d−cb/a) xe ) e a − b f = x  a  Extension of Algebraic Solutions Using the Lambert W Function 17

Derivation (9) Derivation (11)

y = W (x)ex aW (bx) y = + c W (x) x[W (bx) + 1] y = W (x)eW (x)e −W (bx) W (x) abe y W (x)eW (x)e y − c = e = e [W (bx) + 1] y W (x) W (bx) tr3(e ) = e 1 e [W (bx) + 1] y = ln(tr3(e )) = W (x) y − c ab 2 y e eW (bx)+1[W (bx) + 1] ln( tr3(e )) = x = y − c ab eab = eW (bx)+1[W (bx) + 1] y − c eab Derivation (10) W ( y−c ) = W (bx) + 1 W ( eab ) − 1 = W (bx) ln(x) y−c y =   eab W (x) eab [W ( y−c )−1] [W y−c − 1]e = bx u = W (x) eab [W ( eab )−1] [W ( ) − 1]e y−c ln(ueu) y−c y = = x u b eab W ( eab ) ln(u) + u [W ( ) − 1]e y−c y = y−c = x u eb ln(u) eab [W ( eab ) − 1] y−c y = + 1 y−c W ( eab ) u y−c = x eb [1 − 1 ] eab ln(u) W ( eab ) y−c y − 1 = y−c = x u eb v = ln(u) y − 1 = ve−v 1 − y = −ve−v W (1 − y) = −v −W (1 − y) = v −W (1 − y) = ln(u) −W (1 − y) = ln(W (x)) e−W (1−y) = W (x)

−W (1−y) e−W (1−y)ee = x y  1 − y  y−1 = x W (1 − y) Extension of Algebraic Solutions Using the Lambert W Function 18

Derivation (12) Derivation (14)

dx x y = ax + b + ce y = a(x + b) + c ln( d ) x c u = e y x a(x+b) e = c e y = a ln(u) + b + cud d y x a (x+b) e c = e c y − b = a ln(u) + cud d y d a x+ ab ey−b = ea ln(u)+cu de c = xe c c y d − ab a x ey−b = uaecu de c c = xe c y−b c ud ad y − ab a a x e a = ue a e c c = xe c c c d y−b d cd ud e a = u e a ad y − ab a W ( e c c ) = x cd d y−b cd d cd ud c c e a = u e a ad y − ab a a W ( e c c ) c = x cd d y−b dc d a W ( e a ) = u a a a cd d y−b d W ( e a ) = u cd a a cd d y−b 1 ( W ( e a )) d = u Derivation (15), (29) cd a a b a cd d y−b 1 x y = x + ln(x ) ( W ( e a )) d = e cd a a ey = xbex 1 a cd d y−b ln( W ( e a )) = x y xa d cd a e b = xe b a ay a ax e b = x e b a ay a a ax e b = x e b b Derivation (13) ay a a W (e b ) = x b W (x) y = ax + b + ce ay bW (e b ) = xa y − b = ax + ceW (x) a ay 1 W (x) W (x) a y − b = aW (x)e + ce bW (e b ) = x y − b = eW (x)(aW (x) + c) a c y − b = eW (x)(W (x) + )a a y − b c = eW (x)(W (x) + ) a a y − b c c W (x) c e a = e a e (W (x) + ) a a y − b c W (x)+ c c e a = e a (W (x) + ) a a y − b c c W ( e a ) = W (x) + a a y − b c c W ( e a ) − = W (x) a a y − b c c W ( e a ) − = W (x) a a Extension of Algebraic Solutions Using the Lambert W Function 19

Derivation (16), (30) Derivation (17)

y = ln(xa) + W (xb) y = W (ax) − W (bx)

a b ey = eln(x )+W (x ) u = ax, b y = W (u) − W ( b u) ey = xaeW (x ) a b y a W (xb) y W ( a u) e = x e e = bW (u) u = xb −a[W (u) − W ( b u)] ey − a = a . y a W (u) b e = u b e bW (u)  a/b −a[W (u) − W ( b u)] ey = W (u)eW (u) eW (u) ey − a = a b bW (u) y a/b (a/b+1)W (u) e = W (u) e y W (u) − W ( b u) = a by a+b W (u) ey − a b e a = W (u)e a b −a[W (u)−W ( a u)] bW (u) by a+b a+b a a+b a W (u) a e = a W (u)e = bW (u). by a+b a a+b y y W ( a e ) = a W (u) bey −a y e = u = ax a+b by be − a aW ( e a ) a y  y  = W (u) exp = x a + b abey − a2 bey − a by by a+b a+b a aW ( e a ) aW ( a e ) a e a+b = u a + b by by a+b a+b a aW ( e a ) aW ( a e ) a b e a+b = x a + b by by a+b a+b a aW ( e a ) aW ( a e ) 1 a [ ] b e b(a+b) = x a + b Extension of Algebraic Solutions Using the Lambert W Function 20

Derivation (18)

y = W (x) + W (1/x) 1 y2ey = [W (x)2 + 2W (x)W (1/x) + W (1/x)2] W (x)W (1/x) W (x) W (1/x) = + 2 + W (1/x) W (1/x) W (x) W (1/x) y2ey − 2 = + W (1/x) W (x)

1 W (x) At this step, note that the right-hand side is in the form of u + u if u = W (1/x) . The inverse operation of this can therefore be used to return u. h i  W (x) ±1 1 y2ey − 2 ± p(y2ey − 2)2 − 4 = . 2 W (1/x)

There are now two possible inverse relations that will yield slightly different results. Take case 1 to be the positive branch of the square root, then, h i W (x) 1 y2ey − 2 + p(y2ey − 2)2 − 4 = 2 W (1/x) h i 1 W (x) ey 1 y2ey − 2 + p(y2ey − 2)2 − 4 = 2 W (x)W (1/x) W (1/x) 1 = W (1/x)2 h i−1 2e−y y2ey − 2 + p(y2ey − 2)2 − 4 = W (1/x)2 √ h i−1/2 2e−y/2 y2ey − 2 + p(y2ey − 2)2 − 4 = W (1/x)

√ h √ i−1/2 √ h i−1/2 −y/2 2 y 2 y 2 2e−y/2 y2ey − 2 + p(y2ey − 2)2 − 4 e 2e y e −2+ (y e −2) −4 = 1/x ey/2 h i1/2  √ h i−1/2 √ y2ey − 2 + p(y2ey − 2)2 − 4 exp − 2e−y/2 y2ey − 2 + p(y2ey − 2)2 − 4 = x 2 Now to case 2 with the negative branch of the square root: h i W (1/x) 1 y2ey − 2 − p(y2ey − 2)2 − 4 = 2 W (x) h i 1 W (1/x) ey 1 y2ey − 2 − p(y2ey − 2)2 − 4 = 2 W (x)W (1/x) W (x) 1 = W (x)2 e−y h i−1 y2ey − 2 − p(y2ey − 2)2 − 4 = W (x)2 2 e−y/2 h i−1/2 √ y2ey − 2 − p(y2ey − 2)2 − 4 = W (x) 2 e−y/2 h i−1/2 e−y/2 h i−1/2 √ y2ey − 2 − p(y2ey − 2)2 − 4 exp √ y2ey − 2 − p(y2ey − 2)2 − 4 = x 2 2 Extension of Algebraic Solutions Using the Lambert W Function 21

Derivation (19)

y = 2W (ax) − W ([bx]2) u = ax y = W (u) − W ([bu/a]2)  u 2 W ([bu/a]2) ey = W (u) bu 2 [ a ] a2W ([bu/a]2) = b2W (u)2 b2 W ([bu/a]2) ey = a2 W (u)2 b2 W ([bu/a]2) yey = [2W (u) − W ([bu/a]2)] a2 W (u)2 b2 2W (u)W ([bu/a]2) − W ([bu/a]2)2 W (u)2 yey − 1 = − a2 W (u)2 W (u)2 b2 [W (u) − W ([bu/a]2]2 1 − yey = a2 W (u)2 r b2 W (u) − W ([bu/a]] ± 1 − yey = a2 W (u) Case 1: r b2 W (u) − W ([bu/a]2) 1 − yey = a2 W (u) r b2 −W ([bu/a]2) 1 − yey − 1 = a2 W (u) r ! a2 b2 W (u)2 W ([bu/a]2) − e−y 1 − yey − 1 = b2 a2 W ([bu/a]2) W (u) = W (u)   r ! 2 q 2 a2 b2 − a e−y 1− b yey −1 − e−y 1 − yey − 1 e b2 a2 = u b2 a2   r ! 2 q 2 a b2 − a e−y 1− b yey −1 − e−y 1 − yey − 1 e b2 a2 = x b2 a2

Case 2: r b2 W ([bu/a]2) − W (u) 1 − yey = a2 W (u) r b2 W ([bu/a]2) 1 − yey + 1 = a2 W (u)   r ! 2 q 2 a2 b2 a e−y 1− b yey +1 e−y 1 − yey + 1 e b2 a2 = u b2 a2   r ! 2 q 2 a b2 a e−y 1− b yey +1 e−y 1 − yey + 1 e b2 a2 = x b2 a2 Extension of Algebraic Solutions Using the Lambert W Function 22

Derivation (20)

n X y = akx ln(bkx) = a0x ln(b0x) + a1x ln(b1x) + ... + anx ln(bnx) k=0 n X ak ak = x ln(bk x ) = x[a0 ln(b0x) + a1 ln(b1x) + ... + an ln(bnx)] k=0 n ! y X ak ak a0 a0 a1 a1 an an e = exp x ln(bk x ) = exp (x[ln(b0 x ) + ln(b1 x ) + ... + ln(bn x )]) k=0 n k Y a x a0 a0 a1 a1 an an x = [bkx] = [b0 x · b1 x · ... · bn x ] k=0 " n ! #x Y Pn ak k=0 ak a0 a0 a1 a1 an an x = bk x = [b0 x · b1 x · ... · bn x ] k=0 n X u = ak, k=0 n Y ak v = bk , k=0 ey = [vxu]x ey/u = [v1/ux]x  1/u  v y 1/u exp = [v1/ux]v x u    v1/uy  exp W ln exp = v1/ux u   v1/uy  exp W u = x v1/u Pn   Qn ak 1/ k=0 ak  [ k=0 bk ] y exp W Pn k=0 ak Pn = x Qn ak 1/ k=1 ak [ k=0 bk ] Extension of Algebraic Solutions Using the Lambert W Function 23

Derivation (23), (27) Derivation (28)

b cxd cx y = (ax ) y = a tr2(bx) y 1 1 cxd cx y b = (a b x) = tr2(bx) a 1 1 xd b/c y bc = (a b x) hy i bx = tr2(bx) d d d a bc b d x y = (a x ) tr2(bx) tr2(bx) = tr2(bx) d d a b d d d a b xd y bc = (a b x ) hy ib/c d tr3 = tr2(bx)  a b d  d a bc b d sr2 y = a x   b/c 2 hy i tr3 = bx d a W ( a b d ln(y)) 1 d e bc = (a b x) 2   y b/c d tr3 a W ( a b d ln(y))/d 1 = x e bc = a b x b d W ( a b d ln(y))/d e bc 1 = x a b

Derivation (24)

y = a(bxln(c))ln(x) y = eln(a)[eln(b)eln(c) ln(x)]ln(x) 2 y = eln(a)[eln(b) ln(x)eln(c) ln (x)] 2 y = eln(c) ln (x)+ln(b) ln(x)+ln(a) ln(y) = ln(c) ln2(x) + ln(b) ln(x) + ln(a) 0 = ln(c) ln2(x) + ln(b) ln(x) + ln(a) − ln(y) q − ln(b) ± ln2(b) − 4 ln(c)[ln(y) − ln(a)] = ln(x) 2 ln(c)  q  − ln(b) ± ln2(b) − 4 ln(c)[ln(y) − ln(a)] exp = x  2 ln(c)  Extension of Algebraic Solutions Using the Lambert W Function 24

Derivation (26)

c y = axb tr2(x ) c y b tr2(x ) a = x c/b c y
c tr2(x ) a = x c/b c y
c ln(x) tr2(x ) a = e y
c/b c ln(x)eW (c ln(x)) a = e  y
c/b W (c ln(x)) ln a = c ln(x)e y c/b W (c ln(x)) W (c ln(x)) ln(( a ) ) = W (c ln(x))e e y c/b 2W (c ln(x)) 2 ln(( a ) ) = 2W (c ln(x))e y c/b W (2 ln(( a ) )) = 2W (c ln(x)) 1 W (2 ln(( y )c/b)) = W (c ln(x)) 2 a y 1 y c/b 1 W (2 ln(( )c/b)) W (2 ln(( ) ))e 2 a = c ln(x) 2 a 1 y c/b 1 y c 1 W (2 ln(( ) )) 2 (2 ln(( ) a )) 2 = ln(x) 2c a a   1 y c/b y c 1/2 exp √ (W (2 ln(( ) )) ln(( ) a ) = x 2c a a Extension of Algebraic Solutions Using the Lambert W Function 25

Derivation (31) y − xc x + f + a = aeb[y−xc] + d db[x + f + a] = abdeb[y−xc] + b[y − xc]

b[y−xc] edb[x+f+a] = eb[y−xc]eabde

b[y−xc] abdedb[x+f+a] = abdeb[y−xc]eabde   W abdedb[x+f+a] = abdeb[y−cx]   W abdedb[x+f+a] = abdeb[y−cx]   ln(W abdedb[x+f+a] ) = ln(abd) + b[y − cx]     ln abdedb[x+f+a] − W abdedb[x+f+a] =   ln(abd) + db[x + f + a] − W abdedb[x+f+a] =   db[x + f + a] − W abdedb[x+f+a] = b[y − cx]

1  db[x+f+a] d[x + f + a] − b W abde = y − cx

1  db[x+f+a] x[d + c] + d[f + a] − b W abde = y  qR [I + I + I ] nKbT qI0Rsh sh ph 0 I[Rsh + Rs] + Rsh[Iph + I0] − q W nK T exp = V b (nKbT )

Derivation (32)

b h 2 2 i D(t) = 1 + W (−e−1−a t/b ) a aD 2 2 − 1 = W (−e−1−a t/b ) b a 2 2 a [ b D−1] −1−a t/b [ b D − 1]e = −e  a  a2t ln −[ a D − 1]e[ b D−1] = −1 − b b2 a2t −[ a D − 1] + ln( a D − 1) + 1 = − b b b2 b2 a b2 a b2 a2 [ b D − 1)] + a2 ln( b D − 1) + a2 = t Extension of Algebraic Solutions Using the Lambert W Function 26

Derivation (34)

y = AeBx → r sin(θ + φ) = AeBr cos(θ+φ) r sin(θ) cos(φ) + r cos(θ) sin(φ) = AeB[r cos(θ) cos(φ)−r sin(θ) sin(φ)] u = r cos(θ), a = cos(φ), v = r sin(θ), b = sin(φ), av + bu = AeBau−Bbv bu = AeBau−Bbv − av b A Bau−Bbv a u = a e − v Bb2 ABb Bau−Bbv a u = a e − Bbv Bb2 ABb Bau−Bbv a u = a e − Bbv Bb2 ABb Bau−Bbv a u + Bau = a e + Bau − Bbv 2 2 B a +b u ABb eBau−Bbv +Bau−Bbv e a = e a

a2+b2 ABb Bau−Bbv ABb B a u ABb Bau−Bbv a e a e = a e e a2b2 ABb B a u ABb Bau−Bbv W ( a e ) = a e a2+b2 a ABb B a u Bau−Bbv ABb W ( a e ) = e a2+b2 a ABb B a u ln( ABb W ( a e )) = Bau − Bbv a2+b2 a ABb B a u ln( ABb W ( a e )) − Bau = Bbv 2 2 1 a ABb B a +b u a ln( W ( e a )) − u = v Bb ABb a b  2 2  1 cos(φ) AB sin(φ) B cos(φ) +sin(φ) r cos(θ) cos(φ) ln W ( e cos(φ) ) − r cos(θ) = r sin(θ) B sin(φ) AB sin(φ) cos(φ) sin(φ) csc(φ) cot(φ)  ln W (AB tan(φ)eB sec(φ)x) − cot(φ)x = y B AB Extension of Algebraic Solutions Using the Lambert W Function 27

Derivation (35) 1 x y = ln( ) B A 1 r cos(θ) r sin(θ) = ln( ) B A 1 r cos(θ + φ) → r sin(θ + φ) = ln( ) B A 1 r cos(θ) cos(φ) − r sin(θ) sin(φ) r sin(θ) cos(φ) + r cos(θ) sin(φ) = ln( ) B A u = r cos(θ), a = cos(φ), v = r sin(θ), b = sin(φ) 1 ua − vb va + ub = ln( ) B A av+bu au − bv 1 e = ( ) B A au − bv eBav+Bbu = ( ) A AeBav+Bbu = au − bv A a eBav+Bbu = u − v b b A a eBav+Bbu + v = u b b ABa Ba2 eBav+Bbu + Bav = u b b ABa Ba2 eBav+Bbu + Bav = u + Bbu b b ABa Ba2 + Bb2 eBav+Bbu + Bav + Bbu = u b b 2 2 ABa eBav+Bbu+Bav+Bbu Ba +Bb u e b = e b

2 2 ABa Bav+Bbu ABa eBav+Bbu ABa Ba +Bb u e e b = e b b b 2 2 ABa Bav+Bbu ABa Ba +Bb u e = W ( e b ) b b 2 2 Bav+Bbu b ABa Ba +Bb u e = W ( e b ) ABa b 2 2 b ABa Ba +Bb u Bav + Bbu = ln( W ( e b )) ABa b 2 2 b ABa Ba +Bb u Bav = ln( W ( e b )) − Bbu ABa b 2 2 1 b ABa Ba +Bb u b v = ln( W ( e b )) − u Ba ABa b a sec(φ) tan(φ) r sin(θ) = ln( W (AB cot(φ)eB csc(φ)r cos(θ))) − tan(φ)r cos(θ) B AB sec(φ) tan(φ)  y = ln W (AB cot(φ)eB csc(φ)x) − tan(φ)x B AB Extension of Algebraic Solutions Using the Lambert W Function 28

Acknowledgments The author is grateful to the people who supported the author’s passion for this paper. The author would like to publicly thank: David Dollhopf, engineer, for his unparalleled enthusiasm and personal support; David Jeffrey, professor of applied mathematics for his exceptional patience and personal correspondence; Suzanne Woll, chemical engineer, for her time in providing feedback and guidance; Alex Hanhart, professor of mathematics for his continuous guidance in the author’s pursuit of mathematics.

References

[1] Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; Knuth, D. E. (1996). “On the Lambert W function”. Advances in Computational Mathematics. 5: 329-359. doi:10.1007/BF02124750.

[2] Jeffrey, D. J.; Jankowski, J. E. (2014). “Branch Differences of Lambert W”. 16th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing, 16(9), 61-65. doi:10.1109/SYNASC.2014.16. [3] Dennis, M. R.; Glendinning, P.; Martin, P. A.; Santosa, F.; Tanner, J. (2015). The Princeton Companion to Mathematics (pp. 152-156). Princeton, New Jersey: Princeton University Press. ISBN: 1400874475, 9781400874477. [4] Goodstein, R. L. (1947). “Transfinite ordinals in recursive number theory”. Journal of Symbolic Logic. 12(4): 123–129. doi:10.2307/2266486. [5] Choudhury, A. G.; Guha, P.; Pandey, A (2017). “Cheillini integrability and quadratically damped oscillators”. International Journal of Non-Linear Mechanics, 92(6), 153-159. ISSN: 0020-7462. [6] Roberts, Ken; Valluri, S.R. (2015). “On Calculating the Current-Voltage Characteristic of Multi-Diode Models for Organic Solar Cells”. arXiv [eprint], 1(12), 1-21. arXiv:1601.02679 [physics.comp-ph]. [7] Ishkhanyan, A.M. (2016). “The Lambert-W step-potential - an exactly solvable confluent hypergeometric potential”. Physics Letters A, 380(2), 640-644. doi:10.1016/j.physleta.2015.12.004. [8] Hittmeir, M. (2017). “A reduction of integer factorization to modular tetration”. arXiv [eprint], 3(7), 1-18. arXiv:1707.04919 [math.nt]. [9] Lynch, P. (2017). “The Fractal Boundary of the Power Tower Function”. Proceedings of Recreational Mathematics Colloquium, 5(1), 127-138.

[10] Braun, Artur; Wokaun, A; Hermanns, H.-G. (2003). “Analytical solution to a growth problem with two moving boundaries”. Applied Mathematical Modeling, 27(1), 47-52, doi:10.1016/S0307-904X(02)00085-9.