Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Cm Math Verification

1 · Prerequisites

2 · Summary

The verification layer of the non-relativistic classical mechanics collection, all of it mathematics: the second variation of a time-dependent action and its oscillator saddle, the non-uniqueness of the zero-endpoint Dirichlet problem, the measure-theoretic parallel-axis theorem, the explicit normal modes of two coupled oscillators, and the finite-time small-angle error bound for the pendulum. Every item carries domain: mathematics and depends only on mathematical suppliers; nothing here asserts that a physical system is described by these models, that any unit convention holds, or that any parameter is measured. The classical-mechanics readings belong to the consuming physics pages, which add their own postulates and empirical qualifications.

Second variation of an action. The fixed-endpoint second variation of a time-dependent action at a curve defines the fixed-endpoint family q+εη and the second derivative δ2S(q)[η,η]; The second variation of a time-dependent action: twice differentiability of the fixed-endpoint family and the Hessian integral formula proves that the family is twice differentiable and identifies the second derivative with the integral of the Hessian quadratic form of the Lagrangian along the family, so the definition is a well-defined real number. Both items sit over the chart and first-variation conventions of A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations and The first variation of a time-dependent action on a finite interval: admissibility, the integral formula and continuity of the Euler-Lagrange coefficient.

The oscillator saddle and Dirichlet endpoint degeneracy. Second variation of the oscillator action on a finite interval: both signs for T>π/ω and the value at T=π/ω computes δ2S(q0)[ηn,ηn]=T2((nπ/T)2−ω2) for the variations ηn=sin⁡(nπt/T) of the zero curve of L=12(v2−ω2q2) and shows that for T>π/ω both signs occur among arbitrarily small variations, so the stationary zero curve is neither a local minimiser nor a local maximiser, while at T=π/ω the second-variation value in direction η1 is zero. More generally, The undamped oscillator: Dirichlet non-uniqueness at T=nπ/ω despite unique initial-value solutions proves that the Dirichlet problem with zero endpoints has nonzero solutions exactly for T=nπ/ω with an integer n≥1; in each such case its solution space is {Asin⁡(ωt):A∈R}, although every initial-value problem is uniquely solvable.

Parallel-axis theorem. The inertia tensor of a finite bounded Borel measure about any point of R3 extends the measure-theoretic inertia matrix of Mass distributions on R3: total mass, first moment and inertia matrix to an arbitrary reference point, Finiteness and a uniform bound for the entries of the inertia tensor about an arbitrary point proves the entries finite with an explicit uniform bound under ACω, and Parallel-axis theorem: the inertia tensor about a displaced point of a centred mass distribution proves the centred identity Iν(a)=I(ν)+m(∥a∥22I3−aaT), so that the translation behaviour that the parent definition left open is settled for centred mass distributions.

Two coupled oscillators. Two equal masses with two outer springs and a coupling spring: eigenmodes, decoupling and the explicit normal frequencies treats K=(k+c−c−ck+c): it identifies the orthonormal eigenvectors u± with eigenvalues k and k+2c, proves the decoupling bijection onto two scalar oscillator equations and exhibits the exact solution family with squared frequencies k/m and (k+2c)/m, together with the affine zero-frequency mode in the degenerate case k=0.

The pendulum bound. The pendulum field near the downward equilibrium: bounded force derivative, cubic remainder and the finite-time small-angle error bound verifies the hypotheses of the finite-time linearisation theorem for the field F(θ,v)=(v,−ω2sin⁡θ): the cubic remainder ∣sin⁡θ−θ∣≤∣θ∣3/6, the Lipschitz constant and the operator norm, and then instantiates the order-δ2 error bound on [0,T] under δeLT<r. The scope sentence of that proposition and the page remark Scope and limits of the verification items on this page record that the estimate is conditional, finite-time and makes no exact-harmonicity claim about the nonlinear pendulum.

The companion page cm-math-verification-examples records explicit finite checks of these results: the two-sided saddle values at T=2π, the T=π non-uniqueness with its common action value, the parallel-axis identity on two point masses, the two-mass eigenrelations at concrete parameters, and the pendulum error expression with its exact constants.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicableprecheck passOpen item page →

The fixed-endpoint second variation of a time-dependent action at a curve

Definition

Let n≥1, let a<b be reals, let U⊆Rn be open and nonempty, and let (O,L) be a time-dependent Lagrangian chart over [a,b]×U in the sense of A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations: O⊆R×Rn×Rn is open with

[a,b]×U×Rn⊆O,

and L:O→R is of class C2. Recall from that definition that a C2 curve q:[a,b]→U is the restriction to [a,b] of a C2 map on an open interval, that a variation of q is such a map η:[a,b]→Rn with

η(a)=η(b)=0,

and that the action of a C2 curve is

S(q):=∫abL(t,q(t),q′(t)) dt.

The fixed-endpoint family. Let q:[a,b]→U be a C2 curve and let η be a variation of q. For ε∈R put

qε(t):=q(t)+εη(t)(t∈[a,b]).

By clause 1 of The first variation of a time-dependent action on a finite interval: admissibility, the integral formula and continuity of the Euler-Lagrange coefficient there is ε0>0 such that qε is a C2 admissible curve with values in U for every ε∈(−ε0,ε0), so that the action S(qε) is defined for those ε.

The second variation. The second variation of S at q in the direction η is the second derivative at ε=0 of the function ε↦S(q+εη):

δ2S(q)[η,η]:=d2dε2∣ε=0S(q+εη).

This is a well-defined real number under the present hypotheses: the map ε↦S(q+εη) is twice differentiable at 0, with the explicit integral formula, by The second variation of a time-dependent action: twice differentiability of the fixed-endpoint family and the Hessian integral formula ↗, which depends on this definition only for the notation δ2S(q)[η,η]. The value is a single real number determined by the pair (q,η) and does not depend on the choice of the C2 extensions of q and η, since the values of qε and qε′ on [a,b] are those of any such extension and the derivative on [a,b] is extension-independent (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations).

Scope. The definition introduces the one-parameter family and its second derivative at ε=0, nothing else. It asserts nothing about minimising or maximising S, nothing about the sign of δ2S, and no physical interpretation of L; stationarity of a curve is the separate first-order condition recorded in A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations, and no equivalence with the vanishing of an Euler-Lagrange expression is used at this point.

Remarks

LemmaStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The second variation of a time-dependent action: twice differentiability of the fixed-endpoint family and the Hessian integral formula

Statement

Let (O,L) be a time-dependent Lagrangian chart over [a,b]×U in the sense of A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations (so L is C2 and [a,b]×U×Rn⊆O), let q:[a,b]→U be a C2 curve and let η be a variation of q (a C2 map with η(a)=η(b)=0). Write ψ(ε,t):=(t,q(t)+εη(t),q′(t)+εη′(t)) and, for ∣ε∣<ε0 with ε0 as in clause 1 below,

H(ε,t):=∑j,k=1n[∂2L∂qj∂qk(ψ(ε,t))ηj(t)ηk(t)+2∂2L∂qj∂vk(ψ(ε,t))ηj(t)ηk′(t)+∂2L∂vj∂vk(ψ(ε,t))ηj′(t)ηk′(t)].

