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 Approximation — Examples

1 · Prerequisites

2 · Summary

The companion page to cm-math-approximation gathers explicit worked checks for its oscillator, forcing, linearization, averaging and adiabatic results. Everything here is mathematics: the items carry domain: mathematics, depend only on mathematical suppliers, and make no physical or empirical claim.

Five worked initial value problems for the homogeneous oscillator across the discriminant, degenerate and repulsive cases solves five initial value problems of the homogeneous model, one in each discriminant case and in the degenerate k≤0 cases; Resonant forcing at Ω=ω0 and the two regimes of the amplitude maximum checks the resonant solution 12tsin⁡t of x′′+x=cos⁡t and evaluates the amplitude maximum in its two regimes; Constant forcing: free fall and the approach to terminal velocity under linear drag treats constant forcing with the convolution formula, including the approach to terminal velocity under linear drag. Finite-time linearization bound for the planar pendulum with explicit constants verifies the quadratic remainder and growth constants for F(θ,ω)=(ω,−sin⁡θ) on the ball of radius 1/2 and evaluates the resulting error bound for a small datum. Averaging estimate on [0,100] for the scalar equation x′=εcos⁡t computes the average and deviation for f(t,x)=cos⁡t, checks the trapping hypothesis, and compares the theorem's bound with the exact solution on [0,100]; An exactly solvable adiabatic branch: K=I+εIcos⁡θ and comparison with the drift bound exhibits the exactly invariant quantity I(1+εcos⁡θ) for the branch K=I+εIcos⁡θ, verifies the trapping in the stated action window, and evaluates the drift bound. Continuity alone does not give uniqueness: x′=∣x∣ through the origin is the counterexample x′=∣x∣ through the origin, showing that continuity of the field cannot replace the Lipschitz hypothesis in the uniqueness results, and A bounded field whose solution leaves the domain in finite time: x′=v, v′=0 on (0,∞)×(−2,2) shows a solution of a bounded field leaving its domain at a finite time.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

None yet.

5 · Examples, experiments and counterexamples

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Five worked initial value problems for the homogeneous oscillator across the discriminant, degenerate and repulsive cases

Example

In each line, m=1 and the initial datum (x0,v0) is prescribed for the homogeneous model with time domain I=R and zero forcing, x′′+cx′+kx=0 (The damped driven oscillator model: second-order equation, phase-plane system and characteristic data). The solution is read off from the classification theorem and the degenerate-case proposition, and its initial values are checked explicitly.

  • (i) k=4, c=0, (x0,v0)=(1,0): the solution is x(t)=cos⁡(2t).
  • (ii) k=1, c=2, (x0,v0)=(1,0): the solution is x(t)=(1+t)e−t.
  • (iii) k=1, c=3, (x0,v0)=(1,0): with λ±=(−3±5)/2 the solution is x(t)=(−λ−eλ+t+λ+eλ−t)/5.
  • (iv) k=0, c=2, (x0,v0)=(0,1): the solution is x(t)=12(1−e−2t).
  • (v) k=−1, c=0, (x0,v0)=(1,0): the solution is x(t)=12(et+e−t), which is unbounded on [0,∞).

Facts & Assumptions

Given: The homogeneous model of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with m=1 and the parameters of the five lines.

[L1]

For k>0 the solutions are classified by the sign of Δ=c2−4mk as in The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data: distinct real roots in case 1, a double root λ=−c/(2m) in case 2, and the family e−γt(Acos⁡ωt+Bsin⁡ωt) with γ=c/(2m), ω=k/m−γ2 in case 3; clause 4 gives the coefficient formulas A=(v0−λ−x0)/(λ+−λ−), B=(λ+x0−v0)/(λ+−λ−) in case 1; A=x0, B=v0−λx0 in case 2; A=x0, B=(v0+γx0)/ω in case 3. [given]

[L2]

For k≤0 the solutions are those of The degenerate and repulsive oscillator: solutions for k=0 and the explicit growing solution for k<0: for k=0, c>0 they are A+Be−ct/m with A=x0+v0m/c, B=−v0m/c; for k<0 they are Aeλ+t+Beλ−t with λ±=(−c±c2−4mk)/(2m), λ−<0<λ+, and eλ+t→+∞. [given]

