Is there an R function for computing the Jacobean matrix (approx.) from a given vector input?

Question

I\'m curious if there\'s a way to write a code for the computation of a Jacobian matrix from a given vector as an input:

1. Using the midpoint method for approximating