search

haskell / math-functions

Special mathematical functions
http://hackage.haskell.org/package/math-functions
BSD 2-Clause "Simplified" License
40 stars 28 forks source link

Evaluate x**z · exp(-x) / Γ(z+1) precisely #42

Closed Shimuuar closed 4 years ago

Shimuuar commented 8 years ago

At the moment it's evaluated as exp (log x * z - x - logGamma (z+1)) which gives only 12 significant digits for large 100-1000 z.

Shimuuar commented 8 years ago

When this issue fixed test invIGammaIsInverse should be updated too,