[L3]

(cos⁡t)′=−sin⁡t, (sin⁡t)′=cos⁡t, and (eμt)′=μeμt for real μ. [thm-sine-and-cosine-derivatives, thm-derivative-of-exponential]

Verification

technique · direct; compute the discriminant and coefficients in each line and check the initial data
1.1givenL1L3algebra

Line (i): Δ=0−16<0 and γ=0, ω=4=2, so by [L1] case 3 the solution is Acos⁡2t+Bsin⁡2t with A=x0=1, B=(v0+γx0)/ω=0; by [L3] the candidate cos⁡(2t) has value 1 and derivative −2sin⁡(2t) vanishing at 0, so it is the solution.

1.2givenL1L3algebra

Line (ii): Δ=4−4=0, so the double root is λ=−1 and the family is (A+Bt)e−t with A=x0=1, B=v0−λx0=1; the candidate (1+t)e−t has the value 1 and derivative e−t−(1+t)e−t=−te−t equal to 0 at t=0, matching (x0,v0)=(1,0).

1.3givenL1algebra

Line (iii): Δ=5>0 and λ±=(−3±5)/2 with λ+−λ−=5 and λ++λ−=−3; the coefficients of [L1] case 1 are A=(v0−λ−x0)/5=−λ−/5 and B=(λ+x0−v0)/5=λ+/5, so the solution is x(t)=Aeλ+t+Beλ−t=(−λ−eλ+t+λ+eλ−t)/5. Then x(0)=(λ+−λ−)/5=1 and x′(0)=λ+A+λ−B=(−λ+λ−+λ−λ+)/5=0, matching (x0,v0)=(1,0).

1.4givenL2L3algebra

Line (iv): k=0, c=2>0, and [L2] gives x(t)=A+Be−2t with A=0+1/2=1/2, B=−1/2, that is x(t)=12(1−e−2t); its derivative is e−2t by [L3], so x(0)=0 and x′(0)=1.

2.1givenL2L3algebra∎

Line (v): k=−1<0 and λ±=±1 by [L2]; the solution with data (1,0) is 12et+12e−t (the coefficients A=B=1/2 solve A+B=1, A−B=0), and it is unbounded on [0,∞) because et→+∞, in accordance with [L2].

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Resonant forcing at Ω=ω0 and the two regimes of the amplitude maximum

Example

Consider the forced model m x′′+c x′+k x=f0cos⁡(Ωt) on time domain I=R with the values and conclusions below (The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, Sinusoidal forcing: the steady response, its amplitude and phase, and resonance at Ω=k/m, The amplitude of the sinusoidally forced damped oscillator: location and value of its maximum).

  • (i) Resonance: m=k=f0=1, c=0, Ω=ω0=1. Then xp(t)=12tsin⁡t is a solution of x′′+x=cos⁡t, it satisfies xp(0)=0, xp′(0)=0, and ∣xp∣ is unbounded on [0,∞).
  • (ii) Amplitude peak above zero frequency: m=k=f0=1, c=1. Then u∗=1−12=12>0, the maximum of A(Ω)=1/(1−Ω2)2+Ω2 is attained at Ω∗=1/2, and A(Ω∗)=2/3, while A(0)=A(1)=1.
  • (iii) Amplitude peak at zero frequency: m=k=f0=1, c=2. Then u∗=1−2=−1≤0, and the maximum of A(Ω)=1/(1−Ω2)2+4Ω2 on [0,∞) is attained at Ω=0 with value 1, while A(1)=1/2.

Facts & Assumptions

[L1]

In the resonant case c=0, Ω=ω0=k/m the function xp(t)=f0tsin⁡(ω0t)/(2mω0) is a solution of the forced model on R and is unbounded on [0,∞) when f0≠0 (Sinusoidal forcing: the steady response, its amplitude and phase, and resonance at Ω=k/m clause 2). [given]

[L2]

If c>0 the amplitude of the steady response is A(Ω)=∣f0∣/(k−mΩ2)2+c2Ω2; its maximum is attained at Ω∗=u∗ when u∗=k/m−c2/(2m2)>0, with value 2m∣f0∣/(c4mk−c2), and at Ω=0 when u∗≤0 (The amplitude of the sinusoidally forced damped oscillator: location and value of its maximum; the connection to the steady response is clause 1 of Sinusoidal forcing: the steady response, its amplitude and phase, and resonance at Ω=k/m). [given]

