Take the numerical derivative for the first derivative of an ARMA using the 2 point rule.
jacobian_arma(theta, p, q, tau)A vec that contains all the parameter estimates.
A int that indicates the number of AR coefficients
A int that indicates the number of MA coefficients.
A vec containing the scales e.g. \(2^{\tau}\)
A mat that returns the first numerical derivative of the ARMA process.