Natural logarithm

The natural logarithm of a number is its logarithm to the base of the mathematical constant e, where e is an irrational and transcendental number approximately equal to 2.718281828459. The natural logarithm of x is generally written as ln x, log_e x, or sometimes, if the base e is implicit, simply log x.[1][2] Parentheses are sometimes added for clarity, giving ln(x), log_e(x), or log(x). This is done in particular when the argument to the logarithm is not a single symbol, to prevent ambiguity.

Graph of the natural logarithm function. The function slowly grows to positive infinity as x increases and slowly goes to negative infinity as x approaches 0 ("slowly" as compared to any power law of x); the y-axis is an asymptote.

The natural logarithm of x is the power to which e would have to be raised to equal x. For example, ln 7.5 is 2.0149..., because e^2.0149... = 7.5. The natural logarithm of e itself, ln e, is 1, because e¹ = e, while the natural logarithm of 1 is 0, since e⁰ = 1.

The natural logarithm can be defined for any positive real number a as the area under the curve y = 1/x from 1 to a (the area being taken as negative when a < 1). The simplicity of this definition, which is matched in many other formulas involving the natural logarithm, leads to the term "natural". The definition of the natural logarithm can be extended to give logarithm values for negative numbers and for all non-zero complex numbers, although this leads to a multi-valued function: see Complex logarithm.

The natural logarithm function, if considered as a real-valued function of a real variable, is the inverse function of the exponential function, leading to the identities:

${\displaystyle {\begin{aligned}e^{\ln x}&=x\qquad {\text{if }}x>0,\\\ln e^{x}&=x.\end{aligned}}}$

Like all logarithms, the natural logarithm maps multiplication into addition:

${\displaystyle \ln xy=\ln x+\ln y.}$

Logarithms can be defined for any positive base other than 1, not only e. However, logarithms in other bases differ only by a constant multiplier from the natural logarithm, and can be defined in terms of the latter. For instance, the base-2 logarithm (also called the binary logarithm) is equal to the natural logarithm divided by ln 2, the natural logarithm of 2. Logarithms are useful for solving equations in which the unknown appears as the exponent of some other quantity. For example, logarithms are used to solve for the half-life, decay constant, or unknown time in exponential decay problems. They are important in many branches of mathematics and the sciences and are used in finance to solve problems involving compound interest.

History

Main article: History of logarithms

The concept of the natural logarithm was worked out by Gregoire de Saint-Vincent and Alphonse Antonio de Sarasa before 1649.[3] Their work involved quadrature of the hyperbola with equation xy = 1 by determination of the area of hyperbolic sectors. Their solution generated the requisite "hyperbolic logarithm" function having properties now associated with the natural logarithm.

An early mention of the natural logarithm was by Nicholas Mercator in his work Logarithmotechnia published in 1668,[4] although the mathematics teacher John Speidell had already in 1619 compiled a table of what in fact were effectively natural logarithms.[5] It has been said that Speidell's logarithms were to the base e, but this is not entirely true due to complications with the values being expressed as integers.[5]^:152

Notational conventions

The notations ln x and log_e x both refer unambiguously to the natural logarithm of x, and log x without an explicit base may also refer to the natural logarithm. This usage is common in mathematics and some scientific contexts as well as in many programming languages.[nb 1] In some other contexts, especially chemistry, however, log x can be used to denote the common (base 10) logarithm. It may also refer to binary (base 2) logarithm in the context of computer science, particularly in the context of time complexity.

Definitions

ln a as the area of the shaded region under the curve f(x) = 1/x from 1 to a. If a is less than 1, the area taken to be negative.

The area under the hyperbola satisfies the logarithm rule. Here A(s,t) denotes the area under the hyperbola between s and t.

The natural logarithm can be defined in several equivalent ways. The natural logarithm of a positive, real number a may be defined as the area under the graph of the hyperbola with equation y = 1/x between x = 1 and x = a. This is the integral

${\displaystyle \ln a=\int _{1}^{a}{\frac {1}{x}}\,dx.}$

If a is less than 1, this area is considered to be negative.

This function is a logarithm because it satisfies the fundamental property of a logarithm:

${\displaystyle \ln(ab)=\ln a+\ln b.}$

