为什么 scipy 中的 mathieu 函数不是周期性的？

Question

我发现 scipy.special 中的 Mathieu 函数不是周期性的。这是代码和情节：

import scipy.special as sp
import matplotlib.pyplot as plt
import numpy as np
x = np.linspace(0,10*np.pi,100)
plt.plot(x,sp.mathieu_cem(3,3,x)[0])
plt.show()

但是在 Mathematica 中，它显然是周期性的：

Plot[MathieuC[MathieuCharacteristicA[3, 3], 3, x], {x, 0, 10 Pi}]

这两个功能不一样吗？由于 scipy 的文档没有太多细节，我不知道如何获取 ce.

的定期版本

对了，gsl里有similar functions，和scipy有什么关系？他们是同一个人吗？

编辑

转换为度后，我发现 scipy 中的循环比 mathematica 中的要多，有什么解释吗？

Answer 1

scipy.special.mathieu_cem 需要 x 度数。对于上升到 31.4 度的 Python 代码和上升到 31.4 弧度的 Mathematica 代码，此输出看起来是正确的。