23-Chem-B3 Simulation, Modelling, and Optimization · December 2014
Question 6 of 7: Convergence of Newton's Method and Bisection
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams — December 2014 — 04-Chem-B3, Simulation, Modelling & Optimization. Closed-book; any five of the seven questions constitute a complete paper (equal value). All seven are solved below as a study resource.
Reference texts: T. F. Edgar, D. M. Himmelblau & L. S. Lasdon, Optimization of Chemical Processes (2nd ed., McGraw-Hill) — unconstrained & constrained optimization and linear programming; S. S. Rao, Engineering Optimization: Theory and Practice (4th ed., Wiley); S. C. Chapra & R. P. Canale, Numerical Methods for Engineers (7th ed., McGraw-Hill) — root finding (Newton, bisection) and difference equations; D. C. Montgomery & G. C. Runger, Applied Statistics and Probability for Engineers (7th ed., Wiley) — linear regression and hypothesis testing; E. Kreyszig, Advanced Engineering Mathematics (10th ed., Wiley) — the error function and difference equations.
Question 6: Convergence of Newton's Method and Bisection (20 marks)
Given. A function with a single root $\alpha$; along the axis the marked points are ordered $B,\ A,\ C,\ \alpha,\ D$. The curve rises to a maximum at $A$ (so $f'(A)=0$) and then descends through $\alpha$; thus $f\gt0$ at $B$, $A$ and $C$, while $f\lt0$ at $D$.
Find. Whether each iteration reaches the root $\alpha$.
Figure 3 — Reproduced sketch for Q6: single root α; f rises to a maximum at A and falls through α. f(B),f(A),f(C) > 0 and f(D) < 0.
Approach. For Newton, apply $x_{{n+1}} = x_n - f(x_n)/f'(x_n)$ and check the direction of the first step; for bisection, check whether the starting interval brackets the root (opposite signs at its ends).
(a) Newton from $x=A$. $A$ is the maximum, so $f'(A)=0$: the tangent is horizontal and the update $x_{{n+1}} = A - f(A)/0$ is undefined (it shoots to $\pm\infty$). Newton fails / diverges — the root is not reached.
(b) Newton from $x=B$. $B$ lies on the rising branch to the left of the maximum, where $f(B)\gt0$ and $f'(B)\gt0$. The step is $-f(B)/f'(B) \lt 0$, moving the iterate left, away from $\alpha$ (which lies to the right). The iteration diverges — the root is not reached.
(c) Bisection on $\overline{{AD}}$. $f(A)\gt0$ and $f(D)\lt0$ have opposite signs, so the interval brackets $\alpha$ (the only root inside). Bisection converges to $\alpha$.
(d) Bisection on $\overline{{BD}}$. $f(B)\gt0$ and $f(D)\lt0$ again have opposite signs, and $\alpha$ is the single sign change inside $\overline{{BD}}$. Bisection converges to $\alpha$.