fix: ϕ₅

Author: ymtoo

src/AlphaStableDistributions.jl

@@ -172,22 +172,22 @@ const ϕ₃ = [
     1.908  1.908  1.908  1.908  1.908
 ]
 const ϕ₅ = [
-    0.0 0.0    0.0    0.0    0.0
-    0.0 -0.017 -0.032 -0.049 -0.064
-    0.0 -0.030 -0.061 -0.092 -0.123
-    0.0 -0.043 -0.088 -0.132 -0.179
-    0.0 -0.056 -0.111 -0.170 -0.232
-    0.0 -0.066 -0.134 -0.206 -0.283
-    0.0 -0.075 -0.154 -0.241 -0.335
-    0.0 -0.084 -0.173 -0.276 -0.390
-    0.0 -0.090 -0.192 -0.310 -0.447
-    0.0 -0.095 -0.208 -0.346 -0.508
-    0.0 -0.098 -0.223 -0.383 -0.576
-    0.0 -0.099 -0.237 -0.424 -0.652
-    0.0 -0.096 -0.250 -0.469 -0.742
-    0.0 -0.089 -0.262 -0.520 -0.853
-    0.0 -0.078 -0.272 -0.581 -0.997
-    0.0 -0.061 -0.279 -0.659 -1.198
+    0.0  -0.061  -0.279  -0.659  -1.198
+    0.0  -0.078  -0.272  -0.581  -0.997
+    0.0  -0.089  -0.262  -0.52   -0.853
+    0.0  -0.096  -0.25   -0.469  -0.742
+    0.0  -0.099  -0.237  -0.424  -0.652
+    0.0  -0.098  -0.223  -0.383  -0.576
+    0.0  -0.095  -0.208  -0.346  -0.508
+    0.0  -0.09   -0.192  -0.31   -0.447
+    0.0  -0.084  -0.173  -0.276  -0.39
+    0.0  -0.075  -0.154  -0.241  -0.335
+    0.0  -0.066  -0.134  -0.206  -0.283
+    0.0  -0.056  -0.111  -0.17   -0.232
+    0.0  -0.043  -0.088  -0.132  -0.179
+    0.0  -0.03   -0.061  -0.092  -0.123
+    0.0  -0.017  -0.032  -0.049  -0.064
+    0.0   0.0     0.0     0.0     0.0
 ]
 
 """
@@ -221,7 +221,7 @@ function Distributions.fit(::Type{<:AlphaStable}, x::AbstractArray{T}, alg=Quick
     c = (p[6]-p[2]) / itp₃(α, abs(β))
     itp₄ = interpolate((_α, _β), ϕ₅, Gridded(Linear()))
     ζ = p[4] + c * sign(β) * itp₄(α, abs(β))
-    if abs(α - 1.0) < 0.1
+    if abs(α - 1.0) < 0.05
         δ = ζ
     else
         δ = ζ - β * c * tan(π*α/2)