JuliaMath / HypergeometricFunctions.jl

A Julia package for calculating hypergeometric functions
Other
68 stars 8 forks source link

Implement Kummer's transformations #46

Closed rashidrafeek closed 1 year ago

rashidrafeek commented 2 years ago
julia> HypergeometricFunctions.M(0.5, 1.5, -1000)
NaN

This works correctly in scipy:

from scipy.special import hyp1f1

hyp1f1(0.5, 1.5, -1000)
0.02802495608198965

This used to not work in scipy too (https://github.com/scipy/scipy/issues/2282) but was fixed in https://github.com/scipy/scipy/pull/3022.

MikaelSlevinsky commented 1 year ago

Thanks! It seems to me that the issue is solved by implementing https://dlmf.nist.gov/13.7#E2, equivalently https://dlmf.nist.gov/13.2#E41 and https://dlmf.nist.gov/13.2#E42