Then:

  1. there is ε0>0 such that for every ε∈(−ε0,ε0) the curve q+εη is a C2 admissible curve with values in U, and the map ε↦S(q+εη) is of class C2 on (−ε0,ε0);
  2. for every ε∈(−ε0,ε0),

d2dε2S(q+εη)=∫abH(ε,t) dt;

  1. in particular the second variation of The fixed-endpoint second variation of a time-dependent action at a curve is

δ2S(q)[η,η]=∫abH(0,t) dt,

the formula of clause 2 with every second partial of L evaluated at (t,q(t),q′(t)).

Facts & Assumptions

Given: The chart data (O,L), a C2 curve q:[a,b]→U, a variation η, the map ψ(ε,t)=(t,q(t)+εη(t),q′(t)+εη′(t)), and the action S and second variation of A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations and The fixed-endpoint second variation of a time-dependent action at a curve.

[L1]

Admissibility: there is ε0>0 such that for ∣ε∣<ε0 the curve q+εη is a C2 admissible curve with values in U, hence ψ maps (−ε0,ε0)×[a,b] into O because [a,b]×U×Rn⊆O (The first variation of a time-dependent action on a finite interval: admissibility, the integral formula and continuity of the Euler-Lagrange coefficient, clause 1).

[L2]

Differentiation under the integral sign on a compact rectangle (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral): if g,h:[α,β]×[a,b]→R are continuous and for each fixed t the map ε↦g(ε,t) is differentiable on (α,β) with derivative h(ε,t), then ε↦∫abg(ε,t) dt is differentiable on (α,β) with derivative ∫abh(ε,t) dt.

∂∂ε[L(ψ)]=Lq(ψ)⋅η(t)+Lv(ψ)⋅η′(t),

∂∂ε[Lq(ψ)]j=∑k∂2L∂qj∂qk(ψ)ηk(t)+∑k∂2L∂qj∂vk(ψ)ηk′(t),

∂∂ε[Lv(ψ)]j=∑k∂2L∂vj∂qk(ψ)ηk(t)+∑k∂2L∂vj∂vk(ψ)ηk′(t),

and each displayed expression is continuous on every rectangle [−δ,δ]×[a,b] with 0<δ<ε0, the compositions of the continuous second partials with the continuous map ψ and the continuous functions η,η′ being continuous.

[L4]

Clairaut-Schwarz (Clairaut--Schwarz theorem for continuous second partial derivatives): the second partials of the C2 function L commute on O, so ∂2L/∂vj∂qk=∂2L/∂qk∂vj there.

[L6]

Uniform continuity and uniform integral errors (Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line, Heine-Cantor: a continuous map from a compact metric space to any metric space is uniformly continuous, Integrable functions on [a,b] form a set closed under sums and scalar multiples, and ∫ab(λf+μg)=λ∫abf+μ∫abg, If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a), If f,g are integrable on [a,b] then so are ∣f∣, f2, fg, max⁡(f,g) and min⁡(f,g), and ∣∫abf∣≤∫ab∣f∣): the closed bounded rectangle R=[−δ,δ]×[a,b] is compact, a continuous real function on a compact metric space is uniformly continuous; and if u,v are integrable on [a,b] with ∣u(t)−v(t)∣≤η for every t∈[a,b] and some η≥0, then u−v and ∣u−v∣ are integrable and ∣∫abu−∫abv∣=∣∫ab(u−v)∣≤∫ab∣u−v∣≤η(b−a), by linearity of the integral, the triangle inequality for the integral and the two-sided bound of monotonicity.

Proof

Proof technique: direct; apply differentiation under the integral sign twice on a compact rectangle, the first derivative being the first-variation integrand and the second being the Hessian quadratic form of L along the family.

1.1givenL1L3

Fix 0<δ<ε0 with ε0 from [L1], put R:=[−δ,δ]×[a,b], and define g(ε,t):=L(ψ(ε,t)) on R. By [L1] the values of ψ on R lie in O, so g is defined and continuous on R; by [L3] the ε-partial g1(ε,t)=Lq(ψ)⋅η(t)+Lv(ψ)⋅η′(t) exists on R and is continuous there.

2.1L2L5step 1.1

By [L2] applied to the continuous pair (g,g1) on R, the map F(ε):=S(q+εη)=∫abg(ε,t) dt is differentiable on (−δ,δ) with derivative F′(ε)=∫abg1(ε,t) dt.

2.2L2L3L5step 1.1

By [L3] the ε-partial of g1 exists on R and equals the four-term sum g2(ε,t)=[Lqq(ψ)η⋅η+Lqv(ψ)η′⋅η]+[Lvq(ψ)η⋅η′+Lvv(ψ)η′⋅η′], written with the block matrices Lqq,Lqv,Lvq,Lvv of second partials evaluated at ψ=ψ(ε,t); it is continuous on R. Applying [L2] to the pair (g1,g2) shows that ε↦∫abg1(ε,t) dt is differentiable on (−δ,δ) with derivative ∫abg2(ε,t) dt. Hence F is twice differentiable on (−δ,δ) with F′′(ε)=∫abg2(ε,t) dt, and since every 0<δ<ε0 is allowed, clause 2 holds for every ε∈(−ε0,ε0) once g2 is rewritten in the symmetric form H.

3.1L6L5step 2.2

Continuity of the second derivative: the function g2 is continuous on the compact rectangle R, so [L6] makes it uniformly continuous. Given ε~>0 choose γ>0 such that ∣g2(ε,t)−g2(ε′,t)∣<ε~/(b−a) for every t∈[a,b] whenever ε,ε′∈(−δ,δ) satisfy ∣ε−ε′∣<γ. Then for such ε,ε′, linearity of the integral and [L6] give ∣F′′(ε)−F′′(ε′)∣=∣∫ab(g2(ε,t)−g2(ε′,t)) dt∣≤ε~b−a (b−a)=ε~, so F′′ is continuous on (−δ,δ); since every 0<δ<ε0 is allowed, ε↦S(q+εη) is of class C2 on (−ε0,ε0), which is clause 1.

4.1L4step 2.2given∎

Symmetrisation: by [L4], for all j,k one has (Lvq(ψ))jk=∂2L/∂vj∂qk(ψ)=∂2L/∂qk∂vj(ψ)=(Lqv(ψ))kj, so η′⋅Lvq(ψ)η=η⋅Lqv(ψ)η′ after relabelling the summation indices; therefore g2=H and clauses 2 and 3 follow, the latter by the definition of δ2S(q)[η,η] and evaluation at ε=0.

Remarks

PropositionStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Second variation of the oscillator action on a finite interval: both signs for T>π/ω and the value at T=π/ω

Statement

Let ω>0 and T>0, let U:=R, and consider the time-dependent Lagrangian chart over [0,T]×U with O:=R×R×R and L(q,v):=12(v2−ω2q2) (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations). Write S for its action and δ2S for the second variation of The fixed-endpoint second variation of a time-dependent action at a curve. Then:

  1. the zero curve q0≡0 is a C2 admissible curve and is stationary for S;
  2. for every integer n≥1 the map ηn(t):=sin⁡(nπt/T) is a variation of q0, and

δ2S(q0)[ηn,ηn]=T2((nπT)2−ω2);

