Problem 13 definitions — First-passage time of Brownian motion to an exponentially decaying boundary
Definitionrybin2026_p13_exponential_boundaryElementaryExpr is a parameter-free inductive type whose elements are finite expression trees. It has a nullary variable expression; a rational-constant expression for every ; a real-constant expression for every ; binary addition, subtraction, multiplication, and division expressions formed from any two expressions; and unary negation, exponential, logarithm, square-root, sine, cosine, and normal-density expressions formed from any one expression. There are no side conditions on any constructor and no equations identifying distinct constructor trees.
For every elementary expression and every real input , ElementaryExpr.eval assigns a real value, denoted , recursively as follows: the variable evaluates to ; a rational constant evaluates to its image in , independently of ; a real constant evaluates to , independently of ; addition, subtraction, multiplication, and division evaluate to , , , and , respectively; negation evaluates to ; and the exponential, logarithm, square-root, sine, and cosine constructors evaluate to the corresponding real functions applied to . The normal-density constructor evaluates to . All operations are total: division by zero, including , evaluates to ; the real logarithm evaluates to at and to the logarithm of the absolute value at a nonzero negative input; and the real square root evaluates to on every nonpositive input. Thus every syntactically valid expression has a real value at every real input, with no domain-error condition.
For every , boundary is the real number . No signs or nonzeroness conditions are imposed: gives the identically zero function of , gives the constant value , negative is allowed and produces exponential growth in magnitude as increases when , negative is allowed, and negative, zero, and positive values of are all within the definition. The declaration defines this formula but does not itself assert any decay property.
For every , normalDensity is the real number . It is defined for every real ; the declaration introduces only this function and does not itself assert that it integrates to one or attach a probability measure to it.
Given arbitrary and an arbitrary total function , SolvesAbelEquation(b₀,c,f) means that for every real , if , then, writing and , the equality holds, where the integral is the real Lebesgue interval integral with its positive orientation from to . Equivalently, the boundary terms appearing in the equation are exactly and . For , the implication is automatically true and imposes no condition; no equation at is required. For , , while at the integration endpoint one has ; total real division makes both and equal to there, so the pointwise integrand at that endpoint is , and in any case a single endpoint has Lebesgue measure zero. No continuity, measurability, integrability, sign, normalization, or support condition is assumed directly of . Under the total convention used for the integral, a non-integrable weighted integrand has integral ; since the right-hand side is strictly positive for every , any satisfying the predicate must nevertheless make the displayed weighted integrand integrable on for every positive . The predicate is insensitive to on and to changes of on Lebesgue-null subsets of . It includes the degenerate parameter cases and , permits negative and negative , and asserts neither the existence nor the uniqueness of such an .
import Mathlib
open MeasureTheory
namespace RybinAI2026.P13
/-- A fixed, auditable language of elementary real expressions. -/
inductive ElementaryExpr where
| variable
| rational : ℚ → ElementaryExpr
| constant : ℝ → ElementaryExpr
| add : ElementaryExpr → ElementaryExpr → ElementaryExpr
| sub : ElementaryExpr → ElementaryExpr → ElementaryExpr
| mul : ElementaryExpr → ElementaryExpr → ElementaryExpr
| div : ElementaryExpr → ElementaryExpr → ElementaryExpr
| neg : ElementaryExpr → ElementaryExpr
| exp : ElementaryExpr → ElementaryExpr
| log : ElementaryExpr → ElementaryExpr
| sqrt : ElementaryExpr → ElementaryExpr
| sin : ElementaryExpr → ElementaryExpr
| cos : ElementaryExpr → ElementaryExpr
| normalDensity : ElementaryExpr → ElementaryExpr
/-- Semantics of the elementary-expression language. -/
noncomputable def ElementaryExpr.eval : ElementaryExpr → ℝ → ℝ
| .variable, t => t
| .rational q, _ => q
| .constant x, _ => x
| .add a b, t => a.eval t + b.eval t
| .sub a b, t => a.eval t - b.eval t
| .mul a b, t => a.eval t * b.eval t
| .div a b, t => a.eval t / b.eval t
| .neg a, t => -a.eval t
| .exp a, t => Real.exp (a.eval t)
| .log a, t => Real.log (a.eval t)
| .sqrt a, t => Real.sqrt (a.eval t)
| .sin a, t => Real.sin (a.eval t)
| .cos a, t => Real.cos (a.eval t)
| .normalDensity a, t => (2 * Real.pi) ^ (-(1 : ℝ) / 2) *
Real.exp (-(a.eval t) ^ 2 / 2)
/-- The exponentially decaying boundary. -/
noncomputable def boundary (b₀ c t : ℝ) : ℝ := b₀ * Real.exp (-c * t)
/-- The standard normal density. -/
noncomputable def normalDensity (x : ℝ) : ℝ :=
(2 * Real.pi) ^ (-(1 : ℝ) / 2) * Real.exp (-x ^ 2 / 2)
/-- Durbin's generalized Abel equation for an exponential boundary. -/
def SolvesAbelEquation (b₀ c : ℝ) (f : ℝ → ℝ) : Prop :=
∀ t : ℝ, 0 < t →
∫ s in (0 : ℝ)..t,
(1 / Real.sqrt (t - s)) *
normalDensity ((boundary b₀ c t - boundary b₀ c s) / Real.sqrt (t - s)) * f s =
(1 / Real.sqrt t) * normalDensity (boundary b₀ c t / Real.sqrt t)
end RybinAI2026.P13
Read-back
What the Lean code literally says, in plain math · gpt-5.6-sol
ElementaryExpr is a parameter-free inductive type whose elements are finite expression trees. It has a nullary variable expression; a rational-constant expression for every ; a real-constant expression for every ; binary addition, subtraction, multiplication, and division expressions formed from any two expressions; and unary negation, exponential, logarithm, square-root, sine, cosine, and normal-density expressions formed from any one expression. There are no side conditions on any constructor and no equations identifying distinct constructor trees.
For every elementary expression and every real input , ElementaryExpr.eval assigns a real value, denoted , recursively as follows: the variable evaluates to ; a rational constant evaluates to its image in , independently of ; a real constant evaluates to , independently of ; addition, subtraction, multiplication, and division evaluate to , , , and , respectively; negation evaluates to ; and the exponential, logarithm, square-root, sine, and cosine constructors evaluate to the corresponding real functions applied to . The normal-density constructor evaluates to . All operations are total: division by zero, including , evaluates to ; the real logarithm evaluates to at and to the logarithm of the absolute value at a nonzero negative input; and the real square root evaluates to on every nonpositive input. Thus every syntactically valid expression has a real value at every real input, with no domain-error condition.
For every , boundary is the real number . No signs or nonzeroness conditions are imposed: gives the identically zero function of , gives the constant value , negative is allowed and produces exponential growth in magnitude as increases when , negative is allowed, and negative, zero, and positive values of are all within the definition. The declaration defines this formula but does not itself assert any decay property.
For every , normalDensity is the real number . It is defined for every real ; the declaration introduces only this function and does not itself assert that it integrates to one or attach a probability measure to it.
Given arbitrary and an arbitrary total function , SolvesAbelEquation(b₀,c,f) means that for every real , if , then, writing and , the equality holds, where the integral is the real Lebesgue interval integral with its positive orientation from to . Equivalently, the boundary terms appearing in the equation are exactly and . For , the implication is automatically true and imposes no condition; no equation at is required. For , , while at the integration endpoint one has ; total real division makes both and equal to there, so the pointwise integrand at that endpoint is , and in any case a single endpoint has Lebesgue measure zero. No continuity, measurability, integrability, sign, normalization, or support condition is assumed directly of . Under the total convention used for the integral, a non-integrable weighted integrand has integral ; since the right-hand side is strictly positive for every , any satisfying the predicate must nevertheless make the displayed weighted integrand integrable on for every positive . The predicate is insensitive to on and to changes of on Lebesgue-null subsets of . It includes the degenerate parameter cases and , permits negative and negative , and asserts neither the existence nor the uniqueness of such an .
Confirmed by the mission captain (proposal self-audit).