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.

✓ 4 results · all verified · 3 also independently AI-judged
Every result on this page is machine-checked by a proof checker and read in full and owner-audited; the judge is an additional, independent cross-model AI review of the proofs. The 1 not AI-judged were verified by owner audit (typically over a confirmed judge false positive), not failures.

Singular and infinite interval quadrature — Examples

1 · Prerequisites

2 · Summary

The examples compute the logarithmic endpoint integral, an interior-pole principal value, and a weakly singular reference-triangle integral. The Gaussian example gives the sinh–sinh transformed identity and a rigorous finite-sum tail bound. Its transformed complex extension fails the uniform strip-decay hypothesis, so no total quadrature error bound is asserted for that example; this records the exact boundary of the analytic-strip theorem.

3 · Logical flowchart

4 · Definitions, theorems and proofs

ExampleConstruction: Literature-sourcedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

The Gaussian integral on the line with a certified sinh-sinh truncation tail

Statement

For f(x)=e−x2 on R, use the sinh–sinh mode of Finite truncated tanh–sinh, exp–sinh, and sinh–sinh trapezoid rules. Its transformed function is Φ(t)=π2cosh⁡t cosh⁡u(t) e−sinh⁡2u(t),u(t)=π2sinh⁡t. Then the absolutely convergent improper Riemann integral satisfies ∫RΦ(t)dt=π. There is T>0 such that ∣Φ(t)∣≤e−e∣t∣ for ∣t∣≥T, so if N±h≥T, the omitted finite-sum tail is at most e−eN−h+e−eN+h. However, the complex extension of this transformed Gaussian does not satisfy the uniform exponential-decay hypothesis of An analytic-strip error bound for the infinite trapezoid sum on any nonzero horizontal substrip. Thus this example certifies its truncation tail but does not claim a full discretization-plus-tail error bound.

Facts & Assumptions

Given: The real sinh–sinh transformation and finite truncation with h>0 and nonnegative integers N−,N+.

[F1]

Finite truncated tanh–sinh, exp–sinh, and sinh–sinh trapezoid rules: x(t)=sinh⁡((π/2)sinh⁡t) and Φ(t)=x′(t)e−x(t)2.

[F2]

The double-exponential maps and their change-of-variable identities: if the transformed function is absolutely integrable as an improper Riemann integral, its integral equals the whole-line target integral.

[F3]

A double-exponential envelope bounds the omitted trapezoid tails: an envelope Ae−Beκ∣t∣ bounds each continuous tail and the corresponding omitted trapezoid samples by A(κB)−1e−BeκNh.

[F5]

A finite double-exponential trapezoid has separate strip and tail errors: a total finite-rule error estimate requires both the separate analytic-strip and sampled-tail hypotheses.

[F6]

An analytic-strip error bound for the infinite trapezoid sum: its analytic-strip hypothesis requires uniform exponential decay on each horizontal edge of some nonzero closed substrip.

[F7]

Complex sine, cosine, hyperbolic sine, and hyperbolic cosine from the complex exponential defines sinh⁡z=(ez−e−z)/2 and cosh⁡z=(ez+e−z)/2. The addition formulas for complex trigonometric and hyperbolic functions, exp⁡(x+iy)=ex(cos⁡y+isin⁡y), ∣exp⁡(x+iy)∣=ex, and eiπ+1=0, and exp⁡(z+w)=exp⁡z exp⁡w, and the complex exponential extends the real exponential give sinh⁡(U+iV)=sinh⁡Ucos⁡V+icosh⁡Usin⁡V and cosh⁡(U+iV)=cosh⁡Ucos⁡V+isinh⁡Usin⁡V for real U,V, by substituting the exponential definitions.

[F8]

The exponential tends to +∞ at +∞ and to 0 at −∞: the real exponential tends to +∞ at +∞ and to zero at −∞. The exponential addition formula exp⁡(x+y)=exp⁡(x)exp⁡(y) and The exponential is positive and satisfies exp⁡(−x)=1/exp⁡(x) give ev+w=evew and e−v=1/ev>0. The real definitions The six hyperbolic functions and their natural domains and Addition formulas, identities, parity, and derivatives of the hyperbolic functions give (sinh⁡t)′=cosh⁡t, (cosh⁡t)′=sinh⁡t, parity, positivity of cosh⁡, and unbounded growth of both hyperbolic functions at +∞.

[F9]

