Compute Jacobian
jacobian(expr, vars)
'caracas expression'.
variables to take derivative with respect to
if (has_sympy()) {
x <- paste0("x", seq_len(3))
def_sym_vec(x)
y1 <- x1 + x2
y2 <- x1^2 + x3
y <- c(y1, y2)
jacobian(y, x)
u <- 2 + 4*x1^2
jacobian(u, x1)
}
#> c: [8⋅x₁]