pub fn divergence(field: &[Expr], vars: &[&str]) -> Result<Expr, VecCalcError>
Divergence ∇·F = Σ ∂Fᵢ/∂xᵢ. Requires field.len() == vars.len().
∇·F = Σ ∂Fᵢ/∂xᵢ
field.len() == vars.len()