using CalculusWithJulia
using Plots
plotly()
using SymPy
32 L’Hospital’s Rule
This section uses these add-on packages:
Let’s return to limits of the form \(\lim_{x \rightarrow c}f(x)/g(x)\) which have an indeterminate form of \(0/0\) if both are evaluated at \(c\). The typical example being the limit considered by Euler:
\[ \lim_{x\rightarrow 0} \frac{\sin(x)}{x}. \]
We know this is \(1\) using a bound from geometry, but might also guess this is one, as we know from linearization near \(0\) that we have \(\sin(x) \approx x\) or, more specifically:
\[ \sin(x) = x - \sin(\xi)x^2/2, \quad 0 < \xi < x. \]
This would yield:
\[ \lim_{x \rightarrow 0} \frac{\sin(x)}{x} = \lim_{x\rightarrow 0} \frac{x -\sin(\xi) x^2/2}{x} = \lim_{x\rightarrow 0} 1 - \sin(\xi) \cdot x/2 = 1. \]
This is because we know \(\sin(\xi) x/2\) has a limit of \(0\), when \(|\xi| \leq |x|\).
That doesn’t look any easier, as we worried about the error term, but if just mentally replaced \(\sin(x)\) with \(x\) - which it basically is near \(0\) - then we can see that the limit should be the same as \(x/x\) which we know is \(1\) without thinking.
Basically, we found that in terms of limits, if both \(f(x)\) and \(g(x)\) are \(0\) at \(c\), that we might be able to just take this limit: \((f(c) + f'(c) \cdot(x-c)) / (g(c) + g'(c) \cdot (x-c))\) which is just \(f'(c)/g'(c)\).
Wouldn’t that be nice? We could find difficult limits just by differentiating the top and the bottom at \(c\) (and not use the messy quotient rule).
Well, in fact that is more or less true, a fact that dates back to L’Hospital - who wrote the first textbook on differential calculus - though this result is likely due to one of the Bernoulli brothers.
L’Hospital’s rule: Suppose:
- that \(\lim_{x\rightarrow c+} f(c) =0\) and \(\lim_{x\rightarrow c+} g(c) =0\),
- that \(f\) and \(g\) are differentiable in \((c,b)\), and
- that \(g(x)\) exists and is non-zero for all \(x\) in \((c,b)\),
then if the following limit exists: \(\lim_{x\rightarrow c+}f'(x)/g'(x)=L\) it follows that \(\lim_{x \rightarrow c+}f(x)/g(x) = L\).
That is if the right limit of \(f(x)/g(x)\) is indeterminate of the form \(0/0\), but the right limit of \(f'(x)/g'(x)\) is known, possibly by simple continuity, then the right limit of \(f(x)/g(x)\) exists and is equal to that of \(f'(x)/g'(x)\).
The rule equally applies to left limits and limits at \(c\). Later it will see there are other generalizations.
To apply this rule to Euler’s example, \(\sin(x)/x\), we just need to consider that:
\[ L = 1 = \lim_{x \rightarrow 0}\frac{\cos(x)}{1}, \]
So, as well, \(\lim_{x \rightarrow 0} \sin(x)/x = 1\).
This is due to \(\cos(x)\) being continuous at \(0\), so this limit is just \(\cos(0)/1\). (More importantly, the tangent line expansion of \(\sin(x)\) at \(0\) is \(\sin(0) + \cos(0)x\), so that \(\cos(0)\) is why this answer is as it is, but we don’t need to think in terms of \(\cos(0)\), but rather the tangent-line expansion, which is \(\sin(x) \approx x\), as \(\cos(0)\) appears as the coefficient.
In Gruntz, in a reference attributed to Speiss, we learn that L’Hospital was a French Marquis who was taught in \(1692\) the calculus of Leibniz by Johann Bernoulli. They made a contract obliging Bernoulli to leave his mathematical inventions to L’Hospital in exchange for a regular compensation. This result was discovered in \(1694\) and appeared in L’Hospital’s book of \(1696\).
Examples
- Consider this limit at \(0\): \((a^x - 1)/x\). We have \(f(x) =a^x-1\) has \(f(0) = 0\), so this limit is indeterminate of the form \(0/0\). The derivative of \(f(x)\) is \(f'(x) = a^x \log(a)\) which has \(f'(0) = \log(a)\). The derivative of the bottom is also \(1\) at \(0\), so we have:
\[ \log(a) = \frac{\log(a)}{1} = \frac{f'(0)}{g'(0)} = \lim_{x \rightarrow 0}\frac{f'(x)}{g'(x)} = \lim_{x \rightarrow 0}\frac{f(x)}{g(x)} = \lim_{x \rightarrow 0}\frac{a^x - 1}{x}. \]
Why rewrite in the “opposite” direction? Because the theorem’s result – \(L\) is the limit – is only true if the related limit involving the derivative exists. We don’t do this in the following, but did so here to emphasize the need for the limit of the ratio of the derivatives to exist.
- Consider this limit:
\[ \lim_{x \rightarrow 0} \frac{e^x - e^{-x}}{x}. \]
It too is of the indeterminate form \(0/0\). The derivative of the top is \(e^x + e^{-x}\), which is \(2\) when \(x=0\), so the ratio of \(f'(0)/g'(0)\) is seen to be \(2\) By continuity, the limit of the ratio of the derivatives is \(2\). Then by L’Hospital’s rule, the limit above is \(2\).
- Sometimes, L’Hospital’s rule must be applied twice. Consider this limit:
\[ \lim_{x \rightarrow 0} \frac{1-\cos(x)}{x^2} \]
By L’Hospital’s rule if this following limit exists, the two will be equal:
\[ \lim_{x \rightarrow 0} \frac{\sin(x)}{2x}. \]
But if we didn’t guess the answer, we see that this new problem is also indeterminate of the form \(0/0\). So, repeating the process, this new limit will exist and be equal to the following limit, should it exist:
\[ \lim_{x \rightarrow 0} \frac{\cos(x)}{2} = 1/2. \]
As \(L = 1/2\) for this related limit, it must also be the limit of the original problem, by L’Hospital’s rule.
- Our “intuitive” limits can bump into issues. Take for example the limit of \((\sin(x)-x)/x^2\) as \(x\) goes to \(0\). Using \(\sin(x) \approx x\) makes this look like \(0/x^2\) which is still indeterminate. (Because the difference is higher order than \(x\).) Using L’Hospitals, says this limit will exist (and be equal) if the following one does:
\[ \lim_{x \rightarrow 0} \frac{\cos(x) - 1}{2x}. \]
This particular limit is indeterminate of the form \(0/0\), so we again try L’Hospital’s rule and consider
\[ \lim_{x \rightarrow 0} \frac{-\sin(x)}{2} = 0 \]
So as this limit exists, working backwards, the original limit in question will also be \(0\).
- This example comes from the Wikipedia page. It “proves” a discrete approximation for the second derivative.
Show if \(f''(x)\) exists at \(c\) and is continuous at \(c\), then
\[ f''(c) = \lim_{h \rightarrow 0} \frac{f(c + h) - 2f(c) + f(c-h)}{h^2}. \]
This will follow from two applications of L’Hospital’s rule to the right-hand side. The first says, the limit on the right is equal to this limit, should it exist:
\[ \lim_{h \rightarrow 0} \frac{f'(c+h) - 0 - f'(c-h)}{2h}. \]
We have to be careful, as we differentiate in the \(h\) variable, not the \(c\) one, so the chain rule brings out the minus sign. But again, as we still have an indeterminate form \(0/0\), this limit will equal the following limit should it exist:
\[ \lim_{h \rightarrow 0} \frac{f''(c+h) - 0 - (-f''(c-h))}{2} = \lim_{c \rightarrow 0}\frac{f''(c+h) + f''(c-h)}{2} = f''(c). \]
That last equality follows, as it is assumed that \(f''(x)\) exists at \(c\) and is continuous, that is, \(f''(c \pm h) \rightarrow f''(c)\).
The expression above finds use when second derivatives are numerically approximated. (The middle expression is the basis of the central-finite difference approximation to the derivative.)
- L’Hospital himself was interested in this limit for \(a > 0\) (math overflow)
\[ \lim_{x \rightarrow a} \frac{\sqrt{2a^3\cdot x-x^4} - a\cdot(a^2\cdot x)^{1/3}}{ a - (a\cdot x^3)^{1/4}}. \]
These derivatives can be done by hand, but to avoid any minor mistakes we utilize SymPy
taking care to use rational numbers for the fractional powers, so as not to lose precision through floating point roundoff:
@syms a::positive x::positive
f(x) = sqrt(2a^3*x - x^4) - a * (a^2*x)^(1//3)
g(x) = a - (a*x^3)^(1//4)
g (generic function with 1 method)
We can see that at \(x=a\) we have the indeterminate form \(0/0\):
f(a), g(a)
(0, 0)
What about the derivatives?
= diff(f(x),x), diff(g(x),x)
fp, gp fp(x=>a), gp(x=>a)
(-4*a/3, -3/4)
Their ratio will not be indeterminate, so the limit in question is just the ratio:
fp(x=>a) / gp(x=>a)
\(\frac{16 a}{9}\)
Of course, we could have just relied on limit
, which knows about L’Hospital’s rule:
limit(f(x)/g(x), x => a)
\(\frac{16 a}{9}\)
32.1 Idea behind L’Hospital’s rule
A first proof of L’Hospital’s rule takes advantage of Cauchy’s generalization of the mean value theorem to two functions. Suppose \(f(x)\) and \(g(x)\) are continuous on \([c,b]\) and differentiable on \((c,b)\). On \((c,x)\), \(c < x < b\) there exists a \(\xi\) with \(f'(\xi) \cdot (g(x) - g(c)) = g'(\xi) \cdot (f(x) - f(c))\). In our formulation, both \(f(c)\) and \(g(c)\) are zero, so we have, provided we know that \(g(x)\) is non zero, that \(f(x)/g(x) = f'(\xi)/g'(\xi)\) for some \(\xi\), \(c < \xi < x\). That the right-hand side has a limit as \(x \rightarrow c+\) is true by the assumption that the limit of the ratio of the derivatives exists. (The \(\xi\) part can be removed by considering it as a composition of a function going to \(c\).) Thus the right limit of the ratio \(f/g\) is known.
32.2 Generalizations
L’Hospital’s rule generalizes to other indeterminate forms, in particular the indeterminate form \(\infty/\infty\) can be proved at the same time as \(0/0\) with a more careful proof.
The value \(c\) in the limit can also be infinite. Consider this case with \(c=\infty\):
\[\begin{align*} \lim_{x \rightarrow \infty} \frac{f(x)}{g(x)} &= \lim_{x \rightarrow 0} \frac{f(1/x)}{g(1/x)} \end{align*}\]
L’Hospital’s limit applies as \(x \rightarrow 0\), so we differentiate to get:
\[\begin{align*} \lim_{x \rightarrow 0} \frac{[f(1/x)]'}{[g(1/x)]'} &= \lim_{x \rightarrow 0} \frac{f'(1/x)\cdot(-1/x^2)}{g'(1/x)\cdot(-1/x^2)}\\ &= \lim_{x \rightarrow 0} \frac{f'(1/x)}{g'(1/x)}\\ &= \lim_{x \rightarrow \infty} \frac{f'(x)}{g'(x)}, \end{align*}\]
assuming the latter limit exists, L’Hospital’s rule assures the equality
\[ \lim_{x \rightarrow \infty} \frac{f(x)}{g(x)} = \lim_{x \rightarrow \infty} \frac{f'(x)}{g'(x)}, \]
Examples
For example, consider
\[ \lim_{x \rightarrow \infty} \frac{x}{e^x}. \]
We see it is of the form \(\infty/\infty\). Taking advantage of the fact that L’Hospital’s rule applies to limits at \(\infty\), we have that this limit will exist and be equal to this one, should it exist:
\[ \lim_{x \rightarrow \infty} \frac{1}{e^x}. \]
This limit is, of course, \(0\), as it is of the form \(1/\infty\). It is not hard to build up from here to show that for any integer value of \(n>0\) that:
\[ \lim_{x \rightarrow \infty} \frac{x^n}{e^x} = 0. \]
This is an expression of the fact that exponential functions grow faster than polynomial functions.
Similarly, powers grow faster than logarithms, as this limit shows, which is indeterminate of the form \(\infty/\infty\):
\[ \lim_{x \rightarrow \infty} \frac{\log(x)}{x} = \lim_{x \rightarrow \infty} \frac{1/x}{1} = 0, \]
the first equality by L’Hospital’s rule, as the second limit exists.
32.3 Other indeterminate forms
Indeterminate forms of the type \(0 \cdot \infty\), \(0^0\), \(\infty^\infty\), \(\infty - \infty\) can be re-expressed to be in the form \(0/0\) or \(\infty/\infty\) and then L’Hospital’s theorem can be applied.
Example: rewriting \(0 \cdot \infty\)
What is the limit \(x \log(x)\) as \(x \rightarrow 0+\)? The form is \(0\cdot \infty\), rewriting, we see this is just:
\[ \lim_{x \rightarrow 0+}\frac{\log(x)}{1/x}. \]
L’Hospital’s rule clearly applies to one-sided limits, as well as two (our proof sketch used one-sided limits), so this limit will equal the following, should it exist:
\[ \lim_{x \rightarrow 0+}\frac{1/x}{-1/x^2} = \lim_{x \rightarrow 0+} -x = 0. \]
Example: rewriting \(0^0\)
What is the limit \(x^x\) as \(x \rightarrow 0+\)? The expression is of the form \(0^0\), which is indeterminate. (Even though floating point math defines the value as \(1\).) We can rewrite this by taking a log:
\[ x^x = \exp(\log(x^x)) = \exp(x \log(x)) = \exp(\log(x)/(1/x)). \]
Be just saw that \(\lim_{x \rightarrow 0+}\log(x)/(1/x) = 0\). So by the rules for limits of compositions and the fact that \(e^x\) is continuous, we see \(\lim_{x \rightarrow 0+} x^x = e^0 = 1\).
Example: rewriting \(\infty - \infty\)
A limit \(\lim_{x \rightarrow c} f(x) - g(x)\) of indeterminate form \(\infty - \infty\) can be reexpressed to be of the from \(0/0\) through the transformation:
\[\begin{align*} f(x) - g(x) &= f(x)g(x) \cdot (\frac{1}{g(x)} - \frac{1}{f(x)}) \\ &= \frac{\frac{1}{g(x)} - \frac{1}{f(x)}}{\frac{1}{f(x)g(x)}}. \end{align*}\]
Applying this to
\[ L = \lim_{x \rightarrow 1} \big(\frac{x}{x-1} - \frac{1}{\log(x)}\big) \]
We get that \(L\) is equal to the following limit:
\[ \lim_{x \rightarrow 1} \frac{\log(x) - \frac{x-1}{x}}{\frac{x-1}{x} \log(x)} = \lim_{x \rightarrow 1} \frac{x\log(x)-(x-1)}{(x-1)\log(x)} \]
In SymPy
we have (using italic variable names to avoid a problem with the earlier use of f
as a function):
= x * log(x) - (x-1)
𝑓 = (x-1) * log(x)
𝑔 𝑓(1), 𝑔(1)
(0, 0)
L’Hospital’s rule applies to the form \(0/0\), so we try:
= diff(𝑓, x)
𝑓′ = diff(𝑔, x)
𝑔′ 𝑓′(1), 𝑔′(1)
(0, 0)
Again, we get the indeterminate form \(0/0\), so we try again with second derivatives:
= diff(𝑓, x, x)
𝑓′′ = diff(𝑔, x, x)
𝑔′′ 𝑓′′(1), 𝑔′′(1)
(1, 2)
From this we see the limit is \(1/2\), as could have been done directly:
limit(𝑓/𝑔, x=>1)
\(\frac{1}{2}\)
32.4 The assumptions are necessary
Example: the limit existing is necessary
The following limit is easily seen by comparing terms of largest growth:
\[ 1 = \lim_{x \rightarrow \infty} \frac{x - \sin(x)}{x} \]
However, the limit of the ratio of the derivatives does not exist:
\[ \lim_{x \rightarrow \infty} \frac{1 - \cos(x)}{1}, \]
as the function just oscillates. This shows that L’Hospital’s rule does not apply when the limit of the the ratio of the derivatives does not exist.
Example: the assumptions matter
This example comes from the thesis of Gruntz to highlight possible issues when computer systems do simplifications.
Consider:
\[ \lim_{x \rightarrow \infty} \frac{1/2\sin(2x) +x}{\exp(\sin(x))\cdot(\cos(x)\sin(x)+x)}. \]
If we apply L’Hospital’s rule using simplification we have:
u(x) = 1//2*sin(2x) + x
v(x) = exp(sin(x))*(cos(x)*sin(x) + x)
= diff(u(x),x), diff(v(x),x)
up, vp limit(simplify(up/vp), x => oo)
\(0\)
However, this answer is incorrect. The reason being subtle. The simplification cancels a term of \(\cos(x)\) that appears in the numerator and denominator. Before cancellation, we have vp
will have infinitely many zero’s as \(x\) approaches \(\infty\) so L’Hospital’s won’t apply (the limit won’t exist, as every \(2\pi\) the ratio is undefined so the function is never eventually close to some \(L\)).
This ratio has no limit, as it oscillates, as confirmed by SymPy
:
limit(u(x)/v(x), x => oo)
\(\left\langle e^{-1}, e\right\rangle\)
32.5 Questions
Question
This function \(f(x) = \sin(5x)/x\) is indeterminate at \(x=0\). What type?
Question
This function \(f(x) = \sin(x)^{\sin(x)}\) is indeterminate at \(x=0\). What type?
Question
This function \(f(x) = (x-2)/(x^2 - 4)\) is indeterminate at \(x=2\). What type?
Question
This function \(f(x) = (g(x+h) - g(x-h)) / (2h)\) (\(g\) is continuous) is indeterminate at \(h=0\). What type?
Question
This function \(f(x) = x \log(x)\) is indeterminate at \(x=0\). What type?
Question
Does L’Hospital’s rule apply to this limit:
\[ \lim_{x \rightarrow \pi} \frac{\sin(\pi x)}{\pi x}. \]
Question
Use L’Hospital’s rule to find the limit
\[ L = \lim_{x \rightarrow 0} \frac{4x - \sin(x)}{x}. \]
What is \(L\)?
Question
Use L’Hospital’s rule to find the limit
\[ L = \lim_{x \rightarrow 0} \frac{\sqrt{1+x} - 1}{x}. \]
What is \(L\)?
Question
Use L’Hospital’s rule one or more times to find the limit
\[ L = \lim_{x \rightarrow 0} \frac{x - \sin(x)}{x^3}. \]
What is \(L\)?
Question
Use L’Hospital’s rule one or more times to find the limit
\[ L = \lim_{x \rightarrow 0} \frac{1 - x^2/2 - \cos(x)}{x^3}. \]
What is \(L\)?
Question
Use L’Hospital’s rule one or more times to find the limit
\[ L = \lim_{x \rightarrow \infty} \frac{\log(\log(x))}{\log(x)}. \]
What is \(L\)?
Question
By using a common denominator to rewrite this expression, use L’Hospital’s rule to find the limit
\[ L = \lim_{x \rightarrow 0} \frac{1}{x} - \frac{1}{\sin(x)}. \]
What is \(L\)?
Question
Use L’Hospital’s rule to find the limit
\[ L = \lim_{x \rightarrow \infty} \log(x)/x \]
What is \(L\)?
Question
Using L’Hospital’s rule, does
\[ \lim_{x \rightarrow 0+} x^{\log(x)} \]
exist?
Consider \(x^{\log(x)} = e^{\log(x)\log(x)}\).
Question
Using L’Hospital’s rule, find the limit of
\[ \lim_{x \rightarrow 1} (2-x)^{\tan(\pi/2 \cdot x)}. \]
(Hint, express as \(\exp^{\tan(\pi/2 \cdot x) \cdot \log(2-x)}\) and take the limit of the resulting exponent.)