title: Calculation going to infinity
body: <p>I'm trying to calculate the gamma function of a number and it is going to infinity. When I use maple (for instance) it returns the correct answer (2.57e1133 (yes, it's huge)). I tried to use DECIMAL but 0 success. Is there a solution? Thanks in advance.</p>
  <p>the code</p>
  <pre><code>import scipy as sp                             
  from scipy.special import gamma
  from decimal import Decimal
  from scipy import special
  
  
  def teste(k):
      Bk2 = Decimal(gamma((1/(2*k))+(3/4)))
      return Bk2
  
  print(teste(0.001))
  </code></pre>
  <p>Result</p>
  <p>Infinity</p>