Sympy でいくつかの単純な確率を計算しようとしていますが、数値エラーが発生し続けます。mpmath で任意の精度を使用するようにこのコードを変更するにはどうすればよいですか?
from sympy import *
from sympy.stats import Binomial, density, P
from mpmath import mp
mp.dps = 1000
n = 128
X = Binomial('X', n, 1/2, succ=1, fail=0)
for t in range (n+1):
print(t, P(X<t))