Past exam of the year 1 of the computer science course of the University of Oxford / 2019 / Discrete mathematics, continuous mathematics, and probability / 5 / c

For the vector equation, the Jacobian matrix isWrite the updated coordinates as $(X,Y)$. The second equation of the multivariable Newton iteration gives $X-Y=0$. Substituting $X=Y$ into its first linearized equation givesHence $h(x,y)=6x^5-20y^3$, and, whenever $h\ne 0$,

Starting from $(0,1)$ gives $(x_1,y_1)=(21/20,21/20)$. Every subsequent iterate remains on the diagonal, where this is exactly the scalar Newton root-finding iteration from (b). The value $f(21/20)>0$, so $21/20<r$.

For $1\le x<r$, concavity again puts the tangent's zero at or above $r$. Alsoso $r\le N(21/20)\le N(1)=8/7<6/5$. After one scalar step the sequence therefore enters the invariant interval proved in (b). It converges to $(r,r)$ at least quadratically. If the vector error is measured in the Euclidean norm, then $\Vert (x_k,y_k)-(r,r)\Vert =\sqrt{2}|x_k-r|$, so the scalar quadratic estimate also gives a quadratic vector estimate. The initial off-diagonal step causes no difficulty because its Jacobian matrix is invertible and lands exactly on the diagonal.