Module integrals 

Numeric line and surface integrals over parametric curves/surfaces given as closures. Derivatives of the parametrization are taken by central finite difference, so the caller supplies only the position map. Trapezoidal quadrature.

Functions

line_integral_scalar

Scalar line integral ∫_C f ds = ∫ f(r(t)) |r'(t)| dt.

line_integral_work

Work / vector line integral ∫_C F·dr = ∫ F(r(t))·r'(t) dt.

surface_flux

Flux of a 3-D field through a parametric surface r(u,v): ∫∫ F·(r_u × r_v) du dv. The orientation is that of r_u × r_v.