Introductory university-level calculus, linear algebra, abstract algebra, probability, statistics, and stochastic processes.
Newton's method
Get link
Facebook
X
Pinterest
Email
Other Apps
-
Newton's method (the Newton-Raphson method) is a very powerful numerical method for solving nonlinear equations.
Suppose we'd like to solve a non-linear equation \(f(x) = 0\), where \(f(x)\) is a (twice) differentiable non-linear function. Newton's method generates a sequence of numbers \(c_1, c_2, c_3, \cdots\) that converges to a solution of the equation. That is, if \(\alpha\) is a solution (i.e., \(f(\alpha) = 0\)) then,
\[\lim_{n\to\infty}c_n = \alpha,\]
and this sequence \(\{c_n\}\) is generated by a series of linear approximations of the function \(f(x)\).
Theorem (Newton's method)
Let \(f(x)\) be a function that is twice differentiable on an open interval \(I\) that contains the closed interval \([a, b]\) (i.e., \([a,b]\subset I\)) and satisfy the following conditions:
\(f(a) < 0\) and \(f(b) > 0\);
For all \(x\in [a, b]\), \(f'(x) > 0\) and \(f''(x) > 0\).
where \(\alpha\) is the unique solution of \(f(\alpha) = 0\).
Proof. Since \(f(x)\) is (twice) differentiable, it is continuous on \([a,b]\). Since \(f(a) < 0\) and \(f(b) > 0\), clearly \(f(a) \neq f(b)\). Therefore, by the Intermediate Value Theorem, there exists an \(\alpha \in [a,b]\) such that \(f(\alpha) = 0\). Moreover, since \(f'(x) > 0\), \(f(x)\) is strictly monotone increasing on \([a,b]\), this \(\alpha\) is unique (i.e., there is only one such \(\alpha\) that satisfies \(f(\alpha) = 0\)).
Let us pick any \(c\in [a, b]\) and find the tangent line at \(x = c\), which is given by
\[y = f'(c)(x - c) + f(c).\]
The \(x\)-intercept \(x=c'\) of this line is obtained by
\[0 = f'(c)(c' - c) + f(c),\]
and hence
\[c' = c - \frac{f(c)}{f'(c)}.\]
Therefore if \(c = c_n\), we have \(c' = c_{n+1}\).
Moreover, since \(y=f(x)\) is convex (because \(f''(x) > 0\)), the tangent line is always resides below \(y=f(x)\) on \([a,b]\) so we have \(c' > \alpha\). Since \(f'(x) > 0\) and \(f(x) > 0\) (when \(x > \alpha\)), we have
\[c_{n+1} = c_n -\frac{f(c_n)}{f'(c_n)} < c_n.\]
Thus, \(\{c_n\}\) is a strictly monotone decreasing sequence bounded below. By the Monotone Convergence Theorem, \(\{c_n\}\) converges to a certain value. [See also: Monotone sequences and Cauchy sequences for the Monotone Convergence Theorem.]
Let \(\beta = \lim_{n\to\infty}c_n\), then \(\alpha \leq \beta < b\). If we take the limit \(n \to \infty\) on both sides of Eq. (eq:newton), we have
\[\beta = \beta - \frac{f(\beta)}{f'(\beta)},\]
from which we conclude
\[f(\beta) = 0.\]
However, \(\alpha\) is the only solution to \(f(x) = 0\). Hence \(\beta = \alpha\). ■
Example. Let \(d\) be any positive real number. Using Newton's method, let us find an approximation of \(\sqrt{d}\). First, let us define \(f(x) = x^2 - d\). Then the solution to \(f(x) = 0\) is \(x=\sqrt{d}\) (provided \(x \geq 0\)). \(f(0) = -d < 0\), and \(f'(x) = 2x > 0\) for \(x > 0\), and \(f''(0) = 2 > 0\). By choosing sufficiently small \(a\) and sufficiently large \(b\), \(f(x)\) can satisfy the conditions required for Newton's method. If we define the sequence \(\{c_n\}\) by
Defining the birth process Consider a colony of bacteria that never dies. We study the following process known as the birth process , also known as the Yule process . The colony starts with \(n_0\) cells at time \(t = 0\). Assume that the probability that any individual cell divides in the time interval \((t, t + \delta t)\) is proportional to \(\delta t\) for small \(\delta t\). Further assume that each cell division is independent of others. Let \(\lambda\) be the birth rate. The probability of a cell division for a population of \(n\) cells during \(\delta t\) is \(\lambda n \delta t\). We assume that the probability that two or more births take place in the time interval \(\delta t\) is \(o(\delta t)\). That is, it can be ignored. Consequently, the probability that no cell divides during \(\delta t\) is \(1 - \lambda n \delta t - o(\delta t)\). Note that this process is an example of the Markov chain with states \({n_0}, {n_0 + 1}, {n_0 + 2}...
Generational growth Consider the following scenario (see the figure below): A single individual (cell, organism, etc.) produces \(j (= 0, 1, 2, \cdots)\) descendants with probability \(p_j\), independently of other individuals. The probability of this reproduction, \(\{p_j\}\), is known. That individual produces no further descendants after the first (if any) reproduction. These descendants each produce further descendants at the next subsequent time with the same probabilities. This process carries on, creating successive generations. Figure 1. An example of the branching process. Let \(X_n\) be the random variable representing the population size (number of individuals) of generation \(n\). In the above figure, we have \(X_0 = 1\), \(X_1=4\), \(X_2 = 7\), \(X_3=12\), \(X_4 = 9.\) We shall assume \(X_0 = 1\) as the initial condition. Ideally, our goal would be to find how the population size grows through generations, that is, to find the probability \(\Pr(X_n = k)\) for e...
In mathematics, we must prove (almost) everything and the proofs must be done logically and rigorously. Therefore, we need some understanding of basic logic. Here, I will informally explain some rudimentary formal logic. Definitions (Proposition): A proposition is a statement that is either true or false. "True" and "false" are called the truth values, and are often denoted \(\top\) and \(\bot\). Here is an example. "Dr. Akira teaches at UBD." is a statement that is either true or false (we understand the existence of Dr. Akira and UBD), hence a proposition. The following statement is also a proposition, although we don't know if it's true or false (yet): Any even number greater than or equal to 4 is equal to a sum of two primes. See also: Goldbach's conjecture Next, we define several operations on propositions. Note that propositions combined with these operations are again propositions. (Conjunction, logical "and"): Let \(P\)...
Comments
Post a Comment