public MultinomialDistribution(int numberOfTrials, params double[] probabilities) : base(probabilities.Length) { N = numberOfTrials; this.probabilities = probabilities; nfac = Special.Factorial(numberOfTrials); }