This can be demonstrated by splitting the integral that defines ln ab into two parts and then making the variable substitution x = at (so dx = a dt) in the second part, as follows:

${\displaystyle {\begin{aligned}\ln ab=\int _{1}^{ab}{\frac {1}{x}}\,dx&=\int _{1}^{a}{\frac {1}{x}}\,dx+\int _{a}^{ab}{\frac {1}{x}}\,dx\\[5pt]&=\int _{1}^{a}{\frac {1}{x}}\,dx+\int _{1}^{b}{\frac {1}{at}}a\,dt\\[5pt]&=\int _{1}^{a}{\frac {1}{x}}\,dx+\int _{1}^{b}{\frac {1}{t}}\,dt\\[5pt]&=\ln a+\ln b.\end{aligned}}}$

In elementary terms, this is simply scaling by 1/a in the horizontal direction and by a in the vertical direction. Area does not change under this transformation, but the region between a and ab is reconfigured. Because the function a/(ax) is equal to the function 1/x, the resulting area is precisely ln b.

The number e can then be defined to be the unique real number a such that ln a = 1. Alternatively, if the exponential function, denoted e^x or exp x, has been defined first, say by using an infinite series, the natural logarithm may be defined as its inverse function. In other words, ln is that function such that ln(exp x) = x. Since the range of the exponential function is all positive real numbers, and since the exponential function is strictly increasing, this is well-defined for all positive x.

Properties

• ${\displaystyle \ln 1=0}$
• ${\displaystyle \ln e=1}$
• ${\displaystyle \ln(xy)=\ln x+\ln y\quad {\text{for }}\;x>0\;{\text{and }}\;y>0}$
• ${\displaystyle \ln(x^{y})=y\ln x\quad {\text{for }}\;x>0}$
• ${\displaystyle \ln x<\ln y\quad {\text{for }}\;0<x<y}$
• ${\displaystyle \lim _{x\to 0}{\frac {\ln(1+x)}{x}}=1}$
• ${\displaystyle \lim _{\alpha \to 0}{\frac {x^{\alpha }-1}{\alpha }}=\ln x\quad {\text{for }}\;x>0}$
• ${\displaystyle {\frac {x-1}{x}}\leq \ln x\leq x-1\quad {\text{for}}\quad x>0}$
• ${\displaystyle \ln {(1+x^{\alpha })}\leq \alpha x\quad {\text{for}}\quad x\geq 0\;{\text{and }}\;\alpha \geq 1}$

Proof

The statement is true for ${\displaystyle x=0}$, and we now show that ${\displaystyle {\frac {d}{dx}}\ln {(1+x^{\alpha })}\leq {\frac {d}{dx}}(\alpha x)}$ for all ${\displaystyle x}$, which completes the proof by the fundamental theorem of calculus. Hence, we want to show that ${\displaystyle {\frac {d}{dx}}\ln {(1+x^{\alpha })}={\frac {\alpha x^{\alpha -1}}{1+x^{\alpha }}}\leq \alpha ={\frac {d}{dx}}(\alpha x)}$ (Note that we have not yet proved that this statement is true.) If this is true, then by multiplying the middle statement by the positive quantity ${\displaystyle (1+x^{\alpha })/\alpha }$ and subtracting ${\displaystyle x^{\alpha }}$ we would obtain ${\displaystyle x^{\alpha -1}\leq x^{\alpha }+1}$ ${\displaystyle x^{\alpha -1}(1-x)\leq 1}$ This statement is trivially true for ${\displaystyle x\geq 1}$ since the left hand side is negative or zero. For ${\displaystyle 0\leq x<1}$ it is still true since both factors on the left are less than 1 (recall that ${\displaystyle \alpha \geq 1}$). Thus this last statement is true and by repeating our steps in reverse order we find that ${\displaystyle {\frac {d}{dx}}\ln {(1+x^{\alpha })}\leq {\frac {d}{dx}}(\alpha x)}$ for all ${\displaystyle x}$. This completes the proof. An alternate proof is to observe that ${\displaystyle (1+x^{\alpha })\leq (1+x)^{\alpha }}$ under the given conditions. This can be proved, e.g., by the norm inequalities. Taking logarithms and using ${\displaystyle \ln(1+x)\leq x}$ completes the proof.

Derivative

The derivative of the natural logarithm as a real-valued function on the positive reals is given by