[L3]

(sin⁡t)′=cos⁡t, (cos⁡t)′=−sin⁡t; sums and products are differentiated by the algebra of derivatives and the chain rule. [thm-sine-and-cosine-derivatives, thm-algebra-of-derivatives, thm-chain-rule]

Verification

technique · direct; substitute the resonant candidate, and evaluate the amplitude formula in the two regimes of the preceding proposition
1.1givenL1L3algebra

Line (i): with m=ω0=f0=1 the candidate is xp(t)=12tsin⁡t; by [L3], xp′(t)=12(sin⁡t+tcos⁡t) and xp′′(t)=12(2cos⁡t−tsin⁡t)=cos⁡t−12tsin⁡t, so xp′′+xp=cos⁡t, as required, and xp(0)=0=xp′(0). At tk=π/2+2πk one has sin⁡tk=1 and ∣xp(tk)∣=tk/2→∞, so xp is unbounded on [0,∞), in agreement with [L1].

1.2givenL2algebra

Line (ii): here u∗=1−1/2=1/2>0 and 4mk−c2=3, so [L2] locates the maximum at Ω∗=1/2 with value 2/3; directly, (1−Ω∗2)2+Ω∗2=(1/2)2+1/2=3/4, so A(Ω∗)=2/3, while A(0)=1 and A(1)=1/0+1=1, both below 2/3.

2.1givenL2algebra∎

Line (iii): here u∗=1−2=−1≤0, so [L2] locates the maximum at 0 with value A(0)=1/1=1; directly, A(1)=1/(0)2+4=1/2<1, and the same formula gives A(Ω)≤1 for all Ω≥0 because (1−Ω2)2+4Ω2=1+2Ω2+Ω4≥1.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Constant forcing: free fall and the approach to terminal velocity under linear drag

Example

Let m>0, g∈R and consider the two constant-forcing models with time domain I=R (the half-line solution in (ii) is the restriction of its global solution) (The damped driven oscillator model: second-order equation, phase-plane system and characteristic data).

  • (i) Uniform force: x′′=g with data x(0)=x0, x′(0)=v0. The unique solution is x(t)=x0+v0t+12gt2.
  • (ii) Linear drag: m=1, c=2, g=1, that is x′′+2x′=1, with data (x0,v0)=(0,0). The unique solution on [0,∞) is x(t)=12t−14(1−e−2t), and its velocity x′(t)=12(1−e−2t) increases to the limiting value 12.

Facts & Assumptions

Given: The parameters m,g and the data of the two lines; the forced oscillator theorem applies in both since k=0≥0 and the forcing is constant, hence continuous.

[L1]