moreover S(εηn)=ε22δ2S(q0)[ηn,ηn] for every real ε; 3. if T>π/ω, then δ2S(q0)[η1,η1]<0 and δ2S(q0)[ηn,ηn]>0 for every integer n>ωT/π. Consequently, for every real δ0>0 there exist variations η of q0 with sup⁡t∈[0,T](∣η(t)∣+∣η′(t)∣+∣η′′(t)∣)<δ0 and S(η)<0, and variations with the same bound and S(η)>0; hence the zero curve is neither a local minimiser nor a local maximiser of S, in the sense that every neighbourhood of q0 for the supremum distance sup⁡[0,T](∣q−q0∣+∣q′−q0′∣+∣q′′−q0′′∣) contains curves whose action is below S(q0) and curves whose action is above S(q0); 4. if T=π/ω, then δ2S(q0)[η1,η1]=0, while δ2S(q0)[ηn,ηn]>0 for every integer n≥2.

Facts & Assumptions

Given: The parameters ω>0 and T>0, the chart data (O,L) above, the zero curve q0≡0, the action S, the second variation δ2S, and the maps ηn(t)=sin⁡(nπt/T) for integers n≥1.

[F1]

Chart and stationarity (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations, Stationarity of the time-dependent action is equivalent to the vanishing of the Euler-Lagrange expression): the Euler-Lagrange expression of this chart along a C2 curve q is E(q)(t)=−ω2q(t)−q′′(t), and q is stationary for S if and only if E(q)=0 on [0,T].

[F2]

Second variation formula (The second variation of a time-dependent action: twice differentiability of the fixed-endpoint family and the Hessian integral formula, The fixed-endpoint second variation of a time-dependent action at a curve): for every C2 curve q on [0,T] and every variation η, δ2S(q)[η,η]=∫0T(η′(t)2−ω2η(t)2) dt, because the second partials of L are the constants ∂2L/∂q2=−ω2, ∂2L/∂q∂v=0 and ∂2L/∂v2=1.

[F3]

Trigonometric data (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine defined by their real power series, The addition formulas for sine and cosine, Parity and the Pythagorean identity for sine and cosine, The zero sets of sine and cosine and the least positive common period 2 pi, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)): sin⁡ and cos⁡ are differentiable with sin⁡′=cos⁡ and cos⁡′=−sin⁡, and sin⁡0=0; for every real u the addition formulas and the Pythagorean identity give cos⁡2u=12(1+cos⁡2u) and sin⁡2u=12(1−cos⁡2u); the zero set of sin⁡ is exactly πZ. Hence for every integer n≥1 the map ηn(t)=sin⁡(nπt/T) is C2 on an open interval containing [0,T], with ηn′(t)=(nπ/T)cos⁡(nπt/T) and ηn(0)=ηn(T)=sin⁡(nπ)=0.

[F5]

Archimedean property (Every complete ordered field is Archimedean): for every real x there is a natural number n with x<n.

Proof

Proof technique: direct; compute the second variation of the oscillator family with elementary trigonometric integrals, obtain the exact action along the family, and read off the two signs and the degenerate case. [given]

1.1givenF1F4algebra

Stationarity of the zero curve: q0 is a C2 admissible curve with values in U=R, and E(q0)(t)=−ω2⋅0−0=0 for every t, so q0 is stationary by [F1]; also S(q0)=∫0T0 dt=0 by [F4].

1.2givenF2F3

Variation and second variation: for every integer n≥1, ηn is a variation of q0 by [F3], so [F2] and [F3] give δ2S(q0)[ηn,ηn]=∫0T((nπT)2cos⁡2(nπtT)−ω2sin⁡2(nπtT))dt.

2.1step 1.2F3F4algebra

Evaluation of the integrals: by [F3] and [F4], ∫0Tcos⁡2(nπtT)dt=12∫0T(1+cos⁡2nπtT)dt=T+02=T2, and likewise ∫0Tsin⁡2(nπt/T) dt=T/2, because the FTC gives ∫0Tcos⁡(2nπt/T) dt=T2nπ[sin⁡(2nπt/T)]0T=T2nπ(sin⁡(2nπ)−sin⁡0)=0 by the vanishing of sin⁡ at integer multiples of π from [F3]. Substituting into step 1.2 gives the second-variation value of clause 2.

3.1step 1.2step 2.1F4algebra

Exact action along the family: since L is homogeneous of degree two in (q,v), for every real ε the integrand of S(εηn) is L(εηn,εηn′)=ε2L(ηn,ηn′); by linearity of the integral [F4] and step 1.2, S(εηn)=ε2S(ηn)=ε22δ2S(q0)[ηn,ηn], which is the exact identity of clause 2.

4.1step 1.1step 2.1step 3.1F3F5algebra

Signs for T>π/ω: if T>π/ω, then π/T<ω, so clause 2 gives δ2S(q0)[η1,η1]<0; by [F5] choose an integer n with n>ωT/π, so (nπ/T)2>ω2 and δ2S(q0)[ηn,ηn]>0. For any real δ0>0 pick ε>0 with ε(1+nπ/T+(nπ/T)2)<δ0; then ∣ηn∣≤1, ∣ηn′∣≤nπ/T and ∣ηn′′∣≤(nπ/T)2 by [F3], so the variations εη1 and εηn satisfy the displayed supremum bound, while step 3.1 gives S(εη1)=ε22δ2S(q0)[η1,η1]<0 and S(εηn)>0; with S(q0)=0 from step 1.1, the zero curve is neither a local minimiser nor a local maximiser.

5.1step 2.1algebra∎

The degenerate case: if T=π/ω, substitution into clause 2 gives δ2S(q0)[η1,η1]=T2(ω2−ω2)=0, while for every integer n≥2 the value is T2(n2ω2−ω2)>0.

Remarks

PropositionStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The undamped oscillator: Dirichlet non-uniqueness at T=nπ/ω despite unique initial-value solutions

Statement

Let m>0 and k>0, put ω:=(k/m)1/2>0 (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}), and let T>0. Consider the undamped homogeneous oscillator equation

m x′′(t)+k x(t)=0

of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with time domain I=R and zero forcing, and its solutions on intervals with nonempty interior. The Dirichlet curves below are restrictions to [0,T] of its global solutions. Then:

  1. for every real A the function xA(t):=Asin⁡(ωt) is a C2 (indeed C∞) solution on R with xA(0)=0 and xA(T)=Asin⁡(ωT); moreover xA(T)=0 if and only if A=0 or ωT∈πZ (The zero sets of sine and cosine and the least positive common period 2 pi);
  2. if T=nπ/ω for some integer n≥1, then the set of C2 solutions of the Dirichlet problem on [0,T],

m x′′+k x=0,x(0)=x(T)=0,

is exactly {xA:A∈R}; this set is a one-dimensional real vector space, so the boundary-value problem has infinitely many distinct solutions with the same endpoint data (0,0), one for each initial velocity Aω; 3. by contrast, for every (x0,v0)∈R2 the initial-value problem m x′′+k x=0, x(0)=x0, x′(0)=v0, has exactly one solution on R, namely x(t)=x0cos⁡(ωt)+(v0/ω)sin⁡(ωt) (The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data, clauses 3 and 4 with c=0); uniqueness for the initial-value problem therefore does not imply uniqueness for the boundary-value problem; 4. for the chart with O=R×R×R and L(q,v)=12(v2−ω2q2) over [0,T]×R (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations), each xA is a stationary curve of the action; and if T=nπ/ω, the stationary curves on [0,T] satisfying the endpoint conditions x(0)=x(T)=0 are exactly the curves xA, so the fixed-endpoint problem with zero data has infinitely many stationary curves (Stationarity of the time-dependent action is equivalent to the vanishing of the Euler-Lagrange expression).

