ddevd_symbolic

Symbolic verification of the asymptotic claims in Appendix C of the DDEVD paper.

It reproduces:

• the leading order of int b_0(y) b_2(y) dy ~ -C_b * n^{-gamma},
• the leading order of int V_2(y) dy ~ C_V * n^{1 - gamma/2},
• the resulting stability bound m < n^{1 + gamma/2}.

Layout

ddevd_symbolic/
├── main.py                       # entry point — runs all three computations
└── ddevd_symbolic/
    ├── primitives.py             # f_X, alpha, beta, rho, z_1, z_2, Z (in t)
    ├── coefficients.py           # b_0, b_2, V_2 integrand expressions
    ├── watson.py                 # Watson's lemma + modified Watson + n-power extractor
    └── asymptotics.py            # driver applying Watson to each additive term

Reproduce

cd ddevd_symbolic
uv run python main.py            # or: python main.py

Output is the symbolic leading term and the corresponding power of n for each of the two integrals plus the m-bound.