${\displaystyle {\frac {d}{dx}}\ln x={\frac {1}{x}}.}$

How to establish this derivative of the natural logarithm depends on how it is defined firsthand. If the natural logarithm is defined as the integral

${\displaystyle \ln x=\int _{1}^{x}{\frac {1}{t}}\,dt,}$

then the derivative immediately follows from the first part of the fundamental theorem of calculus.

If the natural logarithm is defined as the inverse of the (natural) exponential function, then the derivative for x > 0 can be found by using the properties of the logarithm and a definition of the exponential function. From the definition of the number ${\displaystyle e=\lim _{u\to 0}(1+u)^{1/u},}$ the exponential function can be defined as ${\displaystyle e^{x}=\lim _{u\to 0}(1+u)^{x/u}=\lim _{h\to 0}(1+hx)^{1/h}}$, where ${\displaystyle u=hx,h=u/x.}$ The derivative can then be found from first principles.

${\displaystyle {\begin{aligned}{\frac {d}{dx}}\ln x&=\lim _{h\to 0}{\frac {\ln(x+h)-\ln x}{h}}\\&=\lim _{h\to 0}\left[{\frac {1}{h}}\ln \left({\frac {x+h}{x}}\right)\right]\\&=\lim _{h\to 0}\left[\ln \left(1+{\frac {h}{x}}\right)^{\frac {1}{h}}\right]\quad &&{\text{all above for logarithmic properties}}\\&=\ln \left[\lim _{h\to 0}\left(1+{\frac {h}{x}}\right)^{\frac {1}{h}}\right]\quad &&{\text{for continuity of the logarithm}}\\&=\ln e^{1/x}\quad &&{\text{for the definition of }}e^{x}=\lim _{h\to 0}(1+hx)^{1/h}\\&={\frac {1}{x}}\quad &&{\text{for the definition of the ln as inverse function.}}\end{aligned}}}$

Series

The Taylor polynomials for ln(1 + x) only provide accurate approximations in the range −1 < x ≤ 1. Beyond some x > 1, the Taylor polynomials of higher degree are increasingly worse approximations.

If ${\displaystyle \vert x-1\vert \leq 1{\text{ and }}x\neq 0,}$ then[6]

${\displaystyle {\begin{aligned}\ln x&=\int _{1}^{x}{\frac {1}{t}}\,dt=\int _{0}^{x-1}{\frac {1}{1+u}}\,du\\&=\int _{0}^{x-1}(1-u+u^{2}-u^{3}+\cdots )\,du\\&=(x-1)-{\frac {(x-1)^{2}}{2}}+{\frac {(x-1)^{3}}{3}}-{\frac {(x-1)^{4}}{4}}+\cdots \\&=\sum _{k=1}^{\infty }{\frac {(-1)^{k-1}(x-1)^{k}}{k}}.\end{aligned}}}$

This is the Taylor series for ln x around 1. A change of variables yields the Mercator series:

${\displaystyle \ln(1+x)=\sum _{k=1}^{\infty }{\frac {(-1)^{k-1}}{k}}x^{k}=x-{\frac {x^{2}}{2}}+{\frac {x^{3}}{3}}-\cdots ,}$

valid for |x| ≤ 1 and x ≠ −1.

Leonhard Euler,[7] disregarding ${\displaystyle x\neq -1}$, nevertheless applied this series to x = −1, in order to show that the harmonic series equals the (natural) logarithm of 1/(1 − 1), that is the logarithm of infinity. Nowadays, more formally, one can prove that the harmonic series truncated at N is close to the logarithm of N, when N is large, with the difference converging to the Euler–Mascheroni constant.

At right is a picture of ln(1 + x) and some of its Taylor polynomials around 0. These approximations converge to the function only in the region −1 < x ≤ 1; outside of this region the higher-degree Taylor polynomials evolve to worse approximations for the function.

A useful special case for positive integers n, taking ${\displaystyle x={\tfrac {1}{n}}}$, is:

${\displaystyle \ln \left({\frac {n+1}{n}}\right)=\sum _{k=1}^{\infty }{\frac {(-1)^{k-1}}{kn^{k}}}={\frac {1}{n}}-{\frac {1}{2n^{2}}}+{\frac {1}{3n^{3}}}-{\frac {1}{4n^{4}}}+\cdots }$