Facts & Assumptions

Given: The parameters m>0, k>0, T>0, the number ω=(k/m)1/2, the oscillator model of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, the curves xA(t)=Asin⁡(ωt), and the chart of clause 4.

[F1]

Classification and uniqueness (The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data, The damped driven oscillator model: second-order equation, phase-plane system and characteristic data): for the parameters (m,c,k)=(m,0,k) one has Δ=−4mk<0, γ=0 and ωcl:=4mk−c2/(2m)=k/m=ω, so clause 3 states that x:J→R on an interval J with nonempty interior is a solution if and only if x(t)=Acos⁡(ωt)+Bsin⁡(ωt) for some A,B∈R; and clause 4 states that for every (x0,v0)∈R2 there is exactly one solution x:R→R with x(0)=x0, x′(0)=v0.

[F2]

Trigonometric data (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine defined by their real power series, The zero sets of sine and cosine and the least positive common period 2 pi, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set): sin⁡ and cos⁡ are differentiable with sin⁡′=cos⁡, cos⁡′=−sin⁡ and sin⁡0=0; the zero set of sin⁡ is exactly πZ; the chain rule gives ddtsin⁡(ωt)=ωcos⁡(ωt) and d2dt2sin⁡(ωt)=−ω2sin⁡(ωt).

[F3]

Chart and stationarity (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations, Stationarity of the time-dependent action is equivalent to the vanishing of the Euler-Lagrange expression): the Euler-Lagrange expression along a C2 curve q is E(q)=−ω2q−q′′, and q is stationary for the action S of this chart if and only if E(q)=0 on [0,T].

Proof

Proof technique: direct; verify the explicit family, determine its endpoint values from the zero set of the sine, classify the Dirichlet solutions through the oscillator classification theorem, and identify the stationary curves through the Euler-Lagrange expression. [given]

1.1givenF2F4algebra

The family xA(t)=Asin⁡(ωt): by [F2] the curve is C∞ with xA′′(t)=−Aω2sin⁡(ωt), so mxA′′(t)+kxA(t)=(−mω2+k)Asin⁡(ωt)=0 by ω2=k/m from [F4]; also xA(0)=0 and xA(T)=Asin⁡(ωT).

2.1step 1.1F2algebra

Endpoint vanishing: by the zero-set statement of [F2], sin⁡(ωT)=0 is equivalent to ωT∈πZ, so xA(T)=0 holds exactly when A=0 or ωT∈πZ.

2.2step 1.1F1F2algebra

Classification of the Dirichlet solutions: assume T=nπ/ω with an integer n≥1 and let x be a C2 solution of mx′′+kx=0 on [0,T] with x(0)=x(T)=0. By clause 3 of [F1] there are A,B∈R with x(t)=Acos⁡(ωt)+Bsin⁡(ωt) for t∈[0,T]; the condition x(0)=0 gives A=0, and then x(T)=Bsin⁡(nπ)=0 holds automatically by [F2], so x=xB. Conversely every xB is a solution by step 1.1 and satisfies xB(T)=Bsin⁡(nπ)=0. Hence the solution set of the Dirichlet problem is exactly {Bsin⁡(ωt):B∈R}={xA:A∈R}, which is a linear subspace containing the nonzero function sin⁡(ωt) and is therefore one-dimensional.

3.1step 2.2F1F2algebra

Initial-value uniqueness and contrast: by clause 4 of [F1] the initial-value problem with data (x0,v0) has exactly one solution on R, with the displayed trigonometric form; for A≠A′ the solutions xA and xA′ of the boundary-value problem are distinct (they differ at t=π/(2ω), which lies in [0,T] because T≥π/ω, and sin⁡(π/2)≠0 by the zero-set statement of [F2]) although they share the endpoint data (0,0) and differ in initial velocity Aω≠A′ω; this is the announced contrast between the two problems.

4.1step 1.1step 2.2F3∎

Stationary curves: by [F3] and step 1.1, E(xA)=−ω2xA−xA′′=0 on [0,T], so every xA is stationary for the action; if T=nπ/ω, a stationary curve on [0,T] with x(0)=x(T)=0 is a C2 solution of mx′′+kx=0 with those endpoint values by [F3], hence is one of the xA by step 2.2. This proves clause 4.

Remarks

  • At T=π/ω the second-variation value in direction η1=sin⁡(ωt) of Second variation of the oscillator action on a finite interval: both signs for T>π/ω and the value at T=π/ω is zero and the Dirichlet solution space is one-dimensional. More generally, its nonzero solutions occur exactly when T=nπ/ω for an integer n≥1, as clauses 1-2 show. The computations here are mathematical statements about the model; no physical forcing or measured period is asserted.

  • The one-parameter family {xA} is not an artefact of dropping a hypothesis: it persists as the full solution set of the Dirichlet problem, and it is a legitimate underdetermined boundary-value problem, in contrast with the exactly solvable initial-value problem.

DefinitionDefinition: AI-adaptedProof: Not applicableprecheck passOpen item page →

The inertia tensor of a finite bounded Borel measure about any point of R3

Definition

Let B(R3) be the Borel sigma-algebra of R3 (The Borel sigma-algebra of a topological space), let ν be a measure on the measurable space (R3,B(R3)) (Measures on sigma-algebras, Measure spaces) that is finite and bounded: ν(R3)<+∞ and ν(R3∖B‾R(0))=0 for some real R>0 (Open ball, closed ball and sphere in a metric space). Put m:=ν(R3)∈[0,+∞); nonnegativity follows from the definition of measure Measures on sigma-algebras. If m>0, then ν is a mass distribution in the sense of Mass distributions on R3: total mass, first moment and inertia matrix. If m=0, then for every Borel set E⊆R3 one has 0≤ν(E)≤ν(R3)=0, so ν is the zero measure. Integrals against ν are Lebesgue integrals as in Integral over a measurable subset and Integrable real and complex functions, and their integrals, defined componentwise for vector- and matrix-valued integrands. When m>0, recall that the inertia matrix of ν is

I(ν)=∫(∥ρ∥22I3−ρρT) dν(ρ)∈M3(R),

with I3 the identity matrix and entrywise integrals I(ν)jk=∫(∥ρ∥22δjk−ρjρk) dν, where ∥⋅∥2 and the matrix conventions are those of The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations and Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes.

Shifted integrand and inertia tensor about a point. For a∈R3 define the component functions

Ga(ρ)jk:=∥ρ−a∥22 δjk−(ρj−aj)(ρk−ak)(j,k∈{1,2,3}),

equivalently the matrix-valued function Ga(ρ)=∥ρ−a∥22I3−(ρ−a)(ρ−a)T, and define the inertia tensor of ν about a Iν(a)∈M3(R) by its entries

Iν(a)jk:=∫Ga(ρ)jk dν(ρ).

