josef-pkt/scipy — comment on commit 2a34dcd780e9f089118b830af03e02f275a55192
File:  line 0
By josef-pkt:

The simplest would be to set it to (pi**2)/3 in the range (-5e-6, 5e-6) or something like this.

```
>>> stats.tukeylambda.stats(5e-6)[1] - np.pi**2 / 3
-4.6385342992394385e-05
>>> stats.tukeylambda.stats(-5e-6)[1] - np.pi**2 / 3
3.3550714780616886e-05
```

The difference looks monotonic up to this range and starts to bounce and go wrong closer to zero using this

```
>>> np.array([[k*10**-i, stats.tukeylambda.stats(k*10**-i)[1] - np.pi**2 / 3]
for i in range(1,8) for k in [0.25, 0.5,0.75,1][::-1]])

>>> np.array([[-k*10**-i, stats.tukeylambda.stats(-k*10**-i)[1] - np.pi**2 / 3]
for i in range(1,8) for k in [0.25, 0.5,0.75,1][::-1]])
```