For m>0, c≥0, k≥0 and continuous f, the unique solution of mx′′+cx′+kx=f with data (x0,v0) at t0=0 is x=xh+∫0tg(t−s)f(s)ds, where xh is the homogeneous solution with the data and g is the homogeneous kernel with g(0)=0, g′(0)=1/m; for k=c=0 one has g(u)=u/m, and for k=0, c>0 one has g(u)=(1−e−cu/m)/c (The forced oscillator with k≥0: unique global solution and the convolution (Green's function) formula, The degenerate and repulsive oscillator: solutions for k=0 and the explicit growing solution for k<0). [given]

[L2]

For k=0 the homogeneous solutions are xh(t)=A+Be−ct/m when c>0 and xh(t)=A+Bt when c=0, with A=x0+v0m/c, B=−v0m/c respectively A=x0, B=v0 (The degenerate and repulsive oscillator: solutions for k=0 and the explicit growing solution for k<0). [given]

[L3]

(eμt)′=μeμt for real μ, and sums and scalar multiples are differentiated termwise. [thm-derivative-of-exponential, thm-algebra-of-derivatives]

Verification

technique · direct; read the solutions off the convolution formula in the two degenerate cases and differentiate
1.1givenL1L2algebra

Line (i): here c=k=0, so [L1] gives g(u)=u/m and xh(t)=x0+v0t by [L2]; hence x(t)=x0+v0t+∫0tt−sm mg ds=x0+v0t+g∫0t(t−s) ds=x0+v0t+12gt2, because ∫0t(t−s)ds=t2/2.

1.2givenL1L2algebra

Line (ii): here m=1, c=2, k=0 and f≡1; by [L1], g(u)=(1−e−2u)/2, and the datum (0,0) gives xh=0 by [L2]. Hence x(t)=∫0t1−e−2(t−s)2 ds=t2−14(1−e−2t).

2.1step 1.1step 1.2L3algebra∎

Velocity: by [L3], x′(t)=12−14⋅2e−2t=12(1−e−2t), so x′(0)=0 and x′(t)→12 as t→+∞ because e−2t→0; moreover x′′(t)=e−2t, so x′′+2x′=e−2t+1−e−2t=1, confirming both lines satisfy their equations with the prescribed data.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Finite-time linearization bound for the planar pendulum with explicit constants

Example

Let l=g=1 and consider the planar pendulum field

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

Let A:R2→R2 be the linear map A(θ,ω):=(ω,−θ), and let A=(01−10)∈M2(R) be its standard-coordinate matrix; use the matrix action of A from Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate. On the closed ball B‾(0,12) one has

∥F(z)−Az∥2≤112∥z∥22,∥F(z)∥2≤∥z∥2,

so the linearization theorem applies with r=12, C=112, L=1. For the datum z0=(1100,0) and T=1 one has δeLT=e100<3100<120<12, and the theorem bounds the linearization error by

∥z(t)−y(t)∥2≤δ212∫01e(1−s)e2s ds<10−4(0≤t≤1),

where y solves y′=Ay, y(0)=z0; in particular the linear solution stays within 10−4 of the nonlinear one on [0,1].

Facts & Assumptions

Given: The field F, the matrix A, the radius r=1/2, the datum z0=(1/100,0) and the horizon T=1.

[L1]

The linearization theorem (Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate): if U⊆Rn is open, 0∈U, F:U→Rn is continuous and locally Lipschitz in the state variable, r,T>0, C,L≥0, B‾(0,r)⊆U, δ=∥z0∥2>0, and ∥F(z)−Az∥2≤C∥z∥22 and ∥F(z)∥2≤L∥z∥2 for ∥z∥2≤r and δeLT<r, then the nonlinear solution exists on [0,T] and ∥z(t)−y(t)∥2≤Cδ2∫0te∥A∥(t−s)e2Lsds for t∈[0,T]. [given]

[L2]

Taylor remainder bound: if f has derivatives through order n+1 on the closed interval between a and x, with endpoint continuity, and ∣f(n+1)∣≤M there, then ∣Rn,af(x)∣≤M∣x−a∣n+1/(n+1)! (A uniform derivative bound gives a uniform Taylor remainder bound); with f=sin⁡, n=2, a=0 and M=1 this gives ∣sin⁡θ−θ∣≤∣θ∣3/6, since ∣sin⁡′′′∣=∣cos⁡∣≤1 (The derivatives of sine and cosine are cosine and minus sine, Signs, monotonicity intervals, and ranges of sine and cosine). [given]

[L3]

Operator norm: ∥A∥=sup⁡∥u∥2=1∥Au∥2 (The operator norm as the least bound and as the unit-sphere or unit-ball supremum); for A=(01−10) and u=(u1,u2) one has Au=(u2,−u1), so ∥Au∥2=∥u∥2 and ∥A∥=1. [given]

Verification

technique · direct; estimate the two remainder terms of the pendulum field on the ball and substitute the explicit constants into the theorem
1.1givenL2algebra

Quadratic remainder: for z=(θ,ω) with ∥z∥2≤1/2, F(z)−Az=(0, −sin⁡θ)−(0,−θ)=(0, θ−sin⁡θ), so ∥F(z)−Az∥2=∣sin⁡θ−θ∣≤∣θ∣3/6≤12∣θ∣2/6≤112∥z∥22 by [L2] and ∣θ∣≤∥z∥2≤1/2; hence C=1/12 is admissible.

1.2givenL1L3L4L6algebra

Regularity and growth: U=R2 is open and contains B‾(0,1/2). The differentiable sine and coordinate functions make F continuous. For z=(θ,ω) and w=(η,ν), [L6] and [L4] give ∣sin⁡θ−sin⁡η∣≤∣θ−η∣, hence ∥F(z)−F(w)∥22=(ω−ν)2+(sin⁡θ−sin⁡η)2≤∥z−w∥22. Thus F is globally 1-Lipschitz and in particular locally Lipschitz in the state variable. Growth bound: if θ≠0, the mean value theorem (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)) applied to sin⁡ between 0 and θ, with ∣cos⁡∣≤1 from [L4], gives ∣sin⁡θ∣=∣sin⁡θ−sin⁡0∣≤∣θ∣; this inequality is immediate at θ=0. Thus ∥F(z)∥22=ω2+sin⁡2θ≤ω2+θ2=∥z∥22, so L=1 is admissible. With δ=1/100 and T=1, [L4] gives δeLT=e/100<3/100<1/20<1/2=r, so the hypotheses of [L1] hold.

2.1step 1.1step 1.2L1L3L4L5algebra∎

Substitute: by [L1] and [L3], for 0≤t≤1 the error is at most 112δ2∫0te(1−s)e2sds≤112δ2∫01e(1−s)e2sds, since et−s≤e1−s for s≤t≤1 and the resulting integrand is positive by [L4]. The addition formula gives e(1−s)e2s=e1+s; since 0≤s≤1, monotonicity gives e1+s≤e2<9 by [L4]. Monotonicity of the integral [L5] therefore yields ∫01e1+sds≤9, and the error is at most 112⋅10−4⋅9=34⋅10−4<10−4.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Averaging estimate on [0,100] for the scalar equation x′=εcos⁡t

Example

Take n=1, P=2π, O=R, f(t,x)=cos⁡t (independent of x), x0=0, and the trapping data R=1, ρ=12, A=1. Then f is continuous, C1 in x with Dxf=0, and 2π-periodic in t; on R×B‾(0,1) one may take M=1 and L=0. The average and deviation of Periodic averaging data: the average field, the oscillatory part and the periodic deviation are

fˉ=0,g(t,x)=sin⁡t(0≤t≤2π, extended 2π-periodically).

The averaged system is y′=0, so y(t)≡0 solves it and stays in B‾(0,12). With B=2MP=4π, B1=2LP=0 and K=eLA[B+(LB+B1M)A]=4π, the smallness hypothesis reads ε⋅4π<12; for ε=1100 it holds because εK=π25<12. The averaging theorem then gives

∣x(t)−y(t)∣=∣x(t)∣≤π25(0≤t≤100).

The exact solution is x(t)=1100sin⁡t, whose modulus is at most 1100; the theorem's bound is valid but not sharp here.

Facts & Assumptions

Given: The data above and the scalar ODE x′=εcos⁡t, x(0)=0.

[L1]

For this f, fˉ(x)=12π∫02πcos⁡s ds=0 and g(t,x)=∫0tcos⁡s ds=sin⁡t on [0,2π], extended periodically; this is the construction of Periodic averaging data: the average field, the oscillatory part and the periodic deviation, evaluated with the fundamental theorem for continuous integrands (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G, The derivatives of sine and cosine are cosine and minus sine); the endpoint values sin⁡0=sin⁡(2π)=0 and periodicity used here are supplied by [L4]. [given]

[L2]

The averaging theorem (Averaging for fast periodic oscillations: order-ε agreement on the O(1/ε) time scale): under its hypotheses with a solution y:[0,A/ε]→B‾(x0,R−ρ) of the averaged system and εK<ρ, the solution x of x′=εf(t,x), x(0)=x0 exists on [0,A/ε] and satisfies ∥x(t)−y(t)∥2≤εK there, with K=eLA[B+(LB+B1M)A], B=2MP, B1=2LP. [given]

[L3]
[L4]

Sine and cosine are 2π-periodic, and sin⁡(2π)=sin⁡0=0 by their zero set (The zero sets of sine and cosine and the least positive common period 2 pi). Their ranges are [−1,1] (Signs, monotonicity intervals, and ranges of sine and cosine); in particular ∣cos⁡t∣,∣sin⁡t∣≤1. Their differentiability in [L1] gives continuity.

Verification

technique · direct; evaluate the averaging data, check the hypotheses, and compare the bound with the exact solution
1.1givenL1L3L4algebra

By [L4], the continuous field f(t,x)=cos⁡t is 2π-periodic and bounded by one; its state derivative is zero, so M=1 and L=0 meet the hypotheses on O=R. By [L1], fˉ=0 and g(t,x)=sin⁡t, so the averaged system is y′=0; the constant curve y≡0 is its solution with y(0)=0 on [0,100], and ∥y(t)−0∥=0≤R−ρ=12. The constants are M=1, L=0, B=4π, B1=0, K=e0[4π+0]=4π, and with ε=1100, [L3] gives εK=π25<425<1/2=ρ.

2.1step 1.1L1L2L4algebra∎

By [L2] the solution x of x′=1100cos⁡t, x(0)=0 exists on [0,100] and satisfies ∣x(t)∣≤π25 there. Directly, the primitive sin⁡t gives x(t)=1100sin⁡t by [L1], so ∣x(t)∣≤1100 by [L4]; the stated bound holds and is not attained.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

An exactly solvable adiabatic branch: K=I+εIcos⁡θ and comparison with the drift bound

Example

Take J=R, S=R, h(I,s)=I and k(I,θ,s)=Icos⁡θ in the datum of A one-frequency parameter-dependent action-angle branch, its transformed Hamiltonian and its drift generator, with ε=110, A=2, ω0=1, K=[6/5,19/10], S0=[0,2], s0=0 and initial datum I(0)=3/2, θ(0)=0. The branch system is

I′=εIsin⁡θ,θ′=1+εcos⁡θ,s′=ε,

so s(t)=εt stays in S0=[0,2] on [0,A/ε]=[0,20], and I(t) stays in K for all t≥0. The quantity J:=I(1+εcos⁡θ) is exactly invariant, the solution stays in K×S0 on [0,A/ε], and

∣I(t)−I(0)∣≤13(t≥0),

while the drift bound of Adiabatic action drift bound: order-ε change of the action on the O(1/ε) time scale evaluates in this instance to ∣I(t)−I(0)∣≤1950.

Facts & Assumptions

Given: The datum h(I,s)=I, k(I,θ,s)=Icos⁡θ, ε=1/10, and the sets and data above.

[L1]

The branch vector field is I′=−ε∂θk, θ′=∂Ih+ε∂Ik, s′=ε; here ∂Ih=1, ∂θk=−Isin⁡θ and ∂Ik=cos⁡θ, so the field is as displayed (A one-frequency parameter-dependent action-angle branch, its transformed Hamiltonian and its drift generator). [given]

[L2]

The drift bound (Adiabatic action drift bound: order-ε change of the action on the O(1/ε) time scale): if ∣∂Ih∣≥ω0 and the branch solution stays in the compact sets K×S0, then ∣I(t)−I(0)∣≤ε[2sup⁡K×[0,2π]×S0∣g∣+Asup⁡K×[0,2π]×S0∣∂Ig ∂θk−∂θg ∂Ik−∂sg∣], where g=(k−⟨k⟩)/∂Ih and ⟨k⟩ is the θ-average of k. [given]

[L4]

∣sin⁡θ∣≤1 and ∣cos⁡θ∣≤1 for all real θ (Signs, monotonicity intervals, and ranges of sine and cosine); the mean value theorem then gives ∣sin⁡θ−sin⁡ϑ∣≤∣θ−ϑ∣ and ∣cos⁡θ−cos⁡ϑ∣≤∣θ−ϑ∣ (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)).