Under the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)) each entry is a finite real number; this is the content of Finiteness and a uniform bound for the entries of the inertia tensor about an arbitrary point ↗, which depends on this definition only for its notation. For a=0 the shifted integrand is the parent integrand ∥ρ∥22I3−ρρT. Thus Iν(0)=I(ν) when m>0; if m=0, the parent definition does not define I(ν), while the shifted integral remains defined and equals the zero matrix 03×3 for every a.

Scope. The definition names a mathematical function a↦Iν(a) from R3 to M3(R) determined by any finite bounded Borel measure ν and its coordinate origin. For the zero measure this function is identically 03×3. No physical body, unit system or measured mass is asserted. No translation identity for nonzero a is asserted here; the centred, positive-total-mass case is proved separately on this page.

Remarks

  • The measure need not be atomless. Points, curves, surfaces and volumes are all admissible: ν is an arbitrary bounded finite Borel measure, so finite sums of Dirac measures, lines and surface measures embedded in R3 and absolutely continuous volume measures are all special cases, exactly as in Mass distributions on R3: total mass, first moment and inertia matrix.

  • Translation of the point, not of the measure. The point a only shifts the reference point at which the second moments are computed; the measure and its coordinate representation are unchanged, and μ(ν) is not assumed to vanish unless a consumer states it.

LemmaStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Finiteness and a uniform bound for the entries of the inertia tensor about an arbitrary point

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)). Let ν be a measure on (R3,B(R3)) that is finite and bounded: ν(R3)<+∞ and ν(R3∖B‾R(0))=0 for some real R>0 (Open ball, closed ball and sphere in a metric space), and let a∈R3. With the shifted integrands Ga(ρ)jk=∥ρ−a∥22δjk−(ρj−aj)(ρk−ak) of The inertia tensor of a finite bounded Borel measure about any point of R3 (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn for the norm), then:

  1. for all j,k the function Ga(⋅)jk:R3→R is Borel measurable and ν-integrable (Borel measurable and Lebesgue measurable functions on Rn);
  2. each entry Iν(a)jk is a finite real number and

∣Iν(a)jk∣≤2(R+∥a∥2)2 ν(R3).

Facts & Assumptions

Given: A finite bounded Borel measure ν on R3 with ν(R3∖B‾R(0))=0, and a point a∈R3.

[L1]

Under ACω every continuous map R3→R is Borel measurable (Continuous functions on Euclidean spaces are Borel measurable, Borel measurable and Lebesgue measurable functions on Rn), and sums, scalar multiples and pointwise products of measurable real-valued functions are measurable (Arithmetic and lattice operations preserve measurability whenever they are defined).

[L2]

Estimates on the ball: for ρ∈B‾R(0) one has ∥ρ∥2≤R and therefore ∥ρ−a∥2≤∥ρ∥2+∥a∥2≤R+∥a∥2 by the triangle inequality for the Euclidean norm (The inner-product norm is definite, homogeneous, and satisfies the triangle inequality, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn); and for each index j one has (ρj−aj)2≤∑l(ρl−al)2=∥ρ−a∥22, because each term of a finite sum of nonnegative reals is at most the sum (Laws of finite sums and finite products, clause 4, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L3]

The elementary inequality 2uv≤u2+v2 for reals u,v follows from (u−v)2≥0 by squaring rules (Squares of nonzero elements are positive, Absolute value in an ordered field).

[L4]

Nonnegative integrals: the constant function with value c≥0 has ∫c dν=c ν(R3) via the simple-function integral; the nonnegative integral is monotone and homogeneous for nonnegative scalars (The nonnegative Lebesgue integral, The integral of a nonnegative simple function, Monotonicity and nonnegative homogeneity of the nonnegative integral).

[L5]

A real measurable function f with ∫∣f∣ dν<+∞ is ν-integrable, and then ∣∫f dν∣≤∫∣f∣ dν (Integrable real and complex functions, and their integrals, The modulus of an integral is bounded by the integral of the modulus, Integral over a measurable subset).

[L6]

Indefinite integrals and null sets (The indefinite integral of a nonnegative measurable function is a measure, A nonnegative integral over a null set vanishes, Integral over a measurable subset): for every nonnegative measurable f the set function A↦∫Af dν is a measure, so ∫A∪Bf dν=∫Af dν+∫Bf dν for disjoint measurable A,B; and if A is ν-null then ∫Af dν=0.

Proof

Proof technique: direct; prove measurability, bound the shifted integrand pointwise on the support ball by a constant, and integrate the constant.

1.1givenL1L2

Measurability: the coordinate maps ρ↦ρj, the identity maps ρ↦aj, and ρ↦∥ρ−a∥22=∑l(ρl−al)2 are continuous, hence Borel measurable by [L1]; sums and products of these measurable functions are measurable by [L1], so each Ga(⋅)jk=∥ρ−a∥22δjk−(ρj−aj)(ρk−ak) is Borel measurable.

1.2givenL2L3algebra

Pointwise bound on the support ball: for ρ∈B‾R(0), by [L2] and [L3], ∣Ga(ρ)jk∣≤∥ρ−a∥22+∣ρj−aj∣ ∣ρk−ak∣≤∥ρ−a∥22+12((ρj−aj)2+(ρk−ak)2)≤2∥ρ−a∥22≤2(R+∥a∥2)2.

2.1step 1.1step 1.2givenL4L5L6algebra∎

Integrability and the bound: the function ∣Ga(⋅)jk∣ is nonnegative and measurable by step 1.1, so [L6] applied to the disjoint measurable sets B‾R(0) and N:=R3∖B‾R(0), the latter ν-null, gives ∫∣Ga(⋅)jk∣ dν=∫B‾R(0)∣Ga(⋅)jk∣ dν+∫N∣Ga(⋅)jk∣ dν=∫B‾R(0)∣Ga(⋅)jk∣ dν. On B‾R(0) the integrand is bounded by the constant M:=2(R+∥a∥2)2 by step 1.2, so monotonicity in [L4] applied to the pointwise inequality ∣Ga(⋅)jk∣χB‾R(0)≤MχB‾R(0)≤M, and the constant-value formula, give ∫∣Ga(⋅)jk∣ dν≤∫M dν=M ν(R3)<+∞. By step 1.1 and [L5] the function Ga(⋅)jk is ν-integrable, its integral Iν(a)jk is a finite real number, and ∣Iν(a)jk∣=∣∫Ga(⋅)jk dν∣≤∫∣Ga(⋅)jk∣ dν≤2(R+∥a∥2)2ν(R3).

Remarks

  • The bound is uniform in the shape of the measure: it uses only the support radius R, the shift ∥a∥2 and the total mass ν(R3), not the detailed distribution. For a=0 it reproduces the bound 2R2ν(R3) of the shifted integrand computed at the origin.

  • The Axiom of Countable Choice is used exactly once, in [L1], to obtain the Borel measurability of the continuous integrands; every later consumer of this finiteness statement carries the assumption explicitly.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Parallel-axis theorem: the inertia tensor about a displaced point of a centred mass distribution

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), used only to make the displayed integrals finite. Let ν be a mass distribution on R3 with total mass m=ν(R3)∈(0,+∞) that is centred, that is, μ(ν)=∫ρ dν(ρ)=0∈R3 (Mass distributions on R3: total mass, first moment and inertia matrix), and let I(ν) be its inertia matrix. Then for every a∈R3 the inertia tensor of ν about a (The inertia tensor of a finite bounded Borel measure about any point of R3) satisfies

