21.15109/ARP/BIPBSASajti, Szilárd MihálySzilárd MihálySajti0000-0002-8748-8242Wigner Research Centre for PhysicsARP2024Computer and Information ScienceEngineeringMathematical SciencesPhysicsα-stable distributioncorrelation functionnumerical instabilityWynn's ε-methodSajti, Szilárd MihálySzilárd MihálySajtiWigner Research Centre for Physics2024-09-122024-09-20419255245170469634696875446147699525971910258159976523474473071871447912738847654357474539940761972325218627003189119986378009224372574278869540512176606image/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngimage/pngtext/x-c++srctext/x-chdrapplication/pdftext/plainapplication/pdftext/x-c++srctext/x-chdrtext/plaintext/plaintext/plaintext/plaintext/plaintext/plaintext/plaintext/x-chdr1.0CC BY-NC 4.0The α-stable distributions appear in a broad field. We needed them in X-ray and neutron off-specular scattering, where their characteristic functions are frequently used as correlation functions. In practice, we need the Fourier-transforms of the 2-dimensional (and usually cylindrically symmetric) correlation functions, i.e. the corresponding probability density functions. At first, we used implementations, which we have found sorrily to be numerically unstable. Therefore, we have written a new C++ code using the GNU MPFR (C library for multiple-precision floating-point computations) and boost multiprecision libraries. We present the ‘test’ results generated by this new code here. Additionally, we very shortly summarize the theoretical bases. For further details, we refer the reader to the given references, and the source code and its documentation also are provided.