[L5]

Local Lipschitz continuity in the state variable is the neighborhood-uniform state-Lipschitz property of Local Lipschitz continuity in the state variable, locally uniform in time and parameters. For a continuous field with this property on an open ODE domain, there is a unique maximal solution through each datum, and a solution of an autonomous field whose state remains in a compact subset of its open state domain near a finite endpoint extends past that endpoint (First-order systems, initial value problems, and solutions on intervals, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval, Extensions, maximal solutions, maximal intervals, and global solutions, Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it).

[L6]

If a differentiable vector-valued curve on a compact interval has derivative norm at most C, then its endpoint difference has norm at most C times the interval length (If f:[a,b]→Rm is differentiable with integrable f′ then ∫abf′=f(b)−f(a); and a bounded derivative makes f Lipschitz).

[L7]

The zero-derivative constancy fact used below is A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant: a differentiable real function on an interval whose derivative vanishes is constant.

Verification

technique · direct; exhibit an exact invariant of the branch system, check the trapping, and evaluate the theorem's constant
1.1L1L4L5algebra

Write the branch field as V(I,θ,s)=(110Isin⁡θ,1+110cos⁡θ,110). On each cube CN=[−N,N]3 with N>0, [L4] gives, for p=(I,θ,s) and q=(I′,ϑ,s′), ∣Isin⁡θ−I′sin⁡ϑ∣≤∣I−I′∣+N∣θ−ϑ∣,∣cos⁡θ−cos⁡ϑ∣≤∣θ−ϑ∣. Hence ∥V(p)−V(q)∥2≤110(∣I−I′∣+(N+1)∣θ−ϑ∣)≤N+210∥p−q∥2. Thus V is continuous and locally Lipschitz on its open domain R3; by [L5] there is a unique maximal solution (I,θ,s):(α,β)→R3 through the stated datum.

