public void ProbabilityMassFunctionTest()
{
MultivariateDiscreteDistribution target = new MultinomialDistribution(5, 0.25, 0.25, 0.25, 0.25);
int[] observation = { 1, 1, 1, 2 };
double expected = 0.05859375;
double actual = target.ProbabilityMassFunction(observation);
Assert.AreEqual(expected, actual, 1e-6);
}