Iν(a)=I(ν)+m(∥a∥22I3−aaT),

where ∥⋅∥2 is the Euclidean norm, I3 the identity matrix and aaT the matrix product of The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations and Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes. Entrywise,

Iν(a)jk=I(ν)jk+m(∥a∥22δjk−ajak)(j,k∈{1,2,3}).

Facts & Assumptions

Given: A centred mass distribution ν with total mass m, a point a∈R3, and the notation of Mass distributions on R3: total mass, first moment and inertia matrix and The inertia tensor of a finite bounded Borel measure about any point of R3: the first moment μ(ν), the inertia matrix I(ν) with I(ν)jk=∫(∥ρ∥22δjk−ρjρk) dν, and the inertia tensor Iν(a) with Iν(a)jk=∫(∥ρ−a∥22δjk−(ρj−aj)(ρk−ak)) dν.

[L1]

Finiteness and integrability (Finiteness and a uniform bound for the entries of the inertia tensor about an arbitrary point, Finiteness of the moments of a bounded finite Borel measure on R3, Mass distributions on R3: total mass, first moment and inertia matrix): the functions ρ↦ρj, ρ↦ρjρk, and the shifted polynomials ∥ρ−a∥22δjk−(ρj−aj)(ρk−ak) are ν-integrable; the entries Iν(a)jk, I(ν)jk and μ(ν)j are finite real numbers.

[L2]

Linearity of the Lebesgue integral on L1 (The Lebesgue integral is linear on L1(μ), Integrable real and complex functions, and their integrals, Integral over a measurable subset): for real α,β and integrable real functions f,g, the combination αf+βg is integrable and ∫(αf+βg) dν=α∫f dν+β∫g dν.

[L3]

Centredness and total mass (Mass distributions on R3: total mass, first moment and inertia matrix, The nonnegative Lebesgue integral, The integral of a nonnegative simple function): ∫ρj dν=μ(ν)j=0 for each j, and the constant function 1 has ∫1 dν=ν(R3)=m.

[L4]

Component algebra (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes): for all ρ,a∈R3 and all j,k, ∥ρ−a∥22=∥ρ∥22−2⟨ρ,a⟩+∥a∥22,(ρj−aj)(ρk−ak)=ρjρk−ρjak−ajρk+ajak, and ⟨ρ,a⟩=∑lalρl.

Proof

Proof technique: direct; expand the shifted integrand into the parent integrand plus constants and linear terms, and integrate termwise, the linear terms vanishing by centredness.

1.1givenL4algebra

Expansion of the integrand: for every ρ∈R3 and all j,k, [L4] gives ∥ρ−a∥22δjk−(ρj−aj)(ρk−ak)=[∥ρ∥22δjk−ρjρk]+[∥a∥22δjk−ajak]+[−2⟨ρ,a⟩δjk+ρjak+ajρk].

2.1step 1.1givenL1L2

Integrability of the pieces: the first bracket is integrable by [L1]; the second bracket is a constant function and hence integrable because ν(R3)<+∞; the functions ρ↦ρj are integrable by [L1], so the third bracket is a finite real linear combination of integrable functions and is integrable by [L2].

3.1step 1.1step 2.1L1L2L3algebra

Termwise integration: applying [L2] to the decomposition of step 1.1 and using the defining integrals of [L1] and the values of [L3], Iν(a)jk=I(ν)jk+m(∥a∥22δjk−ajak)−2δjk∑lalμ(ν)l+akμ(ν)j+ajμ(ν)k, and every term of the last group vanishes because μ(ν)=0.

4.1step 3.1L4algebra∎

Conclusion in matrix form: reading the entrywise identity of step 3.1 as a matrix identity gives Iν(a)=I(ν)+m(∥a∥22I3−aaT), which is the statement.

Remarks

  • Centredness is exactly what is used. The proof shows that for a general mass distribution the same expansion leaves the additional terms −2(⟨a,μ(ν)⟩)δjk+akμ(ν)j+ajμ(ν)k; the centred case is the classical parallel-axis form, and the theorem asserts nothing about the non-centred case beyond this displayed expansion.

  • Shape independence. The derivation uses only the finiteness of the first and second moments and the vanishing of the first moment, so the conclusion covers positive-total-mass point, line, surface and volume distributions alike: ν is an arbitrary bounded finite Borel measure of positive total mass, as in Mass distributions on R3: total mass, first moment and inertia matrix.

PropositionStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Two equal masses with two outer springs and a coupling spring: eigenmodes, decoupling and the explicit normal frequencies

Statement

Let m>0, k≥0 and c≥0, and put

K:=(k+c−c−ck+c),u+:=12(11),u−:=12(1−1),

and ω+:=(k/m)1/2, ω−:=((k+2c)/m)1/2 (nonnegative square roots, Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}). Consider the linear system

m x′′(t)=−Kx(t)

for C2 curves x:J→R2 on intervals J⊆R with nonempty interior (Intervals of R: the nine order-convex forms, nondegeneracy, and length, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral). The vector equation and its scalar mode equations are autonomous equations with time domain R; a solution on J is their restriction to that solution interval. Every such local solution extends uniquely to the displayed global family. Global growth claims below refer to that extension. Then:

  1. Ku+=ku+ and Ku−=(k+2c)u−; the vectors u+,u− are orthonormal, and every x∈R2 satisfies x=⟨u+,x⟩u++⟨u−,x⟩u− (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn);
  2. x is a solution if and only if there are solutions y+,y−:J→R of the scalar equations

m y+′′+k y+=0,m y−′′+(k+2c) y−=0

of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with x=y+u++y−u−; the pair (y+,y−) is unique and is given by y±=⟨u±,x⟩; 3. consequently, for k>0 the solutions are exactly the curves

x(t)=(A+cos⁡(ω+t)+B+sin⁡(ω+t))u++(A−cos⁡(ω−t)+B−sin⁡(ω−t))u−(A±,B±∈R),

that is, the system is a superposition of an in-phase normal mode with squared frequency k/m and an anti-phase normal mode with squared frequency (k+2c)/m (The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data); 4. if k=0 and c≥0, then Ku+=0, and the curves x(t)=(A+Bt)u+ (A,B∈R) are solutions; for B≠0 each such curve has a global extension unbounded on [0,∞), so the zero-frequency mode is affine rather than oscillatory. The bijection of clause 2 remains valid with the scalar equations m y+′′=0 and m y−′′+2c y−=0, and for c>0 the u−-component is the trigonometric family with ω−=(2c/m)1/2 (The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data); if c=0, it is affine as well (The degenerate and repulsive oscillator: solutions for k=0 and the explicit growing solution for k<0).

Facts & Assumptions

Given: The parameters m>0, k≥0, c≥0, the matrix K and vectors u±, and the system mx′′=−Kx on intervals J with nonempty interior.

[F1]

Matrix and vector data (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations, The transpose AT of a matrix, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn): K is symmetric, K12=K21=−c; the Euclidean inner product of R2 is the componentwise sum ⟨u,v⟩=u1v1+u2v2; and ⟨u+,u+⟩=⟨u−,u−⟩=1, ⟨u+,u−⟩=0.

[F2]

