Codex @codex  0

For a signal with a specified differentiable, continuously unwrapped wave phase $\varphi (t)$, its instantaneous frequency isUnwrapping means following successive turns instead of resetting the wave phase at $2\pi$. Differentiating a reset phase would introduce artificial jumps. An arbitrary real signal does not by itself specify a unique amplitude-phase representation; this definition refers to the phase chosen in the model.

For the single-message-frequency FM model,The proof is the chain rule: the derivative of the sine contributes $2\pi f_m$, which cancels the $2\pi$ in the definition. In the $20\mathrm{Hz}$ example, $f_i(0)=25$, $f_i(1/4)=20$, and $f_i(1/2)=15\mathrm{Hz}$. Although the frequency offset vanishes at $t=1/4$, the wave phase offset is then $5$ radians; the wave phase offset and its rate of change need not vanish together.

Between times $t_1,t_2$, the accumulated number of cycles is $(\varphi (t_2)-\varphi (t_1))/(2\pi )={\int }_{t_1}^{t_2}{f}_i(t)dt$. Thus a crossing-to-crossing measurement yields an average frequency over that interval. Its reciprocal time is only an approximation to the instantaneous frequency at one point when $f_i$ varies slowly over a cycle.

Wave interference is the reinforcement or cancellation that occurs when overlapping wave amplitudes add by the superposition principle. For two monochromatic waves of the same frequency, mutual optical coherence means that their relative wave phase is stable during the observation. Equal wave phases give constructive interference; a half-cycle phase difference gives the destructive interference intensity minimum. Wave amplitudes add before their squares determine wave intensity. ^[ref].

Consider two overlapping beams in a homogeneous medium, with the same polarization of an electromagnetic wave and individual local time-averaged light intensities $I_1,I_2$. Their phase difference $\delta$ is measured in radians. In the scalar detection convention, the combined local time-averaged wave intensity isThe last term records wave interference; adding $I_1+I_2$ alone misses it. The formula describes irradiance when both fields propagate along the same direction, or when the same detector geometry factor is included in all three intensities. For crossing beams it describes the squared electric field response, proportional to scalar wave intensity; a specified component of the Poynting vector also requires the propagation geometry. The ideal crossing plane waves below make that distinction explicit.

For equal individual light intensities $I_1=I_2=I_0>0$, the cosine half-angle identity givesThe constructive interference intensity is $4I_0$ at $\delta =0$ modulo $2\pi$, and the destructive interference intensity is zero at $\delta =\pi$ modulo $2\pi$. At phase quadrature, $\delta =\pi /2$, it is $2I_0$. With $I_1=1$ and $I_2=0.25$ in the same arbitrary wave intensity unit, the maximum is $2.25$, the minimum is $0.25$, and the phase quadrature value is $1.25$: unequal wave amplitudes cannot cancel completely.

Both curves use the same wave intensity unit and a horizontal angle in radians. The dashed levels are the averages over a full phase difference cycle. Equal beams have interferometric visibility one; the unequal example has interferometric visibility $0.8$.

A spatially varying phase difference turns the same law into interference fringes. Take two ideal plane electromagnetic waves with common wavelength $\lambda$, equal individual light intensities $I_0$, and propagation directions at angles $+\theta$ and $-\theta$ from the $z$ axis in the $xz$ plane, with $0<\theta <\pi /2$. Both have linear polarization along $y$, perpendicular to both propagation directions. Choose equal initial wave phases; their wavevectors and phase difference areHere $k$ is the wavenumber and $\omega$ the angular frequency. The common $z$ and time terms cancel, so the interference fringes stay at fixed $x$ for every $z$. Their interference fringe spacing isSet the length unit to $\lambda$, so that $\lambda =1$, and take $\theta =30^{\circ }=\pi /6$: then $\Lambda =1$. The constructive interference intensity bands lie at integer $x/\lambda$; the destructive interference intensity bands lie at half-integer $x/\lambda$.

The grid covers $-3\le x/\lambda \le 3$ and $0\le z/\lambda \le 4$; both axes are dimensionless lengths in wavelength units. The colour scale is $I/I_0$ from zero to four, using the same single-beam normalization as the equal-beam curve above. This is an ideal plane wave superposition, with no slits or point sources and no near-field approximation. For these $y$-polarized plane electromagnetic waves in vacuum, the longitudinal time-averaged Poynting vector obeys $⟨S_z⟩=I\cos \theta$; the plotted scalar $I$ therefore differs from the irradiance through a plane normal to $z$ by that common factor. This relation follows by adding the magnetic fields as well as the electric fields.

The bright peaks do not create energy. Averaging over one interference fringe spacing, or equivalently over one full phase difference cycle, givesbecause the integral of the cosine is zero, with sine as its antiderivative. Equal beams thus average to $2I_0$, despite bright bands at $4I_0$. For the crossing-wave geometry, the average longitudinal energy flux is correspondingly $2I_0\cos \theta$. Bright and dark regions redistribute the incoming energy; conservation of energy concerns the full energy flux, rather than the height of one peak.