2.1step 1.1L1L3L7algebra

Along the maximal solution of step 1.1, [L1] gives I′=110Isin⁡θ, θ′=1+110cos⁡θ, s′=110. Put J(t):=I(t)(1+110cos⁡θ(t)). By [L3], J′=110Isin⁡θ (1+110cos⁡θ)+I (−110sin⁡θ)(1+110cos⁡θ)=0, so [L7] makes J constant; with I(0)=3/2, θ(0)=0 we get J≡3320.

3.1step 1.1step 2.1L2L4L5L6algebra

Trapping and global forward existence: on [0,β), step 2.1 and [L4] give 9/10≤1+110cos⁡θ(t)≤11/10 and hence 32=3320⋅1011≤I(t)=33/201+110cos⁡θ(t)≤3320⋅109=116. Also ∣θ′(t)∣≤11/10 and ∣s′(t)∣=1/10, so [L6] gives ∣θ(t)∣≤11t/10 and s(t)=t/10. If β<+∞, the state curve near β stays in the compact subset [32,116]×[−11β10,11β10]×[0,β10] of R3; the autonomous compact-continuation clause of [L5] extends the maximal solution past β, a contradiction. Hence β=+∞. Since s(t)=t/10, on [0,20]=[0,A/ε] we have s(t)∈S0, and the bounds above place I(t) strictly inside K=[6/5,19/10]. Thus the solution stays in K×S0 on the required horizon, and [L2] applies.

