Past exam of the mathematics course of the University of Cambridge / 2001 / iii / Paper 53 / 3

At a current approximation $p$, Taylor's theorem replaces $\Phi (p+h)=0$ by the affine equation $\Phi (p)+D\Phi (p)h=0$. Solving this equation gives the Newton iteration in a Banach spaceThe inverse is a bounded linear inverse by the bounded inverse theorem. The displayed operator is defined where the derivative is invertible, in particular on the specified open domain; the hypotheses do not supply a globally defined Newton operator on all of $\mathcal{X}$. If $\Phi (\bar{p})=0$ with $\bar{p}$ in that domain, then $N(\bar{p})=\bar{p}$.

For quadratic convergence, choose a sufficiently small neighbourhood of $\bar{p}$. Continuity of the inverse derivative and of $D^2\Phi$ gives constants $M,L$ there such thatPut $e=p-\bar{p}$ and $B(p)=[D\Phi (p){]}^{-1}$. The quadratic Newton error bound in a Banach space follows from the integral remainder:soTake a closed radius-$r$ ball lying in the neighbourhood with $Cr<1$. This inequality maps the ball into itself and makes errors tend to zero. More explicitly, for $e_n=\Vert p_n-\bar{p}\Vert$, induction gives $C{e}_{n+j}\le (C{e}_n{)}^{2^j}$. Hence every sufficiently close initial point converges to $\bar{p}$ at least quadratically; special maps can converge still faster.

For a fixed-point equation, apply the same method to $F(p)=\Phi (p)-p$. Wherever $D\Phi (p)-I$ is invertible, the exact adapted map isEvery fixed point of $\Phi$ is fixed by $A$, and conversely on this domain. This additional invertibility condition is distinct from invertibility of $D\Phi$; for example, $D\Phi =I$ is invertible but $D\Phi -I=0$ is not.

For certification, a frozen Newton correction for a fixed-point equation is often more useful. Choose a bounded, injective approximate inverse $J$ of $D\Phi (p_0)-I$ and setAgain $A_J(p)=p$ if and only if $\Phi (p)=p$. A well-chosen $J$ makes the corrected operator contract even if the original operator has an expanding direction. The subsequent contraction mapping argument requires the stated derivative and residual bounds; it does not assume that arbitrary inverse approximations provide them.

Let $B=\{p:\Vert p-p_0\Vert \le r\}$ and write $\epsilon =\Vert A(p_0)-p_0\Vert$. Integrating the Fréchet derivative along the line segment between two points in the convex set givesand the operator norm bound impliesIn particular, for $p\in B$,Thus the a posteriori contraction ball is mapped into itself. The ball is a closed subset of a Banach space, hence is a complete metric space. Applying the contraction mapping theorem to $A:B\to B$ gives existence and uniqueness:The estimate also gives $\Vert p_{\ast }-p_0\Vert \le \epsilon /(1-\kappa )$, by applying the same triangle inequality to $A(p_{\ast })=p_{\ast }$. Equality in the residual bound is allowed; strict contraction mapping comes from $\kappa <1$, and a fixed point on the closed ball's boundary is included.

The Lanford contraction proof of the Feigenbaum fixed point replaces an unstable renormalization operator by a contracting Newton correction. Use a coefficient Banach space for normalized even maps, represented by $p(x)=1-x^{2}h(x^2)$ with an absolutely summable coefficient norm on $h$. Its complex domain and a small ball are chosen so that every map is admissible and the renormalization composition is analytic there.

A high-accuracy polynomial $p_0$ gives an approximate fixed point. Use a bounded invertible approximate inverse $J$ for $D\mathcal{T}(p_0)-I$, and defineRigorous estimates certifyThe finite coefficient computations use interval arithmetic with controlled rounding; analytic estimates bound the infinite tail. Thus the certification applies to the full function space, rather than merely a truncated polynomial system.

Part (b) gives a unique fixed point of $A$ in the certified ball. Invertibility of $J$ makes it a unique fixed point $g$ of $\mathcal{T}$ there. The ball also preserves the quadratic critical maximum and nontrivial shape, excluding the constant formal normalized solution. Consequently an actual locally unique analytic Feigenbaum fixed point exists, with the error bound $\Vert g-p_0\Vert \le \epsilon /(1-\kappa )$. The numerical approximation locates it; rigorous residual, derivative and tail bounds establish its existence and uniqueness. This conclusion is local, not a claim of global uniqueness among all analytic solutions.