public void CanEstimateParameters(double mu, double lambda)
{
var original = new InverseGaussian(mu, lambda, new Numerics.Random.MersenneTwister(100));
var estimated = InverseGaussian.Estimate(original.Samples().Take(1000000));
AssertHelpers.AlmostEqualRelative(mu, estimated.Mu, 1);
AssertHelpers.AlmostEqualRelative(lambda, estimated.Lambda, 1);
}
