Question 6 of 7: Root-Finding — Newton, Fixed-Point, and Bisection
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Examinations, December 2019 — 04-BS-5 Advanced Mathematics, 3 hours, closed book (approved Casio/Sharp calculator and one double-sided aid sheet permitted). Any five of the seven questions constitute a complete paper and all questions are of equal value; all seven are answered below as a full study resource.
Reference texts: Kreyszig, Advanced Engineering Mathematics, 10th ed. (Wiley) — Ch. 11 (Sturm–Liouville Problems, Fourier Series, Fourier Integrals and Transforms), Ch. 19 (Numerics in General: interpolation, numerical differentiation, Romberg integration, iterative equation solving), Ch. 20 (Numeric Linear Algebra: LU factorization). Supporting: Chapra & Canale, Numerical Methods for Engineers, 7th ed. — Ch. 5–6 (bracketing and open root-finding methods), Ch. 18 (interpolation), Ch. 22 (Romberg integration).
Find. (A) Three Newton iterates and, from the resulting root, the remaining two roots of the cubic. (B) A convergent fixed-point form and six iterates. (C) Four bisection iterates.
Approach. (A) Standard Newton update $x_{n+1}=x_n-f(x_n)/f'(x_n)$, then deflate the cubic by the found root and solve the resulting quadratic. (B) Test candidate rearrangements' derivatives at the root; pick the one with $|g'(x^\ast)|\lt1$. (C) Standard interval-halving with sign tracking.
Deflate the cubic. Dividing $x^3-26x^2+173x-325$ synthetically by $(x-3.3389600)$ gives the depressed quadratic $x^2+bx+c$ with $b=-26+3.3389600=-22.6610400$ and $c=173+3.3389600\,b=97.3389588$ (remainder $\approx-1.3\times10^{-5}$, confirming $x_3$ is an accurate root). Solving $x^2-22.6610400x+97.3389588=0$ by the quadratic formula,
$$\boxed{x=16.9023231\ \text{and}\ x=5.7587170}$$
(B) Choose a convergent fixed-point form. Rearranging $\ln(x+2)-x^2+6x-5=0$ as $x=g(x)=\dfrac{x^2-\ln(x+2)+5}{6}$ gives
$$g'(x)=\frac{2x-\frac{1}{x+2}}{6}$$
Near the root ($x^\ast\approx0.7605$), $g'(x^\ast)\approx0.193121$, and $|g'(x^\ast)|\lt1$, so this form converges (contraction mapping). By contrast, the alternative $x=\sqrt{\ln(x+2)+6x-5}$ has $g'(x^\ast)\approx4.183000\gt1$ at the same root and would diverge — it is rejected.
Iterate six times from $x_0=1.0$.
$$x_0=1.0000000\to x_1=0.8168980\to x_2=0.7719477\to x_3=0.7627255$$
$$\to x_4=0.7609221\to x_5=0.7605730\to \boxed{x_6=0.7605055}$$
The iterates are converging monotonically toward the true root $\approx0.7604893$, consistent with the small $|g'(x^\ast)|\approx0.19$ (fast linear convergence).
(C) Bisection, four iterations. $f(\alpha)=f(0.80)=0.0421220\gt0$, $f(\beta)=f(0.84)=-0.0138221\lt0$ (root bracketed).
Iter.
$a$
$b$
$c=(a+b)/2$
$f(c)$
New bracket
1
$0.8000000$
$0.8400000$
$0.8200000$
$+0.0142416$
$[0.82,0.84]$
2
$0.8200000$
$0.8400000$
$0.8300000$
$+0.0002326$
$[0.83,0.84]$
3
$0.8300000$
$0.8400000$
$0.8350000$
$-0.0067890$
$[0.83,0.835]$
4
$0.8300000$
$0.8350000$
$0.8325000$
$-0.0032768$
$[0.83,0.8325]$
After four bisections, $\boxed{x\approx0.8325000}$, with the bracket narrowed to $[0.8300000,0.8325000]$ (width $0.0025$, $1/16$ of the original $0.04$).