Matrix mechanics Created 2026-10-18 Updated 2026-10-20
Matrix mechanics represents an observable by an array of matrix elements between chosen quantum states. In an orthonormal basis of energy eigenstates, each a stationary state, the diagonal entries of the Hamiltonian operator are the energy levels; off-diagonal entries of other observables describe couplings between those states. The same array can therefore answer two different questions: which values can a measurement return, and which pairs of energy eigenstates can a specified interaction connect? The eigenvalues answer the first question; the relevant matrix elements, together with the energy differences, answer the second.
This makes linear algebra a practical language for atomic physics. One can calculate energy levels, electric-dipole selection rules and probability amplitudes for transitions without assigning an electron an unobserved classical orbit. The three-state example below separates a missing optical coupling from a missing energy level; the finite oscillator commutator boundary defect explains why a finite array requires care when approximating a quantum harmonic oscillator.
Metric embedding Created 2026-09-28 Updated 2026-10-14
A metric embedding represents the points of one metric space inside another by an injective function . The geometric question is how faithfully the representation preserves pairwise distances. For example, binary strings can be treated as vectors so that tools from linear algebra become available, but the chosen norm determines what their separation means. A drawing of the points alone does not specify a metric embedding: both metrics matter.
An isometric embedding preserves every distance exactly. A metric embedding that is a bi-Lipschitz map allows multiplicative error while keeping its bi-Lipschitz distortion bounded over all pairs. A coarse embedding instead controls separation through possibly nonlinear functions of distance. These are different requirements: the Hamming cube as an L1 and L2 metric below preserves Hamming distance exactly in one representation and takes its square root in another.
The usefulness is that a complicated distance comparison becomes a calculation with familiar vectors and norms. This also helps assess a visualization: nearby points should not acquire arbitrarily large separations, and far-apart points should not collapse together. Matoušek's 2013 lecture notes describe how motivations from computer science, especially approximation algorithms, accelerated the study of metric embeddings from the 1990s onward. [ref]
Module theory Created 2026-09-24 Updated 2026-10-14
Module theory studies abelian groups with scalar multiplication by a ring. It extends linear algebra: a module over a field is a vector space, while ring scalars let the same language describe integer arithmetic and linear operators. The new phenomena are visible before any abstract classification: some nonzero scalars cannot be divided out, some modules have torsion elements of a module, and some have no basis of a module.
Figure 1.
Free integer coordinates and a cyclic module with torsion
. The free module has two independent integer directions. In the cyclic module , repeatedly adding returns to zero after six steps: . The circle records addition by the generator; its geometry is only a drawing convention. Details.
A historical anchor is Emmy Noether's 1921 paper [ref]. Its first section studies rings and ideals under the condition that every ideal is a finitely generated module over its ring, and connects this to the ascending chain condition; later sections develop primary decomposition. This is an early source for the structural viewpoint underlying Noetherian rings. [ref] makes those statements accessible in English.
Required applied partial differential equations and perturbation methods; choose two or three numerical-analysis cores from numerical solution of partial differential equations, numerical linear algebra and continuous optimisation. At least two special-topic units, including modelling/methods and computation. One scientific-computing case study in Michaelmas term and one mathematical-modelling case study in Hilary term. Compulsory four-unit dissertation, normally written in Trinity term and the summer.
Why compare inferred neutrino fluxes rather than raw event counts? For channel , exposure time , target population , cross section , and detection efficiency including reconstruction and selection, a forward model is
Here is expected background; on the left counts events, while counts targets. Distinct cross sections, target numbers, thresholds and detection efficiencies prevent a direct CC/NC count ratio from being an electron neutrino fraction. Channel populations themselves can be correlated estimates from overlapping event patterns.
An original reduced linear algebra model isolates the inference. Assume specified spectral shapes and divide each background-corrected channel rate by its own calibrated electron neutrino response. Put and let be the spectrum- and response-weighted non-electron/electron ES sensitivity ratio. Then
The equal and responses determine their sum, not their separate neutrino fluxes. Changing the assumed spectra generally changes the integrated coefficients; with arbitrary unknown spectra this small matrix is insufficient.
Choose illustrative inputs, not historical measurements: , and , with every neutrino flux in the same arbitrary unit. The forward prediction is . CC and NC independently invert to , ; substituting into ES returns . Alternatively, CC and ES give . Thus the third channel tests the response model and the inferred mixture.
Figure 1.
Chosen vacuum survival curve and complementary neutrino flux constraints
. Original calculations only. Left: stipulated maximal two-flavour mixing and squared-mass splitting, with chosen paths at one GeV. Right: exact normalized-channel equations for chosen fluxes and ES sensitivity; lines are not experimental confidence bands. Details.
CC alone admits any nonnegative . CC and ES lose the second independent constraint when ; for small positive , an ES-minus-CC error is amplified by . With chosen independent errors , their propagated standard deviation is . If CC and NC instead have those independent errors, the difference has standard deviation . Actual analyses must retain covariances, including those introduced by statistical channel separation. Complementarity supplies information; calibration establishes what information each channel carries.
Newton's identities 2026-10-14
Newton's identities let two different summaries of the same polynomial roots talk to each other. Elementary symmetric polynomials record products of distinct root positions; power sums record sums of individual powers, . Either can be easier to obtain in a particular problem. The identities turn one into the other, often avoiding any need to find the individual polynomial roots.
Write and for . For every positive integer , the identity is
It is a polynomial identity over the integers, so it remains true on substitution in any commutative ring. To solve it for by division, however, must have a multiplicative inverse. Over a field of characteristic zero, the first power sums recover every elementary symmetric polynomial. The same recovery works when the positive characteristic of a field is greater than ; small characteristic needs care.[ref]
The first three relations are , , and . When division by two and three is permitted, this becomes
For the three roots , their power sums are , , and . The formulas give and , recovering the polynomial from aggregate data.
Figure 1.
Three power sums reconstruct a cubic without assigning names to its roots
. Left: exact contributions of the stipulated roots to the first three power sums. Middle: Newton's identities recover , , . Right: the resulting monic polynomial has its three roots of a polynomial at the original values. These are algebraic examples, not inferred experimental data. Details.
The beauty is that the summaries do not depend on the order of the polynomial roots. In symmetric polynomials, one can work with quantities that survive every relabeling; in linear algebra, traces of matrix powers give the same kind of aggregate information about eigenvalues. The identities convert those summaries into coefficients of a characteristic polynomial.
A short derivation packages the identities into a formal power series. Let . Since its constant coefficient is one, it has an inverse as a formal power series, and its formal derivative satisfies
Multiply by and compare the coefficient of . This gives Newton's identities without a convergence assumption or division by . The division issue arises only when using them to recover . The logarithmic-derivative approach is also developed in Newton's Identities (January 2019 notes) by Milan Mossé (§4, pp. 6–8).
The alternate name Newton-Girard identities recalls Albert Girard's earlier work. A modern historical account credits formulas through sixth powers to Albert Girard in 1629 and discusses their later appearance in Isaac Newton's Arithmetica universalis.[ref] This citation refers to that account; the original historical editions have not been inspected here.
For three skew lines and a point outside them, let be the plane through and . A line through meeting all three must lie in all three . Therefore the matrix of their normal vectors must have zero determinant. This necessary condition can fail to be sufficient: the common line of the planes can be parallel to, and disjoint from, one of the original lines.
For example, take , , and . The corresponding planes have equations , , . Their common line is , which never meets . This illustrates why a linear algebra condition for a common direction must be supplemented by intersection checks.
Solutions to all three questions of the University of Oxford Part A linear algebra paper, sat on 15 June 2021. Source: released examination PDF, parent PDF pages 253–254; paper pages 1–2. All questions and parts are covered, including both alternatives in questions asking whether an assertion is true.
Solutions to all three questions of the University of Oxford Part A linear algebra paper, sat on 20 June 2023. Source: released examination PDF, parent PDF pages 258–260; paper pages 1–3. Question labels follow the paper; the solutions supply their own reasoning. Part 2(c) includes the exceptional case of characteristic of a field three because that part does not specify its field separately.
Use unrestricted long and short portfolios, as required by the usual arbitrage pricing theory argument. A vector of initial monetary positions with costs zero. Its terminal financial payoff, in the exact one-factor model, is
If also , the financial payoff is the constant . A nonzero value would be an arbitrage after choosing the sign of . Absence of arbitrage thus says that annihilates . By finite-dimensional linear algebra,
Write . Taking expectations, assuming the factor has finite mean, proves
This is exact factor pricing without a traded risk-free asset. The PDF prints an asset-indexed ; the stronger valid result has the same for every asset, so it also satisfies the printed relation by taking all its equal. If every loading is equal, the spanning vectors are dependent and the coefficients need not be unique. For a nondegenerate factor and a zero-exposure unit-cost portfolio, is that portfolio's certain return; if a risk-free asset is explicitly traded, law of one price makes it the risk-free asset's return. No equilibrium preferences are needed for this argument.
For a rational matrix with columns , partition regularity means that every finite colouring of the positive integers admits a positive monochromatic vector with . Repeated coordinates are allowed. The columns condition is an ordered partition into nonempty blocks such that
Rado's theorem says that is partition regular if and only if it satisfies the columns condition. We prove both directions.
For necessity, clear denominators so the columns are integer vectors. For every pair of disjoint index sets with for which is not in the rational linear span of the I-columns, choose an integer linear functional that vanishes on those I-columns but not on that sum. Such a functional exists by finite-dimensional linear algebra; clear its rational denominators. There are only finitely many pairs, so choose a prime number dividing none of the nonzero integers
If there are no offending pairs, any prime number suffices.
Colour a positive integer , , by . This is the last nonzero digit coloring, using colours. Partition regularity gives a positive monochromatic solution. Group its coordinate indices by their distinct P-adic valuations , writing . All units have the same nonzero residue modulo .
For each , let . If its block sum failed the required span condition, apply the chosen to . The earlier terms vanish exactly, not just modulo a power of . Divide by ; terms with later valuations are still divisible by . Reduction modulo then yields
contrary to the choice of . Thus every required span relation holds, and with the first block sum is zero. This is the finite separating-functional proof of the columns condition; it avoids an unjustified passage from congruences to rational equality.
For sufficiency, assume the columns condition. Choose rational coefficients , for and , such that
Choose a positive integer clearing every denominator, and a positive integer , taking if the list is empty. By the allowed monochromatic m-p-c set theorem, there are positive generators such that all numbers
are positive and have one colour. This is a full m-p-c set; positivity is part of its definition, not a rule discarding negative expressions. For , put
These coordinates all belong to that monochromatic m-p-c set. Collecting coefficients of gives
Thus the Rado solution inside an m-p-c set proves partition regularity and completes the theorem.
To obtain the Finite sums theorem from Rado's theorem, introduce a positive variable for every nonempty and impose
Partition the columns of this system into . The -column sum is zero: in a row indexed by containing , the contributions of and cancel, and other rows receive neither. For , the -column sum is
where is the column of . Those columns lie in earlier blocks. Hence the columns condition holds. A monochromatic positive solution gives and , proving
For the system has no rows and the conclusion is immediate. The construction is also visible directly in the sufficiency proof: a monochromatic -set contains every nonempty sum of its generators, classified by the smallest chosen index.
The Frankl-Wilson theorem in its uniform form says: let be a prime number, let have size , and let have , while for every distinct pair of members. Then
The modular-intersection, nonuniform form of the Frankl-Wilson theorem is as follows. Let be a prime number, let have elements, and let satisfy for each , whereas for distinct . Then
This nonuniform Frankl-Wilson theorem allows different sizes and different diagonal residues; in particular it applies when all sizes have one common residue outside .
To prove it by the polynomial method in combinatorics, work over the finite field and associate to each member the intersection polynomial
Replace each positive power of in its expansion by to obtain its multilinear reduction on the Boolean cube . This preserves evaluation on every characteristic vector, whose coordinates are zero or one, and leaves degree at most . At the characteristic vector ,
It is zero for distinct , and nonzero for by the diagonal hypothesis. If , evaluating at each forces . Thus the multilinear polynomials are linearly independent. They belong to the vector space with basis the square-free monomials for , whose dimension is the displayed sum. This proves the nonuniform form, including , where distinct members are impossible and the empty product is one.
For the sharp uniform form, we need the low-degree evaluation rank on a uniform layer. Form the integer matrix whose rows are indexed by with , columns by -subsets , and entries . These are evaluations of the square-free monomials. Over the rational numbers, the row of a -subset satisfies
Indeed, when the left side counts the -subsets between and , and otherwise both sides vanish. Since , the coefficient is a positive integer. Thus the degree- rows span all the rows over the rational numbers, giving . Every square submatrix of larger size has integer determinant zero, so its determinant remains zero modulo . Hence
This reduction of matrix rank is the reason one must not try to divide by directly in .
The evaluations of the on the th uniform layer of the Boolean cube lie in the row linear span of over . They are linearly independent as functions there, because restriction to the columns corresponding to still gives the nonzero diagonal evaluation matrix. Their number is consequently at most . This completes the proof of the uniform Frankl-Wilson theorem as well.
For the odd-cardinality problem the sharper Oddtown theorem follows directly from linear algebra. Let be the characteristic vector. The parity conditions say
For any linear combination , take its dot product with . All off-diagonal terms vanish, leaving . Hence these vectors are linearly independent in an -dimensional vector space, giving
The singleton sets show that the bound can be attained. Notice why the constant term in the general nonuniform Frankl-Wilson theorem is unnecessary here: these dot products already give a diagonal evaluation matrix using only the coordinate functions.
For even , partition into the pairs and take every union of a selection of these pairs. This gives exactly distinct sets, including the empty set. Each has even size; every pairwise intersection is again a union of entire pairs and therefore has even size. This attains the proposed bound.
For any set family satisfying the even-cardinality conditions, let be the linear span of its characteristic vectors over . Now the self-dot products vanish as well as the off-diagonal ones. By bilinearity, for all , so : it is a self-orthogonal binary subspace.
For completeness the standard bilinear form on is nondegenerate, because a vector orthogonal to every coordinate vector has every coordinate zero. Thus the map from onto the dual space is surjective: every linear functional on extends to the whole vector space by extending a basis, and every such functional on is a dot product with a vector. Its kernel is . The rank-nullity theorem consequently gives . The inclusion implies . Finally, distinct members have distinct characteristic vectors, all in , and a -dimensional vector space over has exactly elements. Therefore
This proves the Eventown theorem and also explains the suggested dimension contradiction if the set family had more than members.
A rational matrix is partition regular over the positive integers if every finite colouring of admits a vector whose coordinates all have the same colour and satisfy . Equal coordinates are allowed. Rado's theorem says that this holds exactly when the columns can be partitioned into nonempty blocks with
This is the columns condition. Clearing an overall denominator lets us prove necessity for an integer matrix.
Suppose the columns condition fails. There are only finitely many ordered partitions of the finite column index set. For each ordered partition, choose a failing block: its sum is outside the rational span of preceding columns, with that span interpreted as zero for the first block. Finite-dimensional linear algebra gives a rational linear functional vanishing on that span but not on . For completeness, take a basis of the span, append , extend to a basis of the column space, and prescribe values zero on the first basis vectors and one on . Clearing the functional's denominators gives an integer row functional with .
Choose one prime larger than the absolute values of all these finitely many nonzero integers . Colour each positive integer by its first nonzero base- digit:
Assume a monochromatic solution exists, and partition its coordinates into blocks of equal P-adic valuation, ordered by increasing valuations . Their first nonzero digit is a common . For this ordered partition take its chosen failing block and functional . Apply to . All preceding-block terms disappear exactly. Divide the remaining integer equation by and reduce modulo . Later blocks disappear, while the current block gives
This contradicts the prime's choice. Thus the colouring has no monochromatic solution, proving necessity. This is the finite separating-functional proof of the columns condition; it needs neither a limiting argument nor a bound on the valuations themselves.
For sufficiency, suppose the blocks satisfy the columns condition. Choose rational coefficients , for in preceding blocks, such that
Let be a positive integer clearing their denominators, and choose an integer bounding all . Use the allowed monochromatic m-p-c set theorem with generators. We use the triangular convention in which the resulting positive integer set contains every number
in one colour. This is an m-p-c set, with the generator indices reversed if the alternative lower-triangular convention is used. The permitted theorem applies to either convention, since it applies to every finite number of generators. Its set lies in , so all of the displayed combinations are positive.
For , define
Each coordinate belongs to that monochromatic set. Summing the column contributions and collecting coefficients of gives
For the inner sum is empty and . This proves sufficiency by a Rado solution inside an m-p-c set, and completes the theorem.
For the requested application, choose
The matrix, in the coordinate order , has columns
The block has zero sum. The remaining column satisfies , so it is in the span of the first block. By Rado's theorem, positive monochromatic satisfy the system. In particular,
One can see the positivity directly in the same triangular construction: a monochromatic set containing for , together with , gives
All four are positive because they belong to that positive m-p-c set. This is a partition-regular system forcing an ordered Schur relation.
Write the variables as , with . The hypothesis makes invertible. The implicit function theorem supplies neighbourhoods of and a unique smooth map such that
Its derivative is . Thus the derivative of this parametrization has image exactly .
Since has a constrained maximum at , differentiating gives zero in every direction. Therefore vanishes on . In finite-dimensional linear algebra, is the row space of . Consequently
The full-rank hypothesis is what permits this Lagrange multiplier conclusion; it is not a consequence of constrained maximality alone at a singular constraint.
The regular University of Oxford preliminary examination was held on 20 June 2019, from 14:30 to 17:30. These solutions cover every marked part of all eight questions: digital systems in Questions 1–3, linear algebra in Questions 4–6, and natural deduction in Questions 7–8, including questions beyond the examination's permitted choice. The supplied Thumb instruction set reference and primitive proof rules are part of the held paper.
Source release: Oxford Computer Science examination archive, member 2018/T/a10101w1.pdf. The native cover establishes the 2019 regular sitting despite the archive's 2018/T directory. The printed code is A10101W1. A correction to a printed invariant bound is explained in Question 1(e); all inputs needed for the remaining solutions are present.
Quantum mechanics Created 2026-09-24 Updated 2026-10-18
Quantum mechanics describes a quantum system using a quantum state, with observables represented by self-adjoint operators. For a single particle in one dimension, a pure state can be described by a wavefunction : its squared modulus gives the probability density for position. The Schrödinger equation predicts how this probability amplitude evolves between interventions; measurement in quantum mechanics connects it to experimental outcomes through the Born rule. A wavefunction therefore predicts a distribution of possible positions rather than a single trajectory as in classical mechanics.
The infinite square well makes three central ideas visible with elementary sine functions. Confinement selects discrete energy eigenvalues; wavefunction normalization converts a shape into quantitative probabilities; and a quantum superposition can have a changing probability density even when each constituent energy eigenstate has a stationary one. This combination explains why quantum mechanics uses both linear algebra and differential equations: eigenfunctions describe the building blocks, while quantum phases determine how they interfere. The same tools help explain the discrete energies of chemical atoms and molecules, although their Quantum Hamiltonians require different potentials and often many particles.
In his first 1926 paper, Erwin Schrödinger formulated quantization as an eigenvalue problem for a hydrogen atom in nonrelativistic quantum mechanics. On its opening page he compared the emergence of integers with the number of nodes of a standing wave on a vibrating string. The infinite square well below realizes that appealing idea especially transparently: the allowed integers follow from the boundary conditions. This historical analogy concerns the origin of discrete energies; the Born rule supplies the statistical interpretation of the modern wavefunction. [ref].
Qubit Created 2026-09-25 Updated 2026-10-14
A qubit, or quantum bit, is a quantum system whose state space is a two-dimensional complex Hilbert space. Choose an orthonormal basis , called the computational basis. A pure state has a normalized representative , where are complex probability amplitudes and . General qubit states, including mixed states, require a density matrix. [ref].
The abstract qubit describes the chosen degree of freedom: for spin one-half, can label spin up and down along a specified axis; for photon polarization in a fixed propagation mode, they can label horizontal and vertical polarization, . These implementations share the same two-dimensional linear algebra, but their preparation and measurement in quantum mechanics use different apparatus. A qubit therefore need not be an entire two-level particle. [ref].
The usefulness of a qubit lies in controlling quantum superpositions and their relative quantum phases, then converting them into outcome probabilities with quantum gates and an appropriate measurement. One ideal computational-basis measurement returns one label, zero or one. It does not reveal both probability amplitudes, their relative quantum phase, or all the bits of a larger register's expansion. Quantum state tomography instead estimates a qubit state using many repetitions of the preparation and suitable measurement settings. [ref].
An actual preparation and measurement example is Dederick and Beck's 2014 teaching-laboratory experiment. A 405-nm laser pumped two nonlinear crystals, producing signal and idler photons near 810 nm. In their non-entangled run the source produced ; a quarter-wave plate in the signal arm changed its photon polarization to in their convention. The plate supplies a relative quantum phase between polarization components. An idler detection selected the signal trials; signal detectors measured horizontal/vertical, diagonal, and circular polarization bases on separate runs. Before the polarization analyzer, a half-wave plate at degrees selected the diagonal basis, and a quarter-wave plate at degrees selected the circular basis. Coincidence counts gave the conditional positive-outcome frequency . The reconstructed density matrix for conditioning on detector had squared quantum fidelity with and purity of a density operator : close to the intended pure state, rather than an exact ideal preparation. [ref].
Figure 1.
Ideal polarization-qubit preparation and measurements in three separate settings
. Original model schematic and exact Born-rule probabilities, not measured data or a reproduction of the experimental apparatus. The circular-basis sign uses the explicitly defined state . Details.
The model diagram explains why those settings matter. With and , the intended gives equal probabilities in both the horizontal/vertical computational basis and diagonal Hadamard basis, just as the mixed state does. In the circular basis , however, the ideal always gives the plus outcome, whereas still gives equal probabilities. These are deductions from the Born rule, not the paper's measured counts. Unequal detector efficiencies and finite samples must be accounted for when interpreting experimental frequencies; the ideal comparison assumes equally efficient detection.
A linear subspace is self-orthogonal if for every . The standard bilinear form is nondegenerate, whence and . The condition with means every vector has even Hamming weight. Such a binary linear code gives the linear algebra behind the Eventown theorem.
Sylvester matrix Created 2026-10-14 Updated 2026-10-18
For polynomials of polynomial degrees , the Sylvester matrix has shifted coefficient rows of followed by shifted rows of , using descending-power order. It is a square matrix of size . Its transpose represents the linear map with , ; its determinant is the resultant.
For an original parameter example, choose and . Then
The two lists of coefficients become three rows because is multiplied by both and . Thus the common-root parameter values are and . At , the resultant is , so the two polynomials have no common root.
Figure 1.
A Sylvester matrix tests when a moving root of a polynomial coincides with a fixed one
. Original example , . Left: the vertical dashed line marks the root of a polynomial of ; it does not make zero. Right: the same test for all parameters is the function graph of the resultant . A crossing between two arbitrary function curves would not suffice: both values must be zero. Details.
One can verify the certificate without any determinant expansion:
For , division by gives a Bézout's identity for polynomials. At , the identity instead becomes a nontrivial cancellation between multiples of two polynomials sharing the polynomial factor .
The general proof needs only a short linear algebra outline. A common polynomial factor gives a nonzero cancelling pair within the degree bounds, making the Sylvester matrix singular. Conversely, if are coprime polynomials and , then divides ; the bound forces and then . Thus the linear map has zero kernel of a linear map exactly in the coprime case. Over an algebraic closure, a nonconstant common polynomial factor is equivalent to a common root.
In 1840, James Joseph Sylvester described forming the shifted coefficient rows and taking a signed sum of products selected by permutations. In modern notation, this is the Sylvester matrix and its determinant. The construction turns the nonlinear common-root of a polynomial question into linear algebra. His paper's word derivatives concerns elimination, rather than the formal derivative used below.[ref]
Topological vector space Created 2026-09-24 Updated 2026-10-14
A topological vector space is a vector space over or , equipped with a topology such that vector addition and scalar multiplication are jointly continuous:
The domains carry the product topology, and the scalar field carries its usual topology. Joint continuity permits both inputs to vary at once. Here the definition does not impose the Hausdorff condition; we state it separately when needed. [ref].
Every normed vector space is a topological vector space: the triangle inequality gives , while
controls scalar multiplication near any . More generally, a seminorm topology supplies many measurements of smallness, rather than one norm. Translation by any vector is a homeomorphism, so a neighbourhood basis at zero determines all neighbourhoods by translation.
The compatibility matters: on a nonzero real vector space, the discrete topology makes vector addition continuous, but not scalar multiplication. For , the sequence of scalars tends to zero, whereas the sequence never enters the open set . Thus choosing a topology arbitrarily need not respect the linear algebra.
This framework lets functional analysis choose the appropriate notion of limit of a sequence. Norm topology measures one global error, weak topology tests linear functionals, and the Fréchet space of smooth functions controls every order of derivative. In each case the operations forming linear combinations are continuous, making approximation compatible with the underlying equations.