The exponential dominates every fixed nonnegative integer power at +∞: v/ecv→0 as v→+∞ for every fixed c>0.

[F12]

A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion: a continuous real function on a compact interval is Riemann integrable.

[F13]

The zero sets of sine and cosine and the least positive common period 2 pi and Quarter-turn values and shifts by pi/2 and pi: cos⁡0=1, sin⁡(π/2)=1, and neither function has a zero in (0,π/2); continuity and [F11] therefore give sin⁡y,cos⁡y>0 there.

Proof

1.1F1F8F9algebra

On the real axis, x′(t)=(π/2)cosh⁡tcosh⁡u(t)>0. As t→±∞, x(t)→±∞. We show the tail envelope for t→+∞; evenness gives the left tail. For t≥1, sinh⁡t≥cet with c=(1−e−2)/2>0, so u(t)≥c1et for c1=πc/2. For u≥1, sinh⁡u≥ceu, hence x(t)2≥c2e2u(t)≥c2e2c1et. Also u(t)≤(π/4)et and cosh⁡t≤et, cosh⁡u≤eu, so for a fixed constant C′>0 and all sufficiently large t, (π/2)cosh⁡tcosh⁡u(t)≤eC′et. Put v=et. By [F8], v→∞, and [F9] implies c2e2c1v/v→∞. Thus x(t)2 eventually exceeds both 2C′et and 2et. It follows that Φ(t)≤eC′et−x(t)2≤e−x(t)2/2≤e−et for all sufficiently large t.

1.2F6F7F8F10F11F13algebraconstruct

For any proposed substrip halfwidth a>0, fix 0<y<min⁡(a,π/2); then [F13] gives sin⁡y,cos⁡y>0. Extend u(z)=(π/2)sinh⁡z holomorphically; this extension is entire by [F10]. For z=x+iy, write u=U+iV with U=(π/2)sinh⁡xcos⁡y and V=(π/2)cosh⁡xsin⁡y. The continuous function V(x) has positive derivative for x>0 and tends to infinity. For every sufficiently large integer j, its value (j+1/2)π lies strictly between V(0) and some V(R); [F11] gives a solution xj∈(0,R), and strict increase makes it unique. Defining xj as that unique solution avoids any choice of a sequence. Then xj→∞. At those points, [F7] gives Re⁡(sinh⁡2u)=−cosh⁡2U and ∣cosh⁡u∣=sinh⁡U. Since U(xj)→∞, also ∣cosh⁡(xj+iy)∣2=sinh⁡2xj+cos⁡2y≥sinh⁡2xj by [F7] and cosh⁡2xj−sinh⁡2xj=1 from [F8]. Thus it follows that ∣Φ(xj+iy)∣=π2∣cosh⁡(xj+iy)∣sinh⁡U(xj)ecosh⁡2U(xj)⟶∞. Thus no bound Ce−η∣x∣ can hold on that horizontal line; this disproves the uniform strip hypothesis in [F6] for every proposed positive halfwidth.

2.1F2F3F4F12step 1.1

The envelope in step 1.1 makes Φ absolutely integrable on both tails. On every compact interval, ∣Φ∣ is continuous and Riemann integrable by [F12], so Φ is absolutely integrable on R as an improper Riemann integral. Apply [F2] and then [F4] to obtain ∫Φ=π. Applying [F3] with A=B=κ=1 gives the asserted omitted-sample bound when N±h≥T.

3.1F5step 2.1step 1.2∎

By [F5], the bound of the full-error theorem requires the strip premise as well as the tail premise. The real integral and truncation tail are justified in step 2.1, while step 1.2 disproves that strip premise. Thus no total quadrature error bound is asserted.

ExampleConstruction: Literature-sourcedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

An interior Cauchy principal value reduces to a regular integral

Statement

For f(x)=x on [0,1] and c=1/3, the symmetric Cauchy principal value is 1+13log⁡2.

Facts & Assumptions

Given: f(x)=x, a=0, c=1/3, b=1.

[F1]

The Cauchy principal-value subtraction identity: for f∈C1([a,b]), the principal value is the ordinary integral of (f(x)−f(c))/(x−c) plus f(c)log⁡((b−c)/(c−a)).

[F2]
[F3]

Higher derivatives and the classes Ck and C∞: the linear function is C1 on [0,1].

Proof

1.1givenF3algebra

The function f(x)=x is C1 by [F3] and f(c)=1/3. For every x≠c, (f(x)−f(c))/(x−c)=1, and its continuous extension at c is also 1.

