Computes the Anderson–Darling distribution function at x, with parameter n, using Marsaglia’s et al. algorithm. A port based on the implementation in SSJ.