Now you'll need to use math.factorial(int(Decimal)) or write your own factorial function.
The storyline:
Apropos of recent threads experimenting with high precision (arbitrary precision) numbers as a way to promote interest in "pure math" applications using Python, I discovered this morning that my
Ramanujan Convergence script on repl.it was broken all of a sudden (after working before), and set out to discover why.
Answer: in Python 3.8, to which
repl.it has newly upgraded, math.factorial no longer accepts Decimal type numbers, even if they're integral (integers). You can read the discussion here:
https://bugs.python.org/issue33083
Here's the line of code that broke, and was fixed with explicit castings to int.
term = (fact(int(c1*i))*(c2 + c3*i))/(pow(fact(int(i)),4)*pow(c4,4*i))
By the time of the final division, the numerator and denominator have been coerced back into Decimals and so will divide with the expected precision (set by the context): Where:
```
c1 = Decimal(4)
c2 = Decimal(1103)
c3 = Decimal(26390)
c4 = Decimal(396)
c5 = Decimal(9801)
```
For more context:
(I'm updating it now...)
Kirby