您好我正在尝试集成以下功能
使用∫ 0 1 ∫ 0 1 1 / B(-log x,-log y)dx DY
scipy.integrate.dblquad(lambda y,x: 1/ scipy.special.beta(-math.log(x),-math.log(y) ), 0, 1, lambda x: 0, lambda x: 1)
在Ipython 3.5 Spyder IDE中,其中B()是Beta函数。我得到6.81 ......准确度为10 -6 。但是在Linux上,当我用Python 2.7运行相同的代码时,我得到2.05的正确答案(用Mathematica验证)。我很迷惑。发生了什么?
修改
我在Computational Science SE发布了此问题,我建议将其发布在此处。我使用的是Windows 10,Python 3.5.2 | Anaconda 4.3.0(64位),scipy版本0.18.1,我得到了错误的答案。有人能够重现它们吗?