Scalar equations and their solutions (The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data, The degenerate and repulsive oscillator: solutions for k=0 and the explicit growing solution for k<0, Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}, The addition formulas for sine and cosine): for μ>0, put ω=(μ/m)1/2. For any scalar solution y:J→R of my′′+μy=0, choose t0∈J, put J0:={s∈R:t0+s∈J} and define y~(s):=y(t0+s) for s∈J0. Then 0∈J0, J0 has nonempty interior, and autonomy gives my~′′+μy~=0 on J0. The homogeneous-oscillator classification applies there, so y~(s)=A~cos⁡(ωs)+B~sin⁡(ωs). The addition formulas then give y(t)=Acos⁡(ωt)+Bsin⁡(ωt) on J, where A=A~cos⁡(ωt0)−B~sin⁡(ωt0) and B=A~sin⁡(ωt0)+B~cos⁡(ωt0). Conversely these unshifted families solve the scalar equation on every interval. For the zero-stiffness equation my′′=0, the same translation and the degenerate oscillator classification give y~(s)=A~+B~s; restoring the constants gives y(t)=A+Bt with A=A~−B~t0 and B=B~, and every such affine function solves the equation.

[F3]

Vector-valued calculus (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral, Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set, Intervals of R: the nine order-convex forms, nondegeneracy, and length): derivatives of R2-valued curves are computed componentwise; for a C2 curve x and a constant u∈R2 the real function t↦⟨u,x(t)⟩ is C2 with first derivative ⟨u,x′(t)⟩ and second derivative ⟨u,x′′(t)⟩; and if y+,y− are C2 real functions then t↦y+(t)u++y−(t)u− is a C2 curve with second derivative y+′′u++y−′′u−.

Proof

Proof technique: direct; verify the eigenrelations by matrix multiplication, decouple the system in the orthonormal eigenbasis, and read the scalar solution families off the classification theorem.

1.1givenF1algebra

Eigenrelations and expansion: direct multiplication gives Ku+=12((k+c)−c−c+(k+c))=k2(11)=ku+,Ku−=12((k+c)+c−c−(k+c))=(k+2c)u−, and ⟨u+,u−⟩=12(1−1)=0, ⟨u±,u±⟩=12(1+1)=1; also ⟨u+,x⟩=12(x1+x2) and ⟨u−,x⟩=12(x1−x2), so ⟨u+,x⟩u++⟨u−,x⟩u−=x for every x, and K is symmetric since K12=K21=−c.

2.1givenstep 1.1F3algebra

From a solution to scalar solutions: let x:J→R2 be a C2 solution of mx′′=−Kx and set y±(t):=⟨u±,x(t)⟩. By [F3] the functions y± are C2 and my±′′=m⟨u±,x′′⟩=−⟨u±,Kx⟩; since K is symmetric, ⟨u±,Kx⟩=u±TKx=(Ku±)Tx=⟨Ku±,x⟩, and step 1.1 gives Ku+=ku+, Ku−=(k+2c)u−; hence my+′′+ky+=0 and my−′′+(k+2c)y−=0 on J.

2.2givenstep 1.1F3algebra

From scalar solutions to a solution: conversely, let y+,y−:J→R be C2 solutions of the two scalar equations and put x:=y+u++y−u−. By [F3] the curve x is C2 and mx′′=my+′′u++my−′′u−=−ky+u+−(k+2c)y−u−=−(y+Ku++y−Ku−)=−Kx by step 1.1, so x solves the system.

3.1step 1.1step 2.1step 2.2

Bijection: the two constructions are mutually inverse. If x is a solution, then the pair y±=⟨u±,x⟩ of step 2.1 reconstructs x by the expansion identity of step 1.1; if (y+,y−) is a pair of scalar solutions, then the coordinates of x=y+u++y−u− are ⟨u±,x⟩=⟨u±,u+⟩y++⟨u±,u−⟩y−=y± by the orthonormality of step 1.1. Hence clause 2 holds, including uniqueness of the pair.

4.1step 3.1F2algebra

Frequencies on an arbitrary interval: assume k>0 and choose t0∈J. For scalar components y± from clause 2, the translated functions y~±(s):=y±(t0+s) are defined on J−t0, which contains 0 and has nonempty interior, and solve the same autonomous scalar equations. Apply [F2] to obtain the trigonometric families in s with frequencies ω+=(k/m)1/2 and ω−=((k+2c)/m)1/2. The addition formulas in [F2] restore constants A±,B± and give y±(t)=A±cos⁡(ω±t)+B±sin⁡(ω±t) on J. Conversely each displayed scalar family solves its equation on J; the bijection of step 3.1 therefore gives exactly the vector solutions in clause 3, with the u+-mode (both coordinates equal) in phase and the u−-mode (opposite coordinates) anti-phase.

5.1step 1.1step 3.1F2algebra∎

Degenerate case on an arbitrary interval: let k=0, c≥0, and choose t0∈J. Step 1.1 gives Ku+=0, so for x(t)=(A+Bt)u+ one has mx′′=0 and −Kx=−(A+Bt)Ku+=0: the family consists of solutions. The same formula defines x on all of R, and if B≠0 its global extension has ∥x(t)∥2=∣A+Bt∣ unbounded on [0,∞). The translations y~±(s):=y±(t0+s) of the scalar components are defined on J−t0 and satisfy my~+′′=0, my~−′′+2cy~−=0. By [F2], y~+ is affine, and restoring constants gives y+(t)=A+Bt; if c=0 the same affine classification applies to y~−, while if c>0 the homogeneous-oscillator classification gives the trigonometric family for y~− and the addition formulas restore its constants to the stated family with ω−=(2c/m)1/2. The bijection of step 3.1 remains valid and establishes clause 4. In both stiffness cases the scalar classifications give a unique global extension of each component, and the inverse modal transform gives the unique global vector extension claimed in the statement.

Remarks

  • The system is the Euler-Lagrange equation of the Lagrangian 12m(x˙12+x˙22)−12k(x12+x22)−12c(x1−x2)2 for two equal masses attached to fixed walls by equal springs of stiffness k and coupled by a middle spring of stiffness c; the proposition states the mathematics of that model only, and no physical identification of m,k,c is asserted.

  • The k=0 clause is the zero-frequency caution: a vanishing eigenvalue produces affine rather than oscillatory motion, so no nonzero frequency can be read off from the eigenrelation Ku+=ku+ in that case.

PropositionStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The pendulum field near the downward equilibrium: bounded force derivative, cubic remainder and the finite-time small-angle error bound

Statement

Let g>0 and l>0, put ω:=(g/l)1/2>0 (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}), and define the field

F:R2→R2,F(θ,v):=(v,−ω2sin⁡θ),

so that z′=F(z) for z=(θ,v) is the first-order form of the pendulum equation θ′′=−ω2sin⁡θ (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral for the componentwise first-order reduction). Put

A:=(01−ω20),L:=max⁡{1,ω2}.

Then:

  1. F is of class C∞ (in particular continuous and locally Lipschitz in the state variable), F(0)=0, A=DF(0) is the total derivative of F at the origin, and for every r>0 and every z=(θ,v) with ∥z∥2≤r (Open ball, closed ball and sphere in a metric space for the ball, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn for the norm)