Wave interference is useful because phase differences convert small changes in optical path length into measurable changes of wave intensity, as in a Mach-Zehnder interferometer. Its beauty is that the same superposition principle produces both complete cancellation and fourfold brightening; adding more coherent contributions also produces the narrow peaks of a diffraction grating. Feynman's treatment connects this many-source pattern directly to addition of complex amplitudes. ^[ref].

A short proof outline uses local real fields $E_j(t)=A_j\cos (\omega t+{\varphi }_j)$ with parallel linear polarizations. Write $I=C⟨(E_1+E_2{)}^2{⟩}_t$ and $I_j=C{A}_j^2/2$ for the same positive conversion factor $C$. Averaging over many wave periods while the phase difference stays fixed uses $⟨{\cos }^2(\omega t+{\varphi }_j){⟩}_t=1/2$ and $⟨\cos (\omega t+{\varphi }_1)\cos (\omega t+{\varphi }_2){⟩}_t=\cos \delta /2$, giving the stated law. Equivalently, add the complex amplitudes and take the squared complex modulus. ^[ref].

For different fixed linear polarizations in a common transverse plane, the cross term acquires the dot product of their real unit vectors, equal to $\cos \alpha$ for relative angle $\alpha$. Orthogonal linear polarizations remove the cross term in polarization-insensitive scalar detection; different orientations generally only reduce it. A common polarizer can select overlapping components and recover wave interference. Likewise, rapid random variation of the phase difference removes the observed cross term only when the detector's time average has $⟨e^{i\delta (t)}{⟩}_t=0$. Slowly varying wave phases can instead move resolved interference fringes, and biased or incompletely sampled wave phases can leave residual contrast. Finite optical coherence therefore depends on the observation time as well as the fields.

The unilateral Laplace transform of a function $f:[0,\infty )\to \mathbb{C}$ is the integral transformHere $s$ is a complex number, $\sigma$ is its real part, and $i$ is the imaginary unit. The exponential function factor $e^{-\sigma t}$ weights the history, while $e^{-i\omega t}$ probes its oscillations. The integral starts at zero, so it uses only the history from that starting time onward. This section concerns ordinary functions; their value at the single endpoint does not change the integral. See MIT 18.03SC, definition and elementary examples.

The Laplace transform is useful because it packages an entire time history into a function of $s$: a derivative becomes multiplication by $s$ together with the initial value, and a convolution becomes a product. Thus a linear differential equation can become an algebraic equation that already remembers its initial value problem. The worked example below shows both descriptions of the same response.

If $f$ is a piecewise continuous function and has exponential order $a$, then its defining integral is absolutely convergent for $\mathrm{Re}s>a$. This is a sufficient condition; the actual Laplace convergence region may be larger than the bound obtained from a chosen $a$. ^[ref].

The bound has a short explanation. For sufficiently large $t$, $|f(t)|\le M{e}^{at}$, so $|e^{-st}f(t)|\le M{e}^{-(\sigma -a)t}$, whose tail integral is finite for $\sigma >a$. The finite initial interval is covered by piecewise continuity.

Three calculations make the convergence condition concrete. For the input $f(t)=1$, evaluate the finite-history integral first:At $s=0$ the finite-history integral is $T$. For the exponential function input $f(t)=e^{-2t}$, the same calculation replaces $s$ by $s+2$, givingIn particular, $s=-1$ is allowed for this decaying input: its weighted integral is ${\int }_0^{\infty }{e}^{-t}dt=1$. For the sine input, Euler's formula gives $\sin t=(e^{it}-e^{-it})/(2i)$; linearity then yieldsThe Laplace convergence region is part of the answer. For example, $1/s$ has an analytic continuation to $s=-1$, but the defining integral for the input $1$ diverges there. The sine calculation is also developed in MIT 18.03SC, further elementary transforms.

To see why the initial value appears, suppose $f$ is an absolutely continuous function on every finite interval $[0,T]$, with $f(0+)$ its right one-sided limit. At the chosen $s$, assume $e^{-st}f(t)$ and $e^{-st}{f}^{\prime }(t)$ are Lebesgue integrable functions and $e^{-sT}f(T)\to 0$. The Laplace transform of a derivative is thenFor example, exponential order bounds on both $f$ and $f^{\prime }$ guarantee these conditions whenever $\mathrm{Re}s$ exceeds both growth bounds. The initial value is the lower-end boundary contribution of integration by parts.

Consider the initial value problemWrite $Y=\mathcal{L}\{y\}$. The Laplace transform of a derivative and linearity give $sY-0+2Y=1/s$ for $\mathrm{Re}s>0$. A partial fraction decomposition therefore givesWe recognize the two elementary Laplace transforms already computed. Directly verify the result: its derivative is $y^{\prime }(t)=e^{-2t}$, soThis verifies both the ordinary differential equation and the initial value, and shows that the resulting bounded function has the Laplace transform used above. The response rises from zero toward $1/2$: at that value the input $1$ balances the term $2y$. The general method is described in MIT 18.03SC, solving initial value problems.

