使用四邊形（scipy.integrate）和小數

Question

我需要高度精確地集成一個函數，例如：

from decimal import * 
from scipy.integrate import quad
getcontext().prec = 30
a = Decimal('1.1') + Decimal('1.2') 
print type(a)

def f(X):
    return quad(lambda x: X, 0, 1)[0]

b = f(a)
print "{0:.30f}".format(b), type(b)

我得到

<class 'decimal.Decimal'>
2.300000000000000266453525910038 <type 'float'>

我想使用Quad作為我的集成工具，但不幸的是，似乎將輸出轉換為浮點型。 我希望積分的輸出具有更高的精度。 我該如何實現？

Answer 1

我將看看使用mpmath庫。 它支持高精度集成（ http://mpmath.org/doc/0.19/calculus/integration.html ）。