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