“…Subsequently, we utilize the Optimization.solve function from the Optimation.jl package (a unified optimization package) (Dixit and Rackauckas, 2023). This function employs reverse-mode automatic differentiation through Zygote.jl, allowing efficient computation of gradients and sensitivities required for optimizing the parameters.…”