pub fn gradient(f: &Expr, vars: &[&str]) -> Vec<Expr>
∇f = [∂f/∂x₁, …, ∂f/∂xₙ].
∇f = [∂f/∂x₁, …, ∂f/∂xₙ]