Multidimentional integrals in julia: speed and automatic differentiation

Specific Domains Numerics
4

My working problem is more complicated than that.
However, trying Zygote is an interesting idea. Could maybe sketch an example?
Let’s say, the parameter is a vector of dim2. Thanks in advance!