Using LowLevelCallable for scipy.integrate.nquad

Question

I would like to evaluate the multi-dimensional intergral of the function f. Here f take an array of variables and return a float.

I try to u