DeprecationWarning for scipy: quadpack

[hyanwong]

I'm getting this in the test suite. I guess this is a variational_gamma think @nspope, and would be relatively easy to correct @nspope?

tests/test_functions.py: 1260 warnings
  /Users/yan/miniforge3/lib/python3.10/site-packages/scipy/integrate/_quadpack_py.py:613: DeprecationWarning: Conversion of an array with ndim > 0 to a scalar is deprecated, and will error in future. Ensure you extract a single element from your array before performing this operation. (Deprecated NumPy 1.25.)
    return _quadpack._qagie(func, bound, infbounds, args, full_output,

[nspope]

I think this is a scipy issue (not fixable by us) -- but I'll check.

[nspope]

Oh, actually, I was planning on removing the quadrature here anyway, in favor of exact solutions (the quadrature can live in a separate validation notebook; it's annoyingly slow and unstable to include in the tests). The warning should go away when I do this.