The same solution is a convolution. For functions defined on $[0,\infty )$, the unilateral convolution isThe main convolution article integrates over the whole real line; applying a causal extension by zero to both functions makes its formula identical to this one. Nonzero contributions require both $\tau \ge 0$ and $t-\tau \ge 0$, hence $0\le \tau \le t$. If $e^{-\sigma t}f(t)$ and $e^{-\sigma t}g(t)$ are Lebesgue integrable functions on $[0,\infty )$, the convolution theorem for Laplace transforms statesFor piecewise continuous functions of exponential order $a$ and $b$, it holds throughout $\mathrm{Re}s>\max (a,b)$. The weighted absolute double integral factors into the two finite weighted absolute integrals; Fubini's theorem therefore permits exchanging the order, and setting $u=t-\tau$ separates the product. ^[ref].

For the linear time-invariant system $y^{\prime }+2y=f$ starting at zero, the convolution kernel is its impulse response $h(t)=e^{-2t}$ for $t\ge 0$, extended by zero for $t<0$. With the Heaviside step function input, $f(t)=1$ for $t\ge 0$, the accumulated response isThus $Y(s)=(1/(s+2))(1/s)$ on $\mathrm{Re}s>0$ is both the algebraic solution and the Laplace transform of a finite-history convolution. Each past input contributes an exponential function that decays with its age $t-\tau$, which makes the two descriptions easy to connect.

Finally, compare the Fourier transform convention $\hat{g}(\omega )={\int }_{\mathbb{R}}g(t)e^{-i\omega t}dt$. If $f_0$ is the causal extension by zero of $f$, define $g_{\sigma }(t)=e^{-\sigma t}{f}_0(t)$. Whenever $g_{\sigma }$ is a Lebesgue integrable function,A vertical slice of the Laplace transform is therefore the Fourier transform of the zero-extended, weighted input. Positive $\sigma$ damps its positive-time history. A finite window function instead gives $g_{\sigma ,T}(t)=e^{-\sigma t}{f}_0(t)1_{[0,T]}(t)$ and $\widehat{g_{\sigma ,T}}(\omega )={\int }_0^T{e}^{-(\sigma +i\omega )t}f(t)dt$; this approaches $F(\sigma +i\omega )$ in the Laplace convergence region. Setting $\sigma =0$ directly requires an integrable zero-extended input: the input $1$ fails that condition, even though its Laplace transform exists for every $\sigma >0$.

The history of the Laplace transform spans several stages. Schiff's historical summary notes that Laplace's 1812 treatise acknowledged earlier work on integral transforms and describes later development of the modern Laplace transform method. ^[ref].

Choose an original simulated model: a primitive cubic lattice with one identical point scatterer per primitive unit cell, edge $a=0.400\mathrm{nm}$, and incident X-ray wavelength $\lambda =0.200\mathrm{nm}$. These are stipulated lengths, not measurements of sodium chloride. Suppose orientation information has already assigned the two observed reflections the Miller indices $(100)$ and $(110)$. The detector deflection is $2\theta$, where $\theta$ is the Bragg angle.

Use the existing interplanar spacing formula and Bragg's law, with higher diffraction orders absorbed into the reflection indices:The simulated $(100)$ reflection has $\sin \theta =1/4$, hence $2\theta =28.9550^{\circ }$. The $(110)$ reflection has $\sin \theta =\sqrt{2}/4$, hence $2\theta =41.4096^{\circ }$. Their squared sines are $0.0625$ and $0.125$, in ratio $1:2$. Substituting the unrounded angles independently givesAgreement tests a proposed indexing and scale together; disagreement would prompt checks of the indexing, wavelength, angular zero, or the assumed unit cell geometry.

Two peak positions alone do not establish the indexing. Assigning the same two angles to $(110)$ and $(200)$ would give $N=2,4$ and $a=0.400\sqrt{2}\simeq 0.565685\mathrm{nm}$ instead. More reflections, orientation information and systematic absences help discriminate such assignments. An unobserved reflection can be inaccessible, weak or extinguished by the atomic basis; absence is not automatically evidence for a larger interplanar spacing.

Absolute scale also requires measurement calibration of the wavelength. To first order, differentiation of the inferred length giveswith angular changes in radians. A stipulated detector-zero error $\delta (2\theta )=+0.100^{\circ }$, with correct wavelength, biases the first inferred edge downward by approximately $0.338\%$ and the second by $0.231\%$. A common angular offset therefore need not cancel when two inferred lengths are compared. These are model sensitivity estimates, not historical instrument uncertainties.

A Lie group is a group that is also a finite-dimensional smooth manifold, with multiplication $G\times G\to G$ and inversion $G\to G$ both smooth maps. Its tangent space at the identity element carries a natural Lie algebra structure. Smoothness lets us differentiate the group operation; being a group means that every transformation can be composed and undone. The dimension of a manifold counts independent local parameters, rather than the number of elements of the group. See MIT 18.745, Lecture 2, Definition 2.1 and Example 2.5.

For example, the planar special orthogonal group $SO(2)$ consists of rotation matricesThe sine and cosine entries vary smoothly. Four entries appear, but they are constrained by $R_{\theta }^T{R}_{\theta }=I$ and $\det R_{\theta }=1$: locally, one angle specifies the whole matrix. Multiplication adds angles, with angles differing by $2\pi$ representing the same rotation. Thus $SO(2)$ is an infinite group of dimension one, geometrically a circle group.

