Module integrate 

Numerical integration (trapezoidal, Simpson’s rule). Numerical Integration — deterministic quadrature routines.

Determinism Contract

All routines use KahanAccumulatorF64 from cjc_repro for floating-point summation, guaranteeing bit-identical results across runs for the same inputs.

Functions

• [trapezoid] — composite trapezoidal rule from paired (x, y) arrays
• [simpson] — composite Simpson’s 1/3 rule from paired (x, y) arrays
• [cumtrapz] — cumulative trapezoidal integration, returns array of running integrals

Functions

cumtrapz

Cumulative trapezoidal integration.

simpson

Composite Simpson’s 1/3 rule for numerical integration.

trapezoid

Composite trapezoidal rule for numerical integration.