If ${\displaystyle \operatorname {Re} (x)\geq 1/2,}$ then

${\displaystyle {\begin{aligned}\ln(x)&=-\ln \left({\frac {1}{x}}\right)=-\sum _{k=1}^{\infty }{\frac {(-1)^{k-1}({\frac {1}{x}}-1)^{k}}{k}}=\sum _{k=1}^{\infty }{\frac {(x-1)^{k}}{kx^{k}}}\\&={\frac {x-1}{x}}+{\frac {(x-1)^{2}}{2x^{2}}}+{\frac {(x-1)^{3}}{3x^{3}}}+{\frac {(x-1)^{4}}{4x^{4}}}+\cdots \end{aligned}}}$

Now taking ${\displaystyle x={\tfrac {n+1}{n}}}$ for positive integers n, yields:

${\displaystyle \ln \left({\frac {n+1}{n}}\right)=\sum _{k=1}^{\infty }{\frac {1}{k(n+1)^{k}}}={\frac {1}{n+1}}+{\frac {1}{2(n+1)^{2}}}+{\frac {1}{3(n+1)^{3}}}+{\frac {1}{4(n+1)^{4}}}+\cdots }$

If ${\displaystyle \operatorname {Re} (x)\geq 0{\text{ and }}x\neq 0,}$ then

${\displaystyle \ln(x)=\ln \left({\frac {2x}{2}}\right)=\ln \left({\frac {1+{\frac {x-1}{x+1}}}{1-{\frac {x-1}{x+1}}}}\right)=\ln \left(1+{\frac {x-1}{x+1}}\right)-\ln \left(1-{\frac {x-1}{x+1}}\right).}$

Since

${\displaystyle {\begin{aligned}\ln(1+y)-\ln(1-y)&=\sum _{i=1}^{\infty }{\frac {1}{i}}\left((-1)^{i-1}y^{i}-(-1)^{i-1}(-y)^{i}\right)=\sum _{i=1}^{\infty }{\frac {y^{i}}{i}}\left((-1)^{i-1}+1\right)\\&=y\sum _{i=1}^{\infty }{\frac {y^{i-1}}{i}}\left((-1)^{i-1}+1\right){\overset {i-1\to 2k}{=}}\;2y\sum _{k=0}^{\infty }{\frac {y^{2k}}{2k+1}},\end{aligned}}}$

we arrive at

${\displaystyle {\begin{aligned}\ln(x)&={\frac {2(x-1)}{x+1}}\sum _{k=0}^{\infty }{\frac {1}{2k+1}}{\left({\frac {(x-1)^{2}}{(x+1)^{2}}}\right)}^{k}\\&={\frac {2(x-1)}{x+1}}\left({\frac {1}{1}}+{\frac {1}{3}}{\frac {(x-1)^{2}}{(x+1)^{2}}}+{\frac {1}{5}}{\left({\frac {(x-1)^{2}}{(x+1)^{2}}}\right)}^{2}+\cdots \right).\end{aligned}}}$

Substituting again ${\displaystyle x={\tfrac {n+1}{n}}}$ for positive integers n, yields:

${\displaystyle {\begin{aligned}\ln \left({\frac {n+1}{n}}\right)&={\frac {2}{2n+1}}\sum _{k=0}^{\infty }{\frac {1}{(2k+1)((2n+1)^{2})^{k}}}\\&=2\left({\frac {1}{2n+1}}+{\frac {1}{3(2n+1)^{3}}}+{\frac {1}{5(2n+1)^{5}}}+\cdots \right).\end{aligned}}}$

This is, by far, the fastest converging of the series described here.

The natural logarithm in integration

The natural logarithm allows simple integration of functions of the form g(x) = f '(x)/f(x): an antiderivative of g(x) is given by ln(|f(x)|). This is the case because of the chain rule and the following fact:

${\displaystyle {\frac {d}{dx}}\ln \left|x\right|={\frac {1}{x}}.}$

In other words,

${\displaystyle \int {\frac {1}{x}}\,dx=\ln |x|+C}$

and