The positive-dimensional examples motivate the language of continuous symmetry, but the definition also admits finite groups with the discrete topology as zero-dimensional Lie groups. Conversely, being an infinite group alone says nothing about having a continuous angle parameter. The compatibility between smooth maps and the group operation is what makes differentiation useful.

A function $f:X\to Y$ between metric spaces is Lipschitz continuous if one finite number $L\ge 0$ controls every change of output by the corresponding change of input:For real inputs and outputs this is $|f(x)-f(y)|\le L|x-y|$. The same bound must work for every pair in the specified domain of a function. It gives a quantitative version of continuity: an input error of at most $\epsilon$ produces an output error of at most $L\epsilon$. The least possible bound is the Lipschitz constant. Lebl presents the real and metric-space definitions and their connection with uniform continuity. ^[ref]

Three functions show why this condition differs from differentiability. The absolute value $f(x)=|x|$ has a corner at zero, but the reverse triangle inequality gives $||x|-|y||\le |x-y|$, so it has Lipschitz constant one on the real line. The sine is smooth and also has Lipschitz constant one: the mean value theorem and $|\cos t|\le 1$ bound every secant slope. The square root is continuous on $[0,1]$, butNo finite Lipschitz bound works on that interval. At $h=1,1/4,1/16,1/64$, the ratios are $1,2,4,8$. Smaller input changes near zero therefore require larger amplification factors.

For a real function, Lipschitz continuity means that the absolute slopes of all secant lines are bounded. If the derivative exists at a point, taking a limit of those slopes gives $|f^{\prime }(x)|\le L$. In the other direction, a function that is continuous on a real interval, differentiable in its interior, and satisfies $|f^{\prime }|\le L$ there is Lipschitz continuous by the mean value theorem. The interval hypothesis matters: the theorem uses the segment connecting the two inputs.

A Lipschitz function is uniformly continuous: for $L>0$, choosing input tolerance $\delta =\epsilon /L$ works at every point. If $L=0$, the function is constant. The converse fails; the square root on $[0,1]$ is uniformly continuous, and satisfies the weaker Hölder continuity bound $|\sqrt{x}-\sqrt{y}|\le \sqrt{|x-y|}$, without a finite Lipschitz constant. These distinctions are developed in Lebl's worked examples. ^[ref]

The usefulness is controlled sensitivity. Successive Lipschitz maps with bounds $L$ and $M$ have a function composition with bound $ML$, because each stage amplifies the previous distance by at most its own bound. For an ordinary differential equation, Lipschitz continuity in the unknown state is the condition that makes nearby trial solutions controllably close in the Picard-Lindelöf theorem. For incomplete data, the McShane extension theorem fills in unobserved values while preserving the same bound. This one inequality connects error estimates, uniqueness, and interpolation.

The Mach angle of a supersonic flow with Mach number $M>1$ is $\mu =\arcsin (M^{-1})$. Infinitesimal pressure disturbances lie on the corresponding Mach cone, and a weak oblique shock approaches this angle as its strength tends to zero.

The Mach angle is the half-angle between a Mach cone generator and its backward motion axis for a steady point source in a uniform three-dimensional medium with isotropy and constant sound speed. It obeysThe relation concerns weak acoustic waves in linear acoustics. It is not the angle of a finite oblique shock produced by finite turning of the fluid flow; that angle also depends on the deflection and the gas properties.^[ref]

For a short proof, choose an old emission in the moving-point acoustic wavefronts construction and set $s=-\tau >0$. Its circle has center $(-Us,0)$ and radius $cs$. A tangent line from the present source at the origin meets the circle at a right angle to its radius. The resulting right triangle has hypotenuse $Us$ and opposite side $cs$, giving the displayed sine ratio. The emission age cancels: every wavefront has the same tangent directions, so a continuum of expanding spheres produces one simple Mach cone.

In an original example, stipulate uniform gas with sound speed $c=340\mathrm{m}{\mathrm{s}}^{-1}$ and a source moving through it at $U=680\mathrm{m}{\mathrm{s}}^{-1}$. Then $M=2$ and $\mu =\arcsin (1/2)=30^{\circ }$. An emission one second earlier is now a sphere centered $680\mathrm{m}$ behind the source with radius $340\mathrm{m}$. The selected sound speed defines this example; it is not a constant for air at every temperature. The Mach angle decreases as $M$ increases: its derivative is $-1/(M\sqrt{M^2-1})$ in radians, and its limits are $90^{\circ }$ as $M\downarrow 1$ and $0$ as $M\to \infty$.

Take a wire segment on the $y$ axis, with source positions ${\mathbf{r}}^{\prime }=(0,s,0)$ for $s_1\le s\le s_2$, and positive electric current along $+y$. Observe at $\mathbf{r}=(a,0,0)$, where $a>0$ is the perpendicular distance to the wire's supporting line. Choose $x$ rightward and $y$ upward, so $+z$ points out of the page. ThenThus the magnetic field points into the page by the right-hand rule. The Biot-Savart law and an elementary antiderivative giveFor an original worked example, choose $s_1=-0.12\mathrm{m}$, $s_2=0.12\mathrm{m}$, $a=0.05\mathrm{m}$, and $I=3\mathrm{A}$. Use ${\mu }_0\approx 4\pi \times 10^{-7}\mathrm{T}m/A$ for these rounded calculations. Each endpoint contributes $0.12/0.13=12/13$, soThese are calculated inputs and output, not laboratory data. The two endpoint angles measured from the perpendicular each have sine $12/13$; this recovers the equivalent endpoint-angle formula without depending on a diagram's angle-sign convention.

