Hi Team, gh-16907 (https://github.com/scipy/scipy/pull/16907) changes not only the limits of integration for the computation of expectation values in the vonmises distribution class (vonmises::expect), it also introduces a complex exponential (exp(i x)) as a default function to estimate the mean value. These changes close gh-16655 (https://github.com/scipy/scipy/issues/16655). Associated with this proposed change to the vonmises class, is a change to the integrate.quad method to support the integration of complex valued functions, see gh-17051 (https://github.com/scipy/scipy/pull/17051). The change there required the addition of a new argument, complex_func (which defaults to False), to allow users to opt into complex integration. If you are interested, please take a look at these PRs and let me know if you have feedback that I can incorporate. Thanks and Cheers, Chris