NivaarExam PrepOfficial exam papers ↗

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)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

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$.

xBACαDf(x)Single root α; peak (f'=0) at A; f>0 at B,A,C and f<0 at D
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).

  1. (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.
  2. (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.
  3. (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$.
  4. (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$.
Method & startReaches $\alpha$?Reason
(a) Newton, $x=A$No$f'(A)=0$ (maximum) — step undefined
(b) Newton, $x=B$NoSteps away from $\alpha$ (wrong side of the peak)
(c) Bisection, $\overline{{AD}}$Yes$f(A)\gt0,\ f(D)\lt0$ — brackets the root
(d) Bisection, $\overline{{BD}}$Yes$f(B)\gt0,\ f(D)\lt0$ — brackets the root