I try to integrate this :
integrate(integrate(integrate(2*sin(z)*cos(atan((2*cos(y)-0.5+x)/(2*sin(y)))),y,0,pi/2),x,0,1),z,0,pi/2);
Wolfram find the solution but I would like to control the accuracy. I try with tplquad but there is some error.
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)
The errors are:
Blockquote File "", line 3, in File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 526, in tplquad return dblquad(_infunc2,a,b,gfun,hfun,(func,qfun,rfun,args),epsabs=epsabs,epsrel=epsrel) File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 461, in dblquad return quad(_infunc,a,b,(func,gfun,hfun,args),epsabs=epsabs,epsrel=epsrel) File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 281, in quad retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points) File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 345, in _quad return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit) File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 406, in _infunc a = gfun(x)
Have you an idea where the error could come from ?