Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 129 4 Solution Created 2026-10-03 Updated 2026-10-05
Uniform expectations require nonempty finite sets; assume this throughout the analytic argument. For real-valued functions on , define the box norm byAll four variables are sampled as independent random variables with replacement, so repeated coordinates are included. The defining fourth power equalsIf it vanishes, every summand vanishes. In particular, the terms give for each , so . Absolute homogeneity of a norm follows directly: .
For the triangle inequality, prove the mixed box Cauchy-Schwarz inequality. WriteSeparating the two averages and applying the Cauchy-Schwarz inequality in givesExpanding the square and reversing the order of finite sums,by a second Cauchy-Schwarz inequality. Thus . Expand into its sixteen multilinear terms and apply this bound to each term. Their upper bounds sum to . Taking fourth roots yields the triangle inequality. Together with definiteness of a norm and absolute homogeneity of a norm, this proves the box norm is a norm. Complex functions require conjugates and are outside the printed real-valued convention.
For the bilinear correlation bound for the box norm, use the normalized L2 norm: let and , . The Cauchy-Schwarz inequality in givesThe remaining factor is . By the Cauchy-Schwarz inequality in , its magnitude is at mostTaking square roots proves
For the tripartite graph, assume its parts are nonempty and write for its adjacency indicator functions. Let , so . The normalized triangle count isThe constant-degree hypothesis is exactly for each . Therefore the contribution of the constant isFor fixed , apply the bilinear correlation bound for the box norm with and . Since these are indicator functions, and , the fraction of adjacent to . Averaging and applying the Cauchy-Schwarz inequality in gives the sharper triangle counting with one box-uniform pair and constant opposite degree estimateMultiplying by yieldsThe middle bound also handles or , when the triangle count is zero.
To connect explicitly with the suggested expansion, put , . Both have overall mean zero, and for every . Consequently ; all terms with constant reduce to . The terms containing combine into the single error just estimated. This is where the exact degree assumption is used. If any part is empty the triangle count itself is zero, but the printed uniform expectations and edge density of a bipartite graph values would be undefined; the nonempty-parts convention must therefore be stated rather than dividing by zero.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 131 3 Solution Created 2026-10-03 Updated 2026-10-05
Write and . A Jacobi field is a vector field along the geodesic satisfyingThe curvature convention is that used in Question 1. A geodesic variation is a smooth map with whose -curves are affinely parametrized geodesics. Torsion-freeness gives . Differentiating and commuting covariant derivatives gives the Jacobi equation for .
For the converse, let and . Choose a curve with and . parallel transport along it identifies its tangent spaces with . Set , so and . The geodesics with initial data give a variation . Smooth dependence on initial conditions and compactness of the original time interval ensure that this family is defined for all after shrinking , even if is incomplete. Its variation field has the initial values , so uniqueness for the linear Jacobi equation identifies it with . This proves the realization of Jacobi fields by geodesic variations.
The endpoint-vanishing pointwise normal fields form a vector space, since their conditions and equation are linear. If is nonconstant, the map is injective because and zero initial derivative force the zero solution. Differentiating gives . Hence the endpoint-vanishing normal Jacobi fields have dimension at most n-1 bound isFor a constant geodesic, , and both zero endpoint values force , so the bound still holds.
On the unit round , let . Its speed is . The constant ambient vectors orthogonal to give independent parallel normal fields . Since on normal fields, the fields satisfy the Jacobi equation and vanish at both antipodal endpoints. They attain dimension .
For the last clause, interpret a closed geodesic as a nonconstant smoothly periodic geodesic. A constant loop has length zero and cannot be shortened. Parametrize on at constant nonzero speed. parallel transport once around it fixes . Because the manifold is orientable, this transport preserves orientation; on the normal space of dimension , which is odd, it lies in . An odd-dimensional special orthogonal transformation has a fixed vector: nonreal eigenvalues pair with their conjugates, real eigenvalues are , and determinant one in odd dimension forces an eigenvalue .
Transport such a nonzero normal fixed vector around the loop. It gives a periodic parallel normal field , with and . Its jets also agree at the seam, so it is smooth as a field on the parametrizing circle. For small , is a smooth variation through closed curves, providing their homotopy to .
For energy , the first variation vanishes at the closed geodesic. The permitted second variation of geodesic energy has no endpoint term for this periodic variation, soThus for small nonzero . The Cauchy-Schwarz inequality gives , while constant speed gives . Consequently the instability of a closed geodesic in positive even-dimensional curvature yieldsThe deformation need not remain a geodesic or an embedded curve.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 132 3 a Solution Created 2026-10-03 Updated 2026-10-05
A generalized conformal metric is locally , where is nonnegative and measurable and transforms as a length density under a holomorphic change of coordinate. Its area is . Use the usual convention of using locally rectifiable paths for a path family , and put . ThenThis is extremal length. Zeros and isolated singularities of an admissible density are allowed; requiring a smooth strictly positive Riemannian metric would unnecessarily restrict the definition. Line integrals have their extended nonnegative values; if an arbitrary family is supplied, use its members that are locally rectifiable paths. An empty path family has infinite infimal length, whereas a family containing a constant path has extremal length zero.
Both numerator and denominator scale quadratically when is multiplied by a positive constant. The coordinate transformation of the area element makes the quotient unchanged under conformal equivalence.
The normalization relevant later is worth deriving. On the conformal cylinder , let contain the loops going once around it. For the horizontal loop at height , Cauchy-Schwarz inequality givesIntegrating in shows . The constant density achieves equality, since every winding-one loop has Euclidean length at least one. HenceHere is the height divided by circumference, the conformal modulus of an annulus; the reciprocal is used for the family joining its boundary components.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 201 3 c Solution Created 2026-10-03 Updated 2026-10-05
Work on as a probability space, with Lebesgue measure of total mass one. Let be the filtration generated by the level- half-open dyadic cells together with the separate null cell . Define the dyadic slope martingale byThe endpoint may be assigned any value, since it is a null set for Lebesgue measure. The two child slopes average to their parent slope, by telescoping the two increments. Therefore almost everywhere, so this is a martingale. The Lipschitz condition gives everywhere except possibly at the freely chosen endpoint, where we take zero.
Apply the Lp martingale convergence theorem with . Its limit has almost everywhere and in L1 norm as well, by Cauchy-Schwarz inequality. Choose a measurable representative of and set it to zero on any exceptional null set; it is then a bounded measurable function on the entire interval.
Set . Telescoping at the grid points shows that is the linear interpolation of on the dyadic grid. The Lipschitz condition gives . AlsoThe two uniform limits coincide, giving the absolutely continuous function representationThe chosen bound holds for every after the null-set modification; the integral identity holds for every simultaneously.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 202 1 e Solution Created 2026-10-03 Updated 2026-10-05
Put and . The candidate is continuous, adapted, nondecreasing, starts at zero and has integrable terminal value. We prove is a martingale, rather than assuming the desired quadratic variation of a stochastic integral.
The assumed boundedness of also makes a true martingale on every finite horizon. Indeed, apply its local martingale identity at a localizing sequence; bounded stopped squares and the Fatou lemma show that is integrable. The stopped expressions are then dominated by , so the dominated convergence theorem gives the unrestricted conditional identity on that horizon.
First let be a bounded simple predictable process. On each interval with coefficient , an increment of is . The defining quadratic variation identity for , together with the martingale property of , makesa martingale for . Expanding gives a term and the square term above; their conditional expectations yield the desired identity on that interval. Integrability follows because the variables are square-integrable; concatenating the finitely many intervals gives a martingale.
For general bounded predictable , choose bounded simple predictable processes in , which is possible since they generate the predictable sigma-algebra. Set and . The Itô isometry and Doob L2 maximal inequality imply in ; the Cauchy-Schwarz inequality gives in andPass the conditional expectation identity for to the limit. Thus is a martingale. Uniqueness in the defining quadratic variation property now gives
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 203 2 d Solution Created 2026-10-03 Updated 2026-10-05
Using the usual zero-boundary Sobolev space convention,where weak derivatives define and is the space of test functions. In the Gaussian free field convention, the same notation often denotes the Dirichlet energy space, the completion in the gradient norm alone. On bounded domains the Poincare inequality makes the two definitions equivalent; on unbounded domains one must distinguish them. The following gradient-pairing argument applies in either setting whenever the energy completion is realized as weak functions.
Identify with a vector subspace of by zero extension of H01. Approximating by test functions in shows that this is an isometric embedding in the inhomogeneous Sobolev norm and also in the Dirichlet inner product norm. In particular, arbitrary irregularity of causes no additional boundary term.
DefineThese are weakly harmonic Sobolev functions. For , integration by parts in the weak sense gives . If , choose converging to in , or in energy for the homogeneous convention. The Cauchy-Schwarz inequality givesHence the orthogonality of supported and harmonic Dirichlet functions isBoth are linear vector subspaces. In the inhomogeneous convention they are closed: the first is the isometric image of a complete space, and the second is the intersection of the kernels of the linear functionals . No spanning assertion is needed. The PDF contains this orthogonality statement; the TeX has badly corrupted it into an assertion about openness.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 205 2 Solution Created 2026-10-03 Updated 2026-10-05
Use the normalizationfor the Lasso. A centered sub-Gaussian random variable with parameter satisfies for all ; for a noncentered variable apply this definition to . With fixed , independence of the errors givesThe exponential Markov bound, optimized over separately for the two signs, gives . Take and apply the union bound over columns:The specified positive tuning parameter requires ; at its formula is zero. For small the probability lower bound can be negative and hence uninformative. Although centering makes the centered errors dependent, means , so the preceding proof uses the original independent errors, not independence after centering.
The Karush-Kuhn-Tucker conditions, using the subdifferential of the L1 norm, areConsequently, for nonempty , on ,This establishes the required active-set inequality without any rank condition on .
For , a sufficient Compatibility condition for the Lasso isIndeed, for , comparison of the two Lasso objective values and the score bound on give the Basic inequality for the LassoIt follows that satisfies the Lasso cone condition and . Thus , including the trivial case. This proves the slightly stronger , and in particular the requested . If , the same basic inequality forces on ; no nonempty-support compatibility constant is needed.
For the support-size argument, use the stated prediction bound with constant . Set . By the Cauchy-Schwarz inequality and the definition of the sparse maximum eigenvalue,Combining with the lower bound and squaring gives for every nonempty . If and , choose of size ; this contradicts the strict inequality defining . If , proves the assertion directly. Hence in both cases.
When , minimality and monotonicity of the sparse maximum eigenvalues giveWhen , the first result gives , so the same final bound holds. Thus the finite-index conclusion is .
There is a genuine domain omission in the printed last assertion: the definition of only makes sense for , so is undefined when . This case can occur: take centered orthogonal columns with , , , and . Then for every available and none exceeds , so . A universally defined replacement, from and monotonicity, isAlternatively, explicitly extend the definition by for , including infinity. Under that added convention the printed final expression is meaningful and follows from this replacement bound. Neither convention nor finiteness should be silently assumed.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 210 3 Solution Created 2026-10-03 Updated 2026-10-05
Let and be Radon-Nikodym derivatives. Use the following normalizations of the Kullback-Leibler divergence, total variation distance and unnormalized Hellinger distance:We use and assign infinite Kullback-Leibler divergence if , where denotes absolute continuity of measures. Common domination alone also suffices for the inequalities with the extended-value convention. The negative part of the logarithmic integral is integrable: on , . Thus the extended-value integral defining the Kullback-Leibler divergence is well-defined. These expressions do not depend on which common dominating measure is chosen. The Kullback-Leibler divergence is generally asymmetric and therefore fails a defining property of a metric.
For completeness, the integral formula for total variation distance follows by taking : since , the positive and negative parts of have the same integral, and this set attains the supremum.
Put , by the Cauchy-Schwarz inequality. A second application givesThe logarithmic inequality in the hint implies for . Apply it to on :Multiplying by and integrating provesIf on a set of positive -measure, the right side is finite but the Kullback-Leibler divergence is infinite, so the inequality is immediate. ThusThis total variation–Hellinger–relative entropy inequality requires the stated Hellinger distance normalization for its first constant. With the alternative convention , the correct comparison is and . In particular, and have .
Here is a testing form of the Le Cam two-point lemma. Suppose two parameter values in a metric space have separation , and let be the corresponding observation laws. Every estimator satisfiesTo prove it, classify by the closer of the two parameters, breaking a tie in favour of one. Let . On , the triangle inequality implies , while on it implies . Hence the sum of the two probabilities in the display is at leastAt least one is at least half the sum, proving the claim. It also gives a lower bound for the maximum expected metric loss.
For absolute-error loss, or any loss given by a metric, there is a stronger expectation form. Take densities relative to and use the triangle inequality under their overlap:Therefore the Le Cam lower bound under absolute-error loss isBoth arguments are valid for randomized estimators as well, by adjoining the independent randomization to the observation; doing so does not change total variation distance.
For the normal distribution model, assume the usual nondegenerate convention and use the full product observation laws. ChooseFor one observation, its Gaussian log likelihood ratio isUnder the first law, and . Thus the Kullback-Leibler divergence between normal distributions isFor the independent and identically distributed random variables, the joint log likelihood ratio is the sum of the individual ones, so the Kullback-Leibler divergence of the two full statistical samples is . The comparison just proved implies . Applying the expectation form of the Le Cam two-point lemma givesThis is a Gaussian location minimax lower bound under absolute-error loss. The constant may depend on the fixed known , but neither on the estimator nor on . If a signed nonzero scale were used, the same result holds with in place of . If zero noise were allowed, would estimate exactly, so a positive lower bound would be false; positive variance is necessary.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 215 3 b Solution Created 2026-10-03 Updated 2026-10-06
For a finite irreducible Markov chain satisfying detailed balance, define the directed-edge capacity . For every ordered pair choose a path from to using only positive-capacity edges; take the empty path for . Let count occurrences of the directed edge . The canonical paths Poincare bound says thatis a valid Poincare inequality for a reversible Markov chain constant. Simple paths give ; allowing repetitions requires the multiplicities shown.
To prove the theorem, expand the variance using two independent samples from :Telescoping along and applying the Cauchy-Schwarz inequality givesInterchange the finite sums, then use the definition of :Consequently . All edges here are directed, so the in the Dirichlet form of a Markov chain is retained consistently.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 215 4 c Solution Created 2026-10-03 Updated 2026-10-06
The reference kernel printed in the hint has total mass , so it is not a Markov kernel. We explicitly replace it by the normalized random transposition shuffle: choose two positions independently and uniformly and swap them. Its probabilities areThese sum to one, and its ordinary spectral gap is , the value intended in the hint. This value can also be certified by the Aldous spectral gap theorem: the single-label generator has rate between distinct positions, so on centered functions it acts as multiplication by . Both shuffles have the same uniform distribution on a finite set as stationary distribution.
We use the Canonical paths comparison theorem in the following precise form. For two finite reversible Markov chains with the same positive stationary distribution , route each directed transition of along a positive-capacity path of . Ifthen and . Indeed, telescope each difference along its path, apply the Cauchy-Schwarz inequality, and sum with weights . The load definition gives the Dirichlet form of a Markov chain inequality; the common variance denominator then gives the spectral gap inequality. Identity transitions use empty paths.
Put . A transposition is routed using the wordof length . Each generator occurs at most twice; write its occurrence count as . For a fixed directed adjacent edge and each occurrence of in a word, there is exactly one starting permutation whose translated word crosses this edge at that occurrence. This follows because right multiplication by the prefix is a bijection of the symmetric group. Hence the uniform stationary factors cancel in the comparison load, leavingThe crude bounds , , and imply . ThereforeOnly pairs with contribute, and there are of these. Keeping this restriction gives the stronger and . Neither comparison bound needs the exact adjacent-shuffle spectral gap; the normalization repair of the reference kernel is essential.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 217 3 Solution Created 2026-10-03 Updated 2026-10-06
The Cameron-Martin space of a Gaussian random variable in a Banach space is most conveniently defined using the first Gaussian chaos. Work over the real separable Banach space , with continuous dual space . For a centered Gaussian random variable in a Banach space, defineThe set inside the closure is already a linear subspace. Every is a centered Gaussian random variable. The Bochner integral defining exists: Fernique's theorem states that for some , , and then Cauchy-Schwarz inequality gives .
The map is injective. Indeed, implies for all . Density of these observations in then implies . Define the Reproducing-kernel Hilbert space byThe injectivity makes the inner product unambiguous, and the completeness of makes a Hilbert space. The embedding into is continuous sinceFor put . The reproducing property isThis defines the Banach-space Gaussian RKHS even for a degenerate Gaussian measure.
For , write . There is a measurable coordinate on , obtained as a limit in L2 space of continuous linear observations, and has normal distribution . The Cameron-Martin theorem for a Gaussian measure, in its real separable Banach space form, states that if is the probability law of , then the probability law of is an equivalent probability measure to precisely when . For such a shift,For the two laws are mutually singular measures. The coordinate need not be a continuous functional on , and the sample need not belong to .
Apply this theorem to the shift . With , it givesChoose the measurable coordinate so that on a symmetric set of full -measure. This is possible by taking an almost-surely convergent subsequence of the approximating linear observations and intersecting its convergence set with its negative. Alternatively the joint law of is invariant under simultaneous negation, directly from those observations and their mean-square convergence. Since and the centered Gaussian measure is symmetric, the two exponential integrals with signs and agree. Consequentlybecause . The symmetric Gaussian translation lower bound is thereforeOnly symmetry and Borel measurability of are required; convexity is not needed.
There is one degenerate flaw in the last printed claim. The separable Banach space with has dense in , and all positive-radius balls have probability one. Nevertheless for every . Thus the statement needs the additional hypothesis . If the course convention excludes the zero space, this hypothesis is already implicit.
Under this necessary hypothesis, density of implies that contains a nonzero vector. Given , scale that vector to obtain with , and choose . The triangle inequality givesThe centered closed ball is a symmetric Borel set. Apply the symmetric Gaussian translation lower bound to getHence the Gaussian norm distribution is strictly increasing when is nonzero. This proves the strict increase of a Gaussian norm distribution without assuming that its distribution has a density or that spheres have zero probability.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 326 2 4 a Solution 2026-10-06
Write , , and , with as required by coercivity and by the printed inverse power. The Tikhonov regularization with a coercive penalty operator normal equation givesThis also holds for inconsistent data in , because and therefore . Taking the inner product with givesFor any with , put and use Cauchy-Schwarz inequality to bound the right-hand side by . Since , completing the square in yieldsThe positive root of this quadratic inequality is at most . Infimizing over proves the approximate source bound for coercive Tikhonov regularization, in fact slightly more strongly:Injectivity gives . Given any positive tolerance, approximate by for some finite-norm ; for all sufficiently large it is admissible in the infimum. ThusNo attainment of that infimum is needed.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 326 4 3 Solution Created 2026-10-03 Updated 2026-10-06
Because , a finite sum implies . Thus the printed weighted one-norm on l2 equals the extended nonnegative sum on every . Unbounded weights may make that sum infinite for some ; membership in alone is not a finiteness assertion.
For each , is continuous in the l2 sequence space, and . If in norm, then for every ,Taking the supremum over provesIt is also weakly lower semicontinuous by the same finite-coordinate argument, and convex and proper because zero and finite sequences have finite penalty.
At a finite-penalty base point , coordinate variations give the subdifferential of a weighted one-norm on l2:These conditions are sufficient by summing the coordinate supporting inequalities. Since , an infinitely supported would require infinitely many , contradicting . Thus a finite-penalty point can have empty subdifferential; a Bregman base with an actual subgradient must be finitely supported.
For any such , and the Bregman distance simplifies toEach summand is nonnegative; the sum may be infinite. The pairing converges by Cauchy-Schwarz inequality. This definition does not silently assume a subgradient exists at every .
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 339 1 d Solution Created 2026-10-03 Updated 2026-10-06
A pointed cone has no nonzero line through the origin: equivalently . If , its quadratic form vanishes on the nonnegative orthant. Testing the coordinate vectors gives . Testing then givesSince is a symmetric matrix, every entry is zero. Thus .
The identity matrix gives an interior point. For a symmetric perturbation with operator norm , the Cauchy-Schwarz inequality givesThis open norm ball around lies even in the positive semidefinite cone, hence in . Consequently and the interior is nonempty.
More generally, the interior of the copositive cone consists exactly of strictly copositive matrices. Positivity on the compact nonnegative unit sphere has a positive minimum and persists under small perturbations; a zero there is destroyed by an arbitrarily small negative multiple of .
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 3 i Solution Created 2026-10-03 Updated 2026-10-06
The Dirichlet Laplacian eigenfunctions and positive Dirichlet Laplacian eigenvalues on the interval of length two areA Sturm-Liouville eigenfunction expansion gives , with . Thus all normal modes oscillate precisely when .
To justify a spatially uniform bound for arbitrary appropriate data, assume the usual wave equation energy class and , or smoother compatible data. The energy method conservesThe sharp Poincare inequality is . For , this makes coercive with constant ; for , coercivity is immediate. The one-dimensional Sobolev representative satisfies by the fundamental theorem of calculus and the Cauchy-Schwarz inequality. Hence the conserved energy bounds uniformly in both variables.
At , the smooth solution is unbounded. At , is unbounded. ThereforeThe bound depends on the initial data; it is not a single bound for all arbitrarily rescaled solutions. The energy-class hypothesis supplies pointwise meaning and excludes undefined rough-data interpretations.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 5 i Solution Created 2026-10-03 Updated 2026-10-06
Use the real Hilbert space , the zero-boundary Sobolev space, with norm . The Poincare inequality makes this equivalent to its usual Sobolev space norm. Define the symmetric bilinear form and continuous linear functionalThe Cauchy-Schwarz inequality and the sharp interval Poincare inequality implyThus is a bounded bilinear form and a coercive bilinear form, and is bounded because is square-integrable. The Lax-Milgram theorem supplies a unique weak solution:This is obtained from the differential equation by integration by parts; the test functions have zero boundary trace. Conversely, the weak solution satisfies in the sense of distributions. Since , it belongs to and solves the equation almost everywhere.
Equivalently, the Ritz method minimizesover . Indeed , since the mixed term is . This proves existence and uniqueness of the minimizer directly from the weak solution, as well as its equivalence to the variational problem.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 341 7 Solution Created 2026-10-03 Updated 2026-10-06
For a fully discrete linear partial differential equation evolution, write in a specified discrete norm. Finite-time stability of a numerical method meanswith independent of the allowed meshes. For a multistep method, augment the state with its time-history values and include a stable starter. Bounds may grow as when the continuous problem has growth; demanding decay for all time would be a stronger assertion. The treatment of initial data, forcing, Dirichlet boundary conditions, and the norm is part of the hypothesis, not a detail supplied by an interior calculation.
For constant coefficients on an infinite or periodic uniform grid, Von Neumann stability analysis substitutes the Fourier mode . The discrete Fourier transform and the Parseval identity turn a uniform multiplier estimate into a discrete L2 norm estimate. In a one-step scalar scheme, proves contractivity, while gives finite-time stability. In a multilevel scheme, every root of a polynomial in its amplification polynomial of a multilevel finite difference scheme matters, as does uniform control of the associated companion matrix.
For the Forward Euler diffusion scheme with ,Thus for all frequencies exactly whenSufficiency follows since ; necessity follows from the highest-frequency mode on even periodic grids. This gives the usual parabolic mesh restriction . For the Backward Euler diffusion scheme,so every is stable. The Crank--Nicolson method hasagain of modulus at most one for all , but stiff modes approach rather than zero. These are PDE manifestations of the A-stability and L-stability distinctions for time integration.
For advection , assume and set . The upwind finite difference scheme is , withIt is contractive when . This also has a direct maximum-norm proof: each newest value is a convex combination of two old values. The resulting Courant–Friedrichs–Lewy condition expresses that the numerical domain of dependence covers the physical one. For negative , the upwind direction must be reversed. By contrast, Forward Euler method time stepping with the centered first finite difference has and for nonzero . Under fixed nonzero refinement it is unstable, since a fixed nontrivial mode grows geometrically over steps. This is not a claim of instability under every imaginable mesh coupling: instead bounds its spurious finite-time growth, since .
The leapfrog advection scheme illustrates a multilevel subtlety. Its amplification polynomial of a multilevel finite difference scheme isFor , both roots of a polynomial have unit modulus and separation at least . A uniformly conditioned eigenbasis of the companion matrix then proves power boundedness of a two-level Fourier scheme for arbitrary history data. At and a grid admitting , the roots coincide on the unit circle, and a Jordan block produces growth. Thus the endpoint fails the ordinary arbitrary-history root condition for a multistep method; checking only that both moduli equal one misses the instability. Restricting the starter or filtering the parasitic mode is an additional hypothesis.
The energy method handles variable coefficients and finite boundaries for which Fourier stability analysis may be unavailable. For the Backward Euler diffusion scheme with homogeneous Dirichlet boundary conditions, define . The discrete summation by parts identity isTake the inner product of with . The elementary identity givesThis proves unconditional contractivity including the actual boundary treatment. More generally, a dissipative operator gives of operator norm at most one: with , dissipativity implies , and the Cauchy-Schwarz inequality gives . In finite dimensions this also proves invertibility. For the second-order backward differentiation formula, the BDF2 discrete energy identity used in question 4 controls both time levels and proves unconditional diffusion stability. A repeated root strictly inside the disk is harmless here; the energy argument supplies the uniform bound without a singular eigenvector formula.
A second technique is eigenvalue stability analysis of a finite difference method. Under the method of lines, a time integrator advances by . If is a normal matrix, the scalar linear stability domain criterion on all controls its powers exactly in the corresponding L2 norm. For the centered Dirichlet discrete Laplacian, its eigenvalues lie between and zero; intersecting this interval with a time integrator's stability interval determines its mesh restriction. A uniform bound on the diagonalization of a matrix is needed for nonnormal diagonalizable systems. Eigenvalues alone can be misleading. For example,has both eigenvalues inside the disk for , but the upper-right entry of is . Taking , , this grows like . Hence stable scalar eigenvalues do not give mesh-uniform stability. Boundary closures can introduce precisely the extra growth that an interior Fourier symbol overlooks, so boundary stability of a finite-difference method must be checked separately.
For nonlinear spatial discretizations, a useful replacement for Fourier stability analysis is a convexity argument. If the Forward Euler method map is nonexpansive in a chosen norm for , the second-order strong stability preserving Runge-Kutta methodis also nonexpansive under the same restriction. Indeed, for two inputs, the triangle inequality gives . Its Taylor expansion is , proving order two. Such strong stability preserving Runge-Kutta methods transfer suitable forward-step bounds without relying on a linear spectral calculation.
Finally, consistency of a numerical method explains what stability buys. For a linear Hadamard well-posed problem, the Lax equivalence theorem equates convergence of a consistent discretization with its stability, under the stated approximation-space and norm hypotheses. Directly, if the error satisfies , iteration gives the discrete Duhamel principleThus vanishing normalized local truncation error and initial error yield convergence. With a nonlinear Lipschitz step estimate, the discrete Gronwall inequality plays the same role. A stable but inconsistent stencil, such as the literal defective formulas in questions 3 and 4 under their usual refinement, is not rescued by any amplification bound. The decisive checks are a mesh-uniform evolution bound, a valid boundary closure, and consistency with the actual PDE.
Past exam of the mathematics course of the University of Cambridge 2017 ii Paper 2 34D b iii Solution Created 2026-09-24 Updated 2026-10-05
Differentiating the grand potential givesA reversible adiabatic change in an isolated gas keeps and fixed. To infer constant , one must check that this ratio is one-to-one, rather than merely assuming it.
Put , , and for . Integration by parts gives and , andThe strict inequality is the Cauchy-Schwarz inequality for the positive weight , since is not constant on its support. This proves strict monotonicity of Fermi gas entropy per particle and hence is constant. Then gives constant, and givesThis adiabatic scaling of a Fermi gas with power-law density of states concerns positive temperature; the zero-temperature limit is a separate limiting case.