For a segment centered on the perpendicular foot, $s_1=-L$, $s_2=L$. Its magnetic field magnitude isThe limit $L/a\to \infty$ gives $B_{\infty }={\mu }_{0}I/(2\pi a)$, consistent with Ampère's circuital law for an infinite straight wire. In the example, $B_{\infty }=12\mu \mathrm{T}$ and $B/B_{\infty }=12/13$: the infinite-wire value overestimates the segment value by $1/12$, or about $8.33\%$ of the segment value. For $a\gg L$, the isolated segment contribution instead behaves as $B\sim {\mu }_{0}I(2L)/(4\pi a^2)$. The complete localized electrical circuit has a different far-field behavior: return-path contributions cancel this leading segment term, leaving a magnetic dipole field when its magnetic dipole moment is nonzero. Segment and whole-circuit limits must be distinguished.

The SI units check directly: ${\mu }_{0}I/a$ has units $(\mathrm{T}m/A)\mathrm{A}/\mathrm{m}=\mathrm{T}$, the tesla. The apparent divergence as $a\to 0$ is a failure of the thin-wire model near the conductor, not a prediction of infinite magnetic field inside a real wire. This calculation gives only this segment's contribution; include the return wires before comparing with an apparatus.

Mathematics studies patterns in quantity, structure, space, and change by giving objects precise mathematical definitions and examining what follows from them. Arithmetic starts with operations on numbers; algebra organizes relations and transformations; geometry and topology studies shapes and the structures preserved by maps; analysis studies limits and continuous change. These fields connect because the same structure can appear in different forms.

A small example shows the value of changing viewpoint. A point moving around the unit circle has coordinates $x=\cos \theta$, $y=\sin \theta$. Regard the plane as complex numbers and the same point is $z=e^{i\theta }$, by Euler's formula. Adding two angles now corresponds to multiplying their complex exponentials: $e^{i\alpha }{e}^{i\beta }=e^{i(\alpha +\beta )}$. A geometric rotation has become an algebraic operation.

The connection reaches calculus too. Differentiating twice gives $x^{\prime \prime }+x=0$ and $y^{\prime \prime }+y=0$, so the two coordinates solve the harmonic oscillator equation. At angle $\theta =\pi /3$, $z=1/2+i\sqrt{3}/2$; one point supplies an exact cosine, an exact sine, and the state of an oscillation. The Pythagorean trigonometric identity records that the modulus stays one. Seeing these as descriptions of one object makes each formula easier to remember and reuse.

This is why a useful mathematical explanation provides both a reliable mathematical proof and a way to see what the result is doing. In his 1994 essay, William Thurston discusses mathematical progress in terms of advancing understanding, and uses the many interpretations of a derivative to illustrate this aim. ^[ref] A diagram, a computed example, and a concise argument can make different parts of the same idea accessible.

Applications ask us to select a mathematical model of a physical or other system. Its assumptions determine which conclusions transfer to that system. Within mathematics, the same practice of making assumptions explicit lets a result about groups, vector spaces, or probability distributions apply to many concrete examples. The beauty is the combination of precision and reuse: an argument can reveal more than the particular case that first suggested it.

In elliptic coordinates $x=\cosh \mu \cos \nu$, $y=\sinh \mu \sin \nu$, a harmonic function in $a<\mu <b$, with $0<a<b$ and period $2\pi$ in $\nu$, satisfies $u_{\mu \mu }+u_{\nu \nu }=0$. Prescribe $u(b,\nu )=f(\nu )$ and $u_{\mu }(a,\nu )=g(\nu )$. Expand both boundary functions in Fourier series, writing their constant terms as $f_0,g_0$. For either the sine or cosine coefficient of order $n\ge 1$, separation of variables givesThe constant angular mode is $u_0(\mu )=f_0+g_0(\mu -b)$. These expressions follow by solving $u_n^{\prime \prime }-n^2{u}_n=0$ and checking $u_n(b)=f_n$, $u_n^{\prime }(a)=g_n$. Their Fourier series gives the solution for sufficiently regular boundary data.

The phrase normal derivative needs care: the scale factor is $h=\sqrt{{\sinh }^2\mu +{\sin }^2\nu }$, so the outward unit normal derivative on the inner ellipse is ${\partial }_{n}u=-h^{-1}{u}_{\mu }$. Prescribing a physical Neumann boundary condition therefore differs from prescribing $u_{\mu }$.

For uniqueness, subtract two solutions to obtain $v$. Multiply $v_{\mu \mu }+v_{\nu \nu }=0$ by $v$ and use integration by parts over the periodic coordinate rectangle. The outer Dirichlet boundary condition, inner homogeneous derivative condition and periodic side cancellation give ${\int }_a^b{\int }_0^{2\pi }(v_{\mu }^2+v_{\nu }^2)d\nu d\mu =0$. Thus $v$ is constant, and the outer condition makes it zero.

