I have written the following code in R:
alpha <- rpois(n=15,lambda=10) data.X <- runif(15,2,3) #Defining the function jacobian_func <- function(X, a