title: How to use tplquad?
body: <p>I try to integrate this :</p>
  
  <p>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);</p>
  
  <p>Wolfram find the solution but I would like to control the accuracy. I try with tplquad but there is some error. </p>
  
  <pre><code>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)
  </code></pre>
  
  <p>The errors are:</p>
  
  <blockquote>
    <p>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)</p>
  </blockquote>
  
  <p>Have you an idea where the error could come from ?</p>