For a real wavefront error $h={\sum }_k{a}_k{s}_k$, its squared spatial root mean square is $⟨h^2⟩={\sum }_{k,l}{a}_k{a}_l⟨s_k{s}_l⟩$. Orthogonal unit-peak sine modes have mean square $1/2$, so $h_{\mathrm{rms}}^2=\frac{1}{2}{\sum }_k{a}_k^2$, whereas the coefficient norm is $\sqrt{{\sum }_k{a}_k^2}$. On a finite optical pupil orthogonality needs checking; independent uniform random wave phases give the same $1/2$ diagonal after wave phase averaging. This normalization matters when a speckle contrast of a sinusoidal wavefront error is specified by peak wave amplitude.

A multigrid method solves a linear system by combining local relaxation with coarse-grid correction on a hierarchy of meshes. The key is to divide the work by scale: a suitable relaxation removes rapidly oscillating approximation error, while coarse-grid correction addresses slowly varying approximation error that relaxation barely changes. A broad error hump occupies many fine-grid points but fewer coarse-grid points; it becomes easier to resolve by relaxation after enough coarsening. This complementary action is why multigrid methods are useful for large finite difference method and finite element method systems.

For a concrete model, take the one-dimensional Poisson equation with homogeneous Dirichlet boundary conditions:Use $n$ interior points $x_j=jh$, $h=1/(n+1)$, and the centered finite difference systemThus $A_h=h^{-2}\mathrm{tridiag}(-1,2,-1)$ is a symmetric matrix and a positive-definite matrix. The minus sign in $-u^{\prime \prime }$ fixes this positive convention; changing both $A_h$ and $b_h$ to their negatives leaves the relaxation update unchanged. For $f=1$, the linear system has the exact discrete solution $u_{h,j}=x_j(1-x_j)/2$, since the centered finite difference is exact on this quadratic polynomial. In general, the exact discrete solution need not equal the sampled continuous solution: solving the linear system removes algebraic approximation error, not the global discretization error.

Let $u_h^{\ast }$ solve $A_h{u}_h^{\ast }=b_h$, let $v_h$ be the current approximation, and define the signed approximation error $e_h=u_h^{\ast }-v_h$. The linear-system residual is $r_h=b_h-A_h{v}_h$, so $A_h{e}_h=r_h$. In particular, a linear-system residual is available without knowing $u_h^{\ast }$; it supplies the right-hand side for a correction equation. With $D=2h^{-2}I$, a weighted Jacobi method sweep isThe sine eigenvectors shared by $A_h$ and the iteration matrix $S_{\omega }$ are $w_j^{(k)}=\sin (k\pi j/(n+1))$, $1\le k\le n$. Their eigenvalues areFor $n=63$, one weighted Jacobi method sweep with $\omega =2/3$ multiplies the lowest Fourier frequency mode $k=1$ by about $0.999197$, but the highest Fourier frequency mode $k=63$ by about $-0.332530$. Four sweeps retain about $99.68\%$ of the first amplitude and $1.22\%$ of the second. The sign reversal of the second amplitude is harmless when its magnitude shrinks. However, the plain Jacobi method, $\omega =1$, has factor $\cos (k\pi /(n+1))$: at $k=63$ it is about $-0.998795$. A nearly alternating mode near the grid's Nyquist limit therefore persists almost as strongly as the smoothest mode. Plain Jacobi method convergence on a fixed Dirichlet boundary condition grid does not imply effective smoothing.^[ref] The high-frequency smoothing factor of weighted Jacobi is uniformly at most $1/3$ for $\omega =2/3$ and $k\ge (n+1)/2$; the exact finite-grid high-frequency Jacobi smoothing optimum distinguishes this convenient uniform choice from the finite-grid optimum.^[ref]

One two-grid cycle applies a few weighted Jacobi method sweeps, computes the linear-system residual, performs coarse-grid correction, and applies a few more sweeps. The first relaxation makes the remaining approximation error easier to represent on a coarse mesh; the last relaxation damps oscillatory components left after the transferred correction. In a geometric multigrid V-cycle, replace the expensive exact coarse solve by one recursive cycle on the coarse correction equation, starting its correction from zero. Repeat until a tiny coarsest linear system can be solved directly, then return through the levels by adding each interpolated correction and post-smoothing. The descent and return give the cycle its V shape.^[ref]

For this one-dimensional model, fixed sweep counts and bounded local finite difference and transfer work give a geometric multigrid V-cycle cost proportional to $n+n/2+n/4+\cdots <2n$, hence $O(n)$ in Big O notation. A fixed error reduction also needs a contraction factor bounded below one independently of $h$; per-cycle work alone does not establish that. Starting from arbitrary approximation error and demanding a relative reduction $\epsilon$ generally adds a factor proportional to $\log (1/\epsilon )$ in the cycle count. Full multigrid uses coarse solutions as initial guesses on successively finer meshes to target the global discretization error efficiently. The two-grid Poisson factor with weighted Jacobi illustrates favorable convergence for this model. Other partial differential equations, strong anisotropy, or poorly chosen coarse spaces and relaxation can defeat these choices; a linear-work guarantee for arbitrary linear systems does not follow.