4.1step 2.1step 3.1L2algebra∎

Evaluate the theorem's constant: ⟨k⟩(I,s)=12π∫02πIcos⁡θ dθ=0, so g=Icos⁡θ and sup⁡K×[0,2π]×S0∣g∣=1910. Moreover ∂Ig=cos⁡θ, ∂θg=−Isin⁡θ, ∂sg=0, and with ∂θk=−Isin⁡θ, ∂Ik=cos⁡θ one has ∂Ig ∂θk−∂θg ∂Ik−∂sg=−Icos⁡θsin⁡θ+Isin⁡θcos⁡θ=0. The theorem thus gives ∣I(t)−I(0)∣≤110⋅2⋅1910=1950, while the exact invariant gives ∣I(t)−I(0)∣=110⋅32 ∣cos⁡θ(0)−cos⁡θ(t)∣1+110cos⁡θ(t)≤110⋅32⋅29/10=13, consistent with the bound.

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

Continuity alone does not give uniqueness: x′=∣x∣ through the origin

Statement refuted

Let F:R→R be continuous. Then for every x0∈R the initial value problem x′=F(x), x(0)=x0 has at most one solution on every interval containing 0.

Facts & Assumptions

Given: The field F(x)=∣x∣ and the initial value problem x′=F(x), x(0)=0.

Counterexample

technique · direct; exhibit two distinct solutions through the same initial datum
1.1givenL2

The map F(x)=∣x∣ is continuous on R as the composition of u↦∣u∣ with the square root, by [L2], so the hypothesis of the refuted statement holds.

1.2givenL1

The zero function x0≡0 is a solution: it is differentiable with derivative 0=F(0) everywhere.

1.3givenL3algebra