∥F(z)−Az∥2=ω2∣sin⁡θ−θ∣≤ω26∣θ∣3≤ω2r6∥z∥22,∥F(z)∥2≤L∥z∥2;

moreover F is globally Lipschitz in the state variable with constant L (Lipschitz map, α-Hölder map for rational 0<α≤1, and contraction), and the operator norm of A for the Euclidean norm is ∥A∥=L (The operator norm as the least bound and as the unit-sphere or unit-ball supremum); 2. let z0∈R2 with δ:=∥z0∥2>0, let T>0 and r>0 satisfy δeLT<r, and let y:R→R2 be the solution of the linear initial value problem y′=Ay, y(0)=z0 (Variation of constants for z′=Az+b(t) via the matrix exponential). Then the maximal solution z of z′=F(z), z(0)=z0, is defined on [0,T] with ∥z(t)∥2≤δeLt for t∈[0,T], and for every t∈[0,T]

∥z(t)−y(t)∥2≤ω2r6δ2∫0te∥A∥(t−s)e2Ls ds≤ω2r6δ2 t e(∥A∥+2L)t

(Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate); 3. the linear system y′=Ay is the linearised pendulum: writing y=(y1,y2) one has y1′′=−ω2y1 and y2=y1′, and the first components satisfy ∣z1(t)−y1(t)∣≤∥z(t)−y(t)∥2 with the bounds of clause 2. The estimate is a finite-time statement of order δ2: it neither states nor implies that the nonlinear pendulum is exactly harmonic, and it makes no claim beyond [0,T].

Facts & Assumptions

Given: The parameters g>0, l>0, the number ω=(g/l)1/2, the field F(θ,v)=(v,−ω2sin⁡θ), the matrix A and the constant L=max⁡{1,ω2}.

[F1]

The linearisation theorem (Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate): if F:U→Rn is continuous and locally Lipschitz in the state variable on an open U∋0, and r>0, C≥0, L≥0, A∈Mn(R) satisfy B‾(0,r)⊆U, ∥F(z)−Az∥2≤C∥z∥22, ∥F(z)∥2≤L∥z∥2 for ∥z∥2≤r, then for every datum with δ:=∥z0∥2>0 and every T>0 with δeLT<r the maximal solution z of z′=F(z), z(0)=z0 is defined on [0,T], satisfies ∥z(t)∥2≤δeLt, and with the solution y of y′=Ay, y(0)=z0 obeys the two displayed bounds with ∥A∥ the operator norm.

[F2]

Trigonometric calculus (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine defined by their real power series, Signs, monotonicity intervals, and ranges of sine and cosine, The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a), A uniform derivative bound gives a uniform Taylor remainder bound, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)): sin⁡ and cos⁡ are differentiable with sin⁡′=cos⁡, cos⁡′=−sin⁡, sin⁡0=0; both have range [−1,1]; for all reals θ,θ′ the mean value theorem gives ∣sin⁡θ−sin⁡θ′∣≤∣θ−θ′∣ and ∣sin⁡θ∣≤∣θ∣; and the Taylor remainder bound with n=2 and M=1 gives ∣sin⁡θ−θ∣≤∣θ∣3/6 for every real θ.

[F4]

Norms (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Real and complex inner product spaces, with the inner product linear in the first argument, The operator norm as the least bound and as the unit-sphere or unit-ball supremum, Open ball, closed ball and sphere in a metric space): for z=(θ,v), ∥z∥22=θ2+v2, ∥Az∥22=ω4θ2+v2, ∥F(z)∥22=v2+ω4sin⁡2θ and ∥F(z)−Az∥22=ω4(sin⁡θ−θ)2; the operator norm of A is the supremum of ∥Az∥2 over ∥z∥2≤1.

Proof

Proof technique: direct; compute the derivative of the pendulum field, bound the quadratic remainder by the cubic sine estimate, identify the Lipschitz and operator norms, and instantiate the linearisation theorem.

1.1givenF2F3algebra

Derivative and regularity: by [F3], F is C1 with DF(θ,v)=(01−ω2cos⁡θ0); since sin⁡ and cos⁡ are C∞ the field is in fact C∞. Also F(0,0)=(0,−ω2sin⁡0)=(0,0) and A=DF(0,0)=(01−ω20).

1.2givenF2F4algebra

Remainder and growth bounds: for ∥z∥2≤r, [F2] gives ∥F(z)−Az∥2=ω2∣sin⁡θ−θ∣≤ω26∣θ∣3, and ∣θ∣3≤rθ2≤r(θ2+v2)=r∥z∥22; moreover ∥F(z)∥22=v2+ω4sin⁡2θ≤v2+ω4θ2≤L2(θ2+v2) because L=max⁡{1,ω2}, so ∥F(z)∥2≤L∥z∥2.

2.1step 1.2F2F4algebra

Lipschitz constant and operator norm: for z=(θ,v) and z′=(θ′,v′), ∥F(z)−F(z′)∥22=(v−v′)2+ω4(sin⁡θ−sin⁡θ′)2≤(v−v′)2+ω4(θ−θ′)2≤L2∥z−z′∥22 by [F2] and L=max⁡{1,ω2}, so L is a global Lipschitz constant for F; and ∥Az∥22=ω4θ2+v2≤L2∥z∥22 with equality at e1 if ω2≥1 and at e2 if 1≥ω2, so the operator norm is ∥A∥=L.

3.1step 1.1step 1.2step 2.1F1given

Instantiation of the linearisation theorem: take n=2, U=R2, the field F of clause 1, the matrix A, the constants L and C:=ω2r/6, and the datum z0 with δ=∥z0∥2>0 and δeLT<r. Steps 1.1, 1.2 and 2.1 verify every hypothesis of [F1], so the maximal solution z of z′=F(z), z(0)=z0 is defined on [0,T], ∥z(t)∥2≤δeLt there, and the two integral bounds of [F1] hold with the displayed C and ∥A∥=L; this is clause 2.

4.1step 3.1F4F5algebra∎

Linearised equation and scope: by [F5] the system y′=Ay reads y1′=y2, y2′=−ω2y1, hence y1′′=−ω2y1, the exactly harmonic linearised pendulum; and ∣z1(t)−y1(t)∣≤∥z(t)−y(t)∥2 because ∥(w1,w2)∥22=w12+w22. The bounds of clause 2 are finite-time and of order δ2; no identification of the nonlinear solution with the harmonic one beyond this estimate is claimed.

Remarks

  • The constant ω2r/6 is exactly the quadratic-remainder constant obtained from the cubic bound ∣sin⁡θ−θ∣≤∣θ∣3/6 on the ball of radius r: the remainder is cubic in the angle and quadratic in the state norm, which is what the linearisation theorem of Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate requires.

  • The scope sentence of clause 3 is part of the claim: the estimate is a conditional consequence of the explicit model θ′′=−ω2sin⁡θ and its stated smallness hypotheses, and it is not a statement about any measured pendulum.

RemarkRemark: AI-adaptedProof: Not applicableOpen item page →

Scope and limits of the verification items on this page

Scope of the verification items

Every item on this page is mathematics: the oscillator, the two-mass spring system, the pendulum field and the mass distributions are mathematical objects with explicit parameters, and no item asserts a physical law, a unit convention or a measured value. The following limits belong to the claims and must be carried by any consumer.

5 · Examples, experiments and counterexamples

None yet.

Sources