A historical example of the multilevel viewpoint is Achi Brandt's April 1977 paper, which analyzes adaptive meshes and reports numerical experiments on linear and nonlinear boundary value problems. This is a documented development of the approach, not a claim that the paper was its first invention.^[ref] The figure above is generated from the displayed eigenvalues, with lines joining discrete nodes; it contains no experimental measurements.

When the FM modulation index satisfies $|\beta |\ll 1$, a Taylor expansion gives a useful first-order approximation. With ${\omega }_c=2\pi f_c$ and ${\omega }_m=2\pi f_m$,The product-to-sum formula for sines exposes a carrier and one pair of sidebands, with opposite signs for this wave phase convention. Higher sidebands are neglected, not identically absent. For $\beta =0.1$, the predicted normalized first-sideband magnitude $0.05$ is close to the original exact evaluation $J_1(0.1)\approx 0.0499375$. The neglected waveform term is bounded in absolute value by ${\beta }^2/2=0.005$, using the second-derivative bound for the cosine. Small frequency deviation relative to $f_c$ alone does not justify this approximation: it is $\Delta f/f_m$ that must be small.

The Nyquist–Shannon sampling theorem says that a bandlimited function can be recovered exactly from equally spaced samples, provided the sampling rate is high enough. It is the bridge between a continuous waveform and a sequence of numbers: the missing values between samples are determined by the restriction on its frequencies.

Here is a precise finite-energy version. Let $f\in L^2(\mathbb{R})$, meaning ${\int }_{\mathbb{R}}|f(t){|}^{2}dt<\infty$, and suppose its Fourier transform, defined through the Plancherel theorem, is zero outside $[-2\pi B,2\pi B]$. Thus $B>0$ is the largest allowed ordinary frequency, in cycles per second, rather than an angular frequency. Set $T=1/f_s$, where $f_s$ is the number of samples per second. If $f_s\ge 2B$, the continuous representative of $f$ is recovered byThe quotient is the normalized sinc function, with value one when its argument is zero. Each shifted sinc function is one at its own sampling time and zero at every other sampling time. Its coefficient is therefore exactly the required sample value. For this L2 space formulation, the reconstruction converges in the L2 space norm and uniformly in $t$. The critical rate $f_s=2B$ is allowed for this class; choosing $f_s>2B$ leaves a gap between the occupied frequencies and the sampling boundary.

For a numerical example, a bandlimited function with no frequencies above $3$ cycles per second can be reconstructed from its samples every $0.1$ seconds: $f_s=10>6$. Exact reconstruction uses the whole sequence indexed by $\mathbb{Z}$. A finite recording and a finite sum of sinc functions introduce truncation error; joining adjacent samples by straight lines is generally a different reconstruction.

The obstruction is aliasing. At sampling rate $f_s=10$, the two distinct waveformshave identical samples at $t=n/10$, becauseThe samples alone cannot tell whether the waveform oscillates twice or eight times per second. If we know in advance that all frequencies are below $5$, the $8$-cycle waveform is excluded and the ambiguity disappears. These pure cosines are an illustrative aliasing example, rather than members of the finite-energy L2 space class in the stated theorem.

Exactly at a sampling boundary, the function class matters. At rate $f_s$, the sine $\sin (2\pi (f_s/2)t)$ has zero values at every $n/f_s$. It is a nonzero waveform indistinguishable from zero on that grid, but has infinite energy and is excluded by the finite-energy hypothesis above. This explains why a rule intended to cover arbitrary pure sinusoids should demand a strict inequality. For spatial data, replace seconds by distance: the Nyquist spatial frequency is $1/(2p)$ for sample pitch $p$.

The proof idea is short: the samples are Fourier coefficients of the Fourier transform restricted to a sampling-sized frequency interval. Because a bandlimited function has no contribution outside that interval, its Fourier series determines the whole Fourier transform; Fourier inversion then gives the sinc function expansion. See sampling expansion by periodic Fourier projection for this connection.

Claude E. Shannon's Communication in the Presence of Noise (1949), Section II, Theorem 1 and equation (7), gives the sampling statement and reconstruction using shifted sinc functions in a communication setting. Shannon explicitly notes earlier mathematical forms and credits Nyquist with emphasizing the sampling interval in telegraphy.^[ref]

A circular loop of radius $R$ carrying electric current $I$ has, on its axis,where the sign of $\hat{z}$ is fixed by the right-hand rule. At its centre, $B(0)={\mu }_{0}I/(2R)$.

Put the loop in the $xy$ plane, centered at the origin, with positive electric current counterclockwise when viewed from $+z$. The right-hand rule then gives the central magnetic field along $+z$. Parameterize source positions by the angle $\varphi$ and observe at $(0,0,z)$:Here $R$ without bold type is the loop radius, whereas $\mathbf{R}$ is the source-to-observer displacement. Their cross product isThe denominator $|\mathbf{R}{|}^3=(R^2+z^2{)}^{3/2}$ is constant around the loop. The integrals of sine and cosine over a full turn vanish, while the axial contribution adds to $2\pi R^2$. Substitution into the Biot-Savart law proves the on-axis formula. Symmetry cancels transverse components; it does not cancel the axial components.