Define x1:R→R by x1(t)=0 for t≤1 and x1(t)=(t−1)2/4 for t≥1. Then x1 is continuous, and it is differentiable at t=1 with derivative 0 from both sides: for h>0, (h2/4)/h=h/4→0, while the left difference quotient is 0. For t≠1 the derivative is 0 (for t<1) or (t−1)/2 (for t>1) by [L3]. Hence x1 is differentiable on R with x1(0)=0.

2.1step 1.2step 1.3L1L2algebra∎

For t>1, x1(t)=(t−1)2/4≥0 and (t−1)/2>0, so uniqueness of the nonnegative square root in [L2] gives ∣x1(t)∣=(t−1)/2=x1′(t); for t≤1, x1(t)=0 and x1′(t)=0=∣0∣. Hence x1 is a solution of x′=F(x) through (0,0) by [L1], and it differs from the zero solution on (1,∞): the initial value problem has two distinct solutions on every interval containing [0,2], refuting the statement.

Remarks

  • Why the Lipschitz hypothesis matters. The local uniqueness lemma Continuous and Lipschitz dependence on the initial condition on a Picard-Lindelof cylinder assumes a state-Lipschitz constant; here F fails that hypothesis at the equilibrium: ∣F(x)−F(0)∣/∣x∣=1/∣x∣→∞ as x→0. The counterexample shows that continuity alone cannot replace it, so the Lipschitz hypothesis in the oscillator existence and uniqueness results is not redundant.
  • No physical claim. This is a statement about an abstract scalar ODE; it makes no assertion about any mechanical model.
ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

A bounded field whose solution leaves the domain in finite time: x′=v, v′=0 on (0,∞)×(−2,2)

Example

Let U=(0,∞)×(−2,2)⊆R2 and let F:U→R2, F(x,v)=(v,0), so that the first-order system (x,v)′=F(x,v) is the free particle with velocity v (The local flow of an autonomous first-order system). The solution of z′=F(z), z(0)=(1,−1) is

z(t)=(x(t),v(t))=(1−t, −1),t∈(−∞,1).

Its maximal interval of existence is (−∞,1): at t=1 the position reaches the excluded boundary value x=0, although ∥F(x,v)∥=∣v∣<2 is bounded on U and the field extends smoothly to the plane. Hence a finite maximal interval does not require the field to blow up; the domain can force it.

Facts & Assumptions

Given: The domain U=(0,∞)×(−2,2), the field F(x,v)=(v,0), and the datum z0=(1,−1).

[L1]

F is continuous and locally Lipschitz in the state variable; the flow, maximal solutions and maximal intervals are those of The local flow of an autonomous first-order system and Extensions, maximal solutions, maximal intervals, and global solutions; every initial value problem has a unique maximal solution, and every other solution through the same data is a restriction of it (Every Picard–Lindelöf initial value problem has one maximal solution on an open interval, First-order systems, initial value problems, and solutions on intervals). A solution on an interval I must satisfy (t,z(t))∈R×U for every t∈I (Intervals of R: the nine order-convex forms, nondegeneracy, and length). [given]

Verification

technique · direct; verify the explicit solution, then show no extension past $t=1$ can meet the domain condition
1.1givenL1algebra

The curve z(t)=(1−t,−1) is defined for t<1 and lies in R×U there because 1−t>0 and −1∈(−2,2); it is differentiable with z′(t)=(−1,0)=F(z(t)) for every t<1, and z(0)=(1,−1). So it is a solution on (−∞,1).

2.1step 1.1L1algebra∎

Maximality: suppose a solution ζ through (1,−1) were defined on an interval containing some t0>1. By [L1] it agrees with z on the overlap of their intervals and is continuous at t=1, so ζ(1)=lim⁡t↑1(1−t,−1)=(0,−1), i.e. its position is 0 at time 1; but then (1,ζ(1))=(1,(0,−1)) does not lie in R×U, contradicting the domain condition for a solution on that interval. Hence no solution through the datum is defined at t0>1, while z is defined at every t<1; the maximal interval is exactly (−∞,1). In particular the maximal solution escapes the domain at the finite time 1 with ∥F(x,v)∥=∣v∣<2 bounded on U, so boundedness of the field alone does not guarantee a global solution.

Sources