${\displaystyle \int {{\frac {f'(x)}{f(x)}}\,dx}=\ln |f(x)|+C.}$

Here is an example in the case of g(x) = tan(x):

${\displaystyle {\begin{aligned}&\int \tan x\,dx=\int {\frac {\sin x}{\cos x}}\,dx\\[6pt]&\int \tan x\,dx=\int {\frac {-{\frac {d}{dx}}\cos x}{\cos x}}\,dx.\end{aligned}}}$

Letting f(x) = cos(x):

${\displaystyle \int \tan x\,dx=-\ln \left|\cos x\right|+C}$

${\displaystyle \int \tan x\,dx=\ln \left|\sec x\right|+C}$

where C is an arbitrary constant of integration.

The natural logarithm can be integrated using integration by parts:

${\displaystyle \int \ln x\,dx=x\ln x-x+C.}$

Let:

${\displaystyle u=\ln x\Rightarrow du={\frac {dx}{x}}}$

${\displaystyle dv=dx\Rightarrow v=x}$

then:

${\displaystyle {\begin{aligned}\int \ln x\,dx&=x\ln x-\int {\frac {x}{x}}\,dx\\&=x\ln x-\int 1\,dx\\&=x\ln x-x+C\end{aligned}}}$

Numerical value

For ln(x) where x > 1, the closer the value of x is to 1, the faster the rate of convergence. The identities associated with the logarithm can be leveraged to exploit this:

${\displaystyle {\begin{aligned}\ln 123.456&=\ln(1.23456\cdot 10^{2})\\&=\ln 1.23456+\ln(10^{2})\\&=\ln 1.23456+2\ln 10\\&\approx \ln 1.23456+2\cdot 2.3025851.\end{aligned}}}$

Such techniques were used before calculators, by referring to numerical tables and performing manipulations such as those above.

Natural logarithm of 10

The natural logarithm of 10, which has the decimal expansion 2.30258509...,[8] plays a role for example in the computation of natural logarithms of numbers represented in scientific notation, as a mantissa multiplied by a power of 10:

${\displaystyle \ln(a\cdot 10^{n})=\ln a+n\ln 10.}$

This means that one can effectively calculate the logarithms of numbers with very large or very small magnitude using the logarithms of a relatively small set of decimals in the range ${\displaystyle [1,10)}$.

High precision

To compute the natural logarithm with many digits of precision, the Taylor series approach is not efficient since the convergence is slow. Especially if x is near 1, a good alternative is to use Halley's method or Newton's method to invert the exponential function, because the series of the exponential function converges more quickly. For finding the value of y to give exp(y) − x = 0 using Halley's method, or equivalently to give exp(y/2) − x exp(−y/2) = 0 using Newton's method, the iteration simplifies to

${\displaystyle y_{n+1}=y_{n}+2\cdot {\frac {x-\exp(y_{n})}{x+\exp(y_{n})}}}$

which has cubic convergence to ln(x).

Another alternative for extremely high precision calculation is the formula[9] [10]

${\displaystyle \ln x\approx {\frac {\pi }{2M(1,4/s)}}-m\ln 2,}$

where M denotes the arithmetic-geometric mean of 1 and 4/s, and

${\displaystyle s=x2^{m}>2^{p/2},}$

with m chosen so that p bits of precision is attained. (For most purposes, the value of 8 for m is sufficient.) In fact, if this method is used, Newton inversion of the natural logarithm may conversely be used to calculate the exponential function efficiently. (The constants ln 2 and π can be pre-computed to the desired precision using any of several known quickly converging series.)

Based on a proposal by William Kahan and first implemented in the Hewlett-Packard HP-41C calculator in 1979 (referred to under "LN1" in the display, only), some calculators, operating systems (for example Berkeley UNIX 4.3BSD[11]), computer algebra systems and programming languages (for example C99[12]) provide a special natural logarithm plus 1 function, alternatively named LNP1,[13][14] or log1p[12] to give more accurate results for logarithms close to zero by passing arguments x, also close to zero, to a function log1p(x), which returns the value ln(1+x), instead of passing a value y close to 1 to a function returning ln(y).[12][13][14] The function log1p avoids in the floating point arithmetic a near cancelling of the absolute term 1 with the second term from the Taylor expansion of the ln, thereby allowing for a high accuracy for both the argument and the result near zero.[13][14]

In addition to base e the IEEE 754-2008 standard defines similar logarithmic functions near 1 for binary and decimal logarithms: ${\displaystyle \log _{2}(1+x)}$ and ${\displaystyle \log _{10}(1+x)}$.

Similar inverse functions named "expm1",[12] "expm"[13][14] or "exp1m" exist as well, all with the meaning of expm1(x) = exp(x) - 1.[nb 2]

An identity in terms of the inverse hyperbolic tangent,

${\displaystyle \mathrm {log1p} (x)=\log(1+x)=2~\mathrm {artanh} \left({\frac {x}{2+x}}\right)\,,}$

gives a high precision value for small values of x on systems that do not implement log1p(x).

Computational complexity

Main article: Computational complexity of mathematical operations

The computational complexity of computing the natural logarithm (using the arithmetic-geometric mean) is O(M(n) ln n). Here n is the number of digits of precision at which the natural logarithm is to be evaluated and M(n) is the computational complexity of multiplying two n-digit numbers.

Continued fractions

While no simple continued fractions are available, several generalized continued fractions are, including:

${\displaystyle {\begin{aligned}\ln(1+x)&={\frac {x^{1}}{1}}-{\frac {x^{2}}{2}}+{\frac {x^{3}}{3}}-{\frac {x^{4}}{4}}+{\frac {x^{5}}{5}}-\cdots \\[5pt]&={\cfrac {x}{1-0x+{\cfrac {1^{2}x}{2-1x+{\cfrac {2^{2}x}{3-2x+{\cfrac {3^{2}x}{4-3x+{\cfrac {4^{2}x}{5-4x+\ddots }}}}}}}}}}\end{aligned}}}$

${\displaystyle {\begin{aligned}\ln \left(1+{\frac {x}{y}}\right)&={\cfrac {x}{y+{\cfrac {1x}{2+{\cfrac {1x}{3y+{\cfrac {2x}{2+{\cfrac {2x}{5y+{\cfrac {3x}{2+\ddots }}}}}}}}}}}}\\[5pt]&={\cfrac {2x}{2y+x-{\cfrac {(1x)^{2}}{3(2y+x)-{\cfrac {(2x)^{2}}{5(2y+x)-{\cfrac {(3x)^{2}}{7(2y+x)-\ddots }}}}}}}}\end{aligned}}}$

These continued fractions—particularly the last—converge rapidly for values close to 1. However, the natural logarithms of much larger numbers can easily be computed by repeatedly adding those of smaller numbers, with similarly rapid convergence.

For example, since 2 = 1.25³ × 1.024, the natural logarithm of 2 can be computed as:

${\displaystyle {\begin{aligned}\ln 2&=3\ln \left(1+{\frac {1}{4}}\right)+\ln \left(1+{\frac {3}{125}}\right)\\[8pt]&={\cfrac {6}{9-{\cfrac {1^{2}}{27-{\cfrac {2^{2}}{45-{\cfrac {3^{2}}{63-\ddots }}}}}}}}+{\cfrac {6}{253-{\cfrac {3^{2}}{759-{\cfrac {6^{2}}{1265-{\cfrac {9^{2}}{1771-\ddots }}}}}}}}.\end{aligned}}}$

Furthermore, since 10 = 1.25¹⁰ × 1.024³, even the natural logarithm of 10 similarly can be computed as:

${\displaystyle {\begin{aligned}\ln 10&=10\ln \left(1+{\frac {1}{4}}\right)+3\ln \left(1+{\frac {3}{125}}\right)\\[10pt]&={\cfrac {20}{9-{\cfrac {1^{2}}{27-{\cfrac {2^{2}}{45-{\cfrac {3^{2}}{63-\ddots }}}}}}}}+{\cfrac {18}{253-{\cfrac {3^{2}}{759-{\cfrac {6^{2}}{1265-{\cfrac {9^{2}}{1771-\ddots }}}}}}}}.\end{aligned}}}$

Complex logarithms

Main article: Complex logarithm

The exponential function can be extended to a function which gives a complex number as e^x for any arbitrary complex number x; simply use the infinite series with x complex. This exponential function can be inverted to form a complex logarithm that exhibits most of the properties of the ordinary logarithm. There are two difficulties involved: no x has e^x = 0; and it turns out that e^2πi = 1 = e⁰. Since the multiplicative property still works for the complex exponential function, e^z = e^z+2πki, for all complex z and integers k.

So the logarithm cannot be defined for the whole complex plane, and even then it is multi-valued – any complex logarithm can be changed into an "equivalent" logarithm by adding any integer multiple of 2πi at will. The complex logarithm can only be single-valued on the cut plane. For example, ln(i) = πi/2 or 5πi/2 or -3πi/2, etc.; and although i⁴ = 1, 4 log(i) can be defined as 2πi, or 10πi or −6πi, and so on.

• Plots of the natural logarithm function on the complex plane (principal branch)
• 
  z = Re(ln(x + yi))
• 
  z = abs(Im(ln(x + yi)))
• 
  z = abs(ln(x + yi))
• 
  Superposition of the previous three graphs

See also

• Approximating natural exponents (log base e)
• Napierian logarithm
• Logarithm of a matrix
• Logarithmic integral function
• Nicholas Mercator – first to use the term natural logarithm
• Polylogarithm
• Von Mangoldt function
• The number e

Notes

1. Including C, C++, SAS, MATLAB, Mathematica, Fortran, and some BASIC dialects
2. For a similar approach to reduce round-off errors of calculations for certain input values see trigonometric functions like versine, vercosine, coversine, covercosine, haversine, havercosine, hacoversine, hacovercosine, exsecant and excosecant.

References

1. G.H. Hardy and E.M. Wright, An Introduction to the Theory of Numbers, 4th Ed., Oxford 1975, footnote to paragraph 1.7: "log x is, of course, the 'Naperian' logarithm of x, to base e. 'Common' logarithms have no mathematical interest".
2. Mortimer, Robert G. (2005). Mathematics for physical chemistry (3rd ed.). Academic Press. p. 9. ISBN 0-12-508347-5. Extract of page 9
3. Burn, R. P. (2001). Alphonse Antonio de Sarasa and Logarithms. Historia Mathematica. pp. 28:1–17.
4. O'Connor, J. J.; Robertson, E. F. (September 2001). "The number e". The MacTutor History of Mathematics archive. Retrieved 2009-02-02.
5. Cajori, Florian (1991). A History of Mathematics (5th ed.). AMS Bookstore. p. 152. ISBN 0-8218-2102-4.
6. "Logarithmic Expansions" at Math2.org
7. Leonhard Euler, Introductio in Analysin Infinitorum. Tomus Primus. Bousquet, Lausanne 1748. Exemplum 1, p. 228; quoque in: Opera Omnia, Series Prima, Opera Mathematica, Volumen Octavum, Teubner 1922
8. OEIS: A002392
9. Sasaki, T.; Kanada, Y. (1982). "Practically fast multiple-precision evaluation of log(x)". Journal of Information Processing. 5 (4): 247–250. Retrieved 2011-03-30.
10. Ahrendt, Timm (1999). "Fast Computations of the Exponential Function". Stacs 99. Lecture Notes in Computer Science. 1564: 302–312. doi:10.1007/3-540-49116-3_28. ISBN 978-3-540-65691-3.
11. Beebe, Nelson H. F. (2017-08-22). "Chapter 10.4. Logarithm near one". The Mathematical-Function Computation Handbook - Programming Using the MathCW Portable Software Library (1 ed.). Salt Lake City, UT, USA: Springer International Publishing AG. pp. 290–292. doi:10.1007/978-3-319-64110-2. ISBN 978-3-319-64109-6. LCCN 2017947446. In 1987, Berkeley UNIX 4.3BSD introduced the log1p() function
12. Beebe, Nelson H. F. (2002-07-09). "Computation of expm1 = exp(x)−1" (PDF). 1.00. Salt Lake City, Utah, USA: Department of Mathematics, Center for Scientific Computing, University of Utah. Retrieved 2015-11-02.
13. HP 48G Series – Advanced User's Reference Manual (AUR) (4 ed.). Hewlett-Packard. December 1994 [1993]. HP 00048-90136, 0-88698-01574-2. Retrieved 2015-09-06.
14. HP 50g / 49g+ / 48gII graphing calculator advanced user's reference manual (AUR) (2 ed.). Hewlett-Packard. 2009-07-14 [2005]. HP F2228-90010. Retrieved 2015-10-10. Searchable PDF