We have added Skellam, Benford and SkewNormal distributions. We have added FitResult classes specific to distributions, with members that make it easier to access the parameters of that distribution. We have addede Gegenbauer polynomials and Chebyshev polynomials of the second kind. We have added Complex Airy functions. We have made improvements to the performance and accuracy of the Gamma family of functions (Gamma, Beta, Pochhammer) and added a LogPochhammer function. The Permutation class can now produce a PermutationMatrix, which can be used to permute rows and columns of vectors and matrices.