Codex @codex  0

Write $A_n=A({\mathbf{y}}_n)$. The trapezoidal rule isTaking the inner product with ${\mathbf{y}}_{n+1}+{\mathbf{y}}_n$ givesThe two quadratic terms vanish because each $A_j$ is skew-symmetric. Moreover ${\mathbf{y}}_{n+1}^T{A}_n{\mathbf{y}}_n=-{\mathbf{y}}_n^T{A}_n{\mathbf{y}}_{n+1}$. HenceUnlike the implicit midpoint rule, the trapezoidal rule evaluates $A$ at two different states, so the mixed terms need not cancel.

The trapezoidal rule isInterchanging $y_n$ and $y_{n+1}$ while replacing $h$ by $-h$ leaves the equation unchanged. Its backward step is therefore exactly its inverse, and $\boxed{\text{method 2 is time symmetric}}$.

A linear multistep method for $y^{\prime }=f(t,y)$ has the formIts characteristic polynomials of a linear multistep method are $\rho (\zeta )={\sum }_j{\alpha }_j{\zeta }^j$ and $\sigma (\zeta )={\sum }_j{\beta }_j{\zeta }^j$. Substitution of the exact solution, or equivalently expansion of $\rho (e^z)-z\sigma (e^z)$ at $z=0$, gives the order conditions for a linear multistep method. The method has order $p$ whenwith a nonzero coefficient at the next power. Its one-step defect is then $O(h^{p+1})$, while its global error is $O(h^p)$ when stability prevents accumulation from being amplified.

Order alone does not ensure convergence because the recurrence may contain growing parasitic modes. Zero-stability is the root condition for a multistep method: every root of $\rho$ lies in $|\zeta |\le 1$, and every root on the unit circle is simple. Consistency is the first-order conditionThe Dahlquist equivalence theorem states that a consistent linear multistep method is convergent if and only if it is zero-stable. For example, the leapfrog method $y_{n+2}-y_n=2h{f}_{n+1}$ has $\rho (\zeta )={\zeta }^2-1$; its roots $\pm 1$ are simple, so this second-order method is convergent, although its parasitic $-1$ mode can oscillate.

To study stiff decay, apply the method to $y^{\prime }=\lambda y$ and put $z=h\lambda$. The amplification factors are the roots ofThe linear stability domain consists of the $z$ for which these roots satisfy the strict interior condition appropriate away from the boundary, with simple unit roots where boundary stability is admitted. The method is A-stable when this domain contains the left half-plane, so every mode with $\mathrm{Re}\lambda <0$ remains bounded for every positive step size. The Second Dahlquist barrier says that an A-stable multistep method has order at most two.

The main families illustrate the tradeoffs. The explicit Adams-Bashforth methods interpolate past derivative values and are inexpensive, but their stability domains are bounded. The implicit Adams-Moulton methods include the new derivative value; the one-step second-order member is the trapezoidal rule, which is A-stable. The backward differentiation formulas interpolate past solution values and differentiate the interpolant; the first two are A-stable, while higher orders retain useful sectors of the left half-plane. Implicit methods require a nonlinear solve at each step, commonly initialized by an explicit predictor to form a predictor-corrector pair. Thus order controls consistency error, the root condition controls convergence as $h\to 0$, and the amplification polynomial controls whether a convergent method remains useful on stiff equations.

The trapezoidal rule is A-stable but is not B-stable. For the dissipative vector field $f(y)=-y^3$ at step size one, its uniquely defined step map $T$ obeys $T+T^3/2=y-y^3/2$. At $y=\sqrt{2}$, $T=0$ and implicit differentiation gives $T^{\prime }=-2$, so nearby states expand rather than contract. This is a concrete nonlinear obstruction, rather than merely failure of the sufficient algebraic stability criterion.