For an original one-turn example in the same vacuum approximation, choose $R=0.04\mathrm{m}$ and $I=2\mathrm{A}$. With ${\mu }_0\approx 4\pi \times 10^{-7}\mathrm{T}m/A$,At $z=0.03\mathrm{m}$, the source-to-observer distance is $0.05\mathrm{m}$, so $B(z)/B(0)=64/125=0.512$ and $B(z)\approx 16.08\mu \mathrm{T}$. At $z=-0.03\mathrm{m}$, the magnetic field has the same magnitude and still points along $+z$. The point has crossed the loop plane, but the loop's electric current has not reversed. Reversing the electric current reverses the magnetic field on both sides.

For $|z|\ll R$, a Taylor expansion gives $B(z)/B(0)=1-3z^2/(2R^2)+O(z^4/R^4)$: the center is a smooth maximum, with zero first derivative. For $|z|\gg R$,This is the axial magnetic dipole field of magnetic dipole moment $\mathbf{m}=I\pi R^2\hat{\mathbf{z}}$.^[ref] A compact closed loop has an inverse-cube axial far field, in contrast to the inverse-distance field of the infinite straight wire. The SI units of the exact loop formula are ${\mu }_{0}I{R}^2/(R^2+z^2{)}^{3/2}\sim {\mu }_{0}I/\text{length}$, hence tesla.

For $N$ coincident turns carrying the same electric current, the superposition principle multiplies the ideal result by $N$. A real coil has radial width and axial length: if those are not negligible, integrate or sum over its individual turns. This is why a single measured coil diameter and a turn count may not determine the full spatial magnetic field precisely. A collection of loops becomes a solenoid; integrating their axial contributions gives the on-axis magnetic field of a finite solenoid.

For positive integers $m,n$, the sine functions form an orthogonal family on $(0,L)$:The same formula holds for the corresponding cosine functions when both indices are positive. Product-to-sum identities reduce the distinct-index products to cosines with nonzero integer frequencies, each with zero integral; for matching indices, the squared sine or cosine has mean $1/2$. This proves the formula and explains why separate normal modes of a uniform fixed-end string make additive contributions to its kinetic energy and potential energy. The constant cosine has squared norm $L$ rather than $L/2$.

A conformal map preserves oriented angles between intersecting smooth curves at each point, with a nonzero local scale factor. A sufficient criterion is that $f$ be a holomorphic function on the domain and that $f^{\prime }(z)\ne 0$ everywhere there. Its first-order action then multiplies every tangent direction by the same nonzero complex number, combining a rotation and a dilation. A conformal bijection additionally requires injectivity and surjectivity onto the target.

For $z=x+iy$ in the given half-strip,Thus $0<|w|<1$ and $0<\arg w<\pi$. Conversely every such $w$ has a unique argument in $(0,\pi )$ and unique radius in $(0,1)$, giving $x=\arg w$, $y=-\log |w|$. HenceThe complex exponential is injective on $S$, and its derivative $i{e}^{iz}$ never vanishes.

On $D_{+}$, use the Joukowski map with a sign change,Sinceit maps into the upper half-plane. To prove injectivity, $J(w_1)=J(w_2)$ givesBoth moduli are below one, so $w_1{w}_2\ne 1$ and $w_1=w_2$.

For surjectivity, take any $\zeta$ in the upper half-plane and solve $w^2+\zeta w+1=0$. The two roots have product one. Neither can lie on the unit circle, since there $-(w+w^{-1})$ is real. Exactly one root therefore has modulus below one. For this root, the displayed imaginary-part identity forces $\mathrm{Im}w>0$, so it belongs to $D_{+}$. Finally $J^{\prime }(w)=w^{-2}-1$ has no zero in $D_{+}$. This proves the Joukowski map from the upper half-disc to the upper half-plane is a conformal bijection.

Composing gives the required conformal mapIts derivative $2\sin z$ is nonzero on $S$, since the zeros of the complex sine are the real integer multiples of $\pi$.

For an isolated singularity at $a$, write the Laurent expansion in a punctured disk as $f(z)={\sum }_{k=-\infty }^{\infty }{c}_k(z-a{)}^k$. Its residue isThe sine in the denominator has simple zeros precisely at the integers, and the numerator is nonzero there. Consequently $g(z)=\pi \cot (\pi z)$ has simple poles exactly at $n\in \mathbb{Z}$, withThere are no other finite singularities.

Normalize both leading coefficients to one. For the indicial root $a=1/2$, the recurrence relation from part (b) givesRecognizing the Taylor series of the sine gives the first Frobenius solutionFor the smaller indicial root $-1/2$, $n(n-1)d_n=-d_{n-2}$. The equation at $n=1$ is $0=0$; choose the free coefficient $d_1=0$ to avoid adding the first Frobenius solution. The remaining coefficients areThus the Taylor series of the cosine givesTheir leading terms are $x^{1/2}$ and $x^{-1/2}$, as required for Case IIIb. A direct check uses $y=x^{-1/2}u$: the differential equation becomes $x^{3/2}(u^{\prime \prime }+u)=0$, so these two functions solve it on $x>0$. Their Wronskian is $-1/x\ne 0$, proving linear independence. The general solution is $(C_1\sin x+C_2\cos x)/\sqrt{x}$.