我尝试整合这个:
积分(积分(积分(2 * sin(z)* cos(atan((2 * cos(y)-0.5 x)/(2 * sin(y)))),y,0,pi / 2),x ,0,1)中,z,0,π/ 2);
Wolfram找到解决方案,但我想控制准确性 . 我尝试使用tplquad,但有一些错误 .
def f(x,y,z):
return 2*sin(z)*cos(atan((2*cos(y)-0.5+x)/(2*sin(y))))
tplquad(f,0,1,0,pi/2,0,pi/2)
错误是:
Blockquote文件“”,第3行,在文件“/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py”,第526行,在tplquad中返回dblquad(_infunc2,a,b,gfun, hfun,(func,qfun,rfun,args),epsabs = epsabs,epsrel = epsrel)文件“/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py”,第461行,在dblquad中返回quad(_infunc,a,b,(func,gfun,hfun,args),epsabs = epsabs,epsrel = epsrel)文件“/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py”,第281行,在quad retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points)文件“/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py” ,第345行,在_quad中返回_quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit)文件“/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py” ,406行,_infunc a = gfun(x)
你知道错误可能来自哪里吗?