final double mu = 4.5;
final double sigma = 0.86;
final ProbabilityDistribution<Double> p1 = new NormalDistribution(mu, sigma, new MersenneTwister64(MersenneTwister.DEFAULT_SEED));
final double[] x = new double[n];
for (int i = 0; i < n; i++) {
x[i] = p1.nextRandom();
}
final NormalDistribution p2 = (NormalDistribution) CALCULATOR.evaluate(x);
assertEquals(p2.getMean(), mu, 2.5e-2);
assertEquals(p2.getStandardDeviation(), sigma, 2.5e-2);
}