add news

Author: heltonmc

NEWS.md

@@ -18,6 +18,7 @@ For bug fixes, performance enhancements, or fixes to unexported functions we wil
 
 ### Fixed
  - fix cutoff in `bessely` to not return error for integer orders and small arguments (#33)
+ - fix NaN return for small arguments (issue [#35]) in bessely (#40)
 
 # Version 0.1.0
 

src/bessely.jl

@@ -351,7 +351,7 @@ function bessely_power_series(v, x::T) where T
     out2 = zero(T)
     a = (x/2)^v
     # check for underflow and return limit for small arguments
-    iszero(a) && return -T(Inf)
+    iszero(a) && return (-T(Inf), a)
 
     b = inv(a)
     a /= gamma(v + one(T))