2.1F1F2step 1.1algebra∎

Apply [F1]. The ordinary integral of this quotient is ∫011dx=1, and (b−c)/(c−a)=(2/3)/(1/3)=2. Thus the principal value equals 1+(1/3)log⁡2.

ExampleConstruction: Literature-sourcedVerification: AI-generatedprecheck passOpen item page →

A 1 over r singularity on the reference triangle

Statement

On T={(x,y):0≤y≤x≤1}, the weakly singular integral ∫Tdx dyx2+y2=log⁡(1+2). Under The reference-triangle Duffy map from a vertex, its transformed amplitude is G(r,s)=(1+s2)−1/2 with radial weight r0. For every n≥1, Weighted product Gauss quadrature after the Duffy transform applies with α=1, and The Duffy weighted product rule has an explicit bound using pure coordinate derivatives bounds its error by K0∥∂r2nG∥∞+K0∥∂s2nG∥∞ with the constants specified there; here ∂r2nG=0.

Facts & Assumptions

Given: The reference triangle T and integrand (x2+y2)−1/2 away from its vertex.

[F1]

The Duffy map turns a vertex singularity into an integrable radial weight: for exponent α=1 and amplitude q=1, the transformed radial weight is r0 and G(r,s)=(1+s2)−1/2.

[F2]

The Duffy weighted product rule has an explicit bound using pure coordinate derivatives: the weighted product rule error is bounded by the two pure derivative norms with its stated constants.

Proof

1.1F1F3algebra

For G(r,s)=(1+s2)−1/2, the transformed integral in [F1] is ∫01∫01G(r,s)dsdr. An antiderivative in s is the logarithm in [F3], so the integral equals log⁡(1+2)−log⁡1=log⁡(1+2).

2.1F2F4F5givenalgebra∎

The function G is smooth on the compact square by [F4], and its derivatives are bounded there by [F5]. The Duffy weighted product rule therefore has the error bound [F2]; since G is independent of r, its positive-order r derivative vanishes.

ExampleConstruction: Literature-sourcedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

Gauss–Jacobi regularization for x to the minus one half times log x

Statement

The improper integral is ∫01x−1/2log⁡x dx=−4. For example, Finite endpoint power substitution followed by Gauss–Jacobi quadrature applies with n=1, α=−1/2, β=0, ma=1, mb=0, h(x)=1, qa=7, and qb=4. The transformed functions are C2 and Error bound after endpoint power regularization gives an explicit error bound C1(∥FL′′∥∞+∥FR′′∥∞) for the returned two-node sum.

Facts & Assumptions

Given: The integrand f(x)=x−1/2log⁡x on (0,1] and the specified endpoint-regularization inputs.

[F1]

The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t: log⁡ is differentiable on (0,∞) with derivative 1/x.

[F2]

A sufficiently high power regularizes finite logarithmic endpoint factors: tδP(log⁡t) extends to Cr([0,1]) when δ>r and P is a polynomial.

[F3]

Finite endpoint power substitution followed by Gauss–Jacobi quadrature: the left and right substitution orders qa,qb and the shifted Gauss–Jacobi sums define the stated finite output.

[F4]

Error bound after endpoint power regularization: when the substituted functions are C2n, the improper integral exists and the error is bounded by the displayed derivative constant.

Proof

1.1F1F2F5F6algebra

A primitive on (0,1] is P(x)=2xlog⁡x−4x by [F1] and [F6]. The endpoint term xlog⁡x tends to zero by the power-regularity lemma [F2] with r=0, δ=1/2, and polynomial z. Hence P(x)→0 as x↓0. Applying [F5] on [ε,1] and taking ε↓0 gives the improper integral P(1)−0=−4.

1.2F2F3givenalgebra

For the left substitution, δL=qa(α+1)−1=7/2−1=5/2>2; its transformed function is a constant times t5/2 times a polynomial of degree one in log⁡t and a smooth factor bounded with two continuous derivatives. For the right substitution, δR=qb(β+1)−1=3>2 and the remaining factors are smooth on [0,1]. Therefore both transformed functions are C2 by [F2], and [F3] applies.

2.1F3F4step 1.1step 1.2∎

Apply [F4] with n=1 to the functions in step 1.2. Its error constant is C1 and the output is the stated sum, so the claimed error bound follows. The value −4 from step 1.1 is the integral being approximated.

5 · Examples, counterexamples and false statements

None yet.

Sources