internal static double probability(double x, double mean, double stdev, bool cumulative)
{
var normalDistribution = new MathNet.Numerics.Distributions.Normal(mean, stdev);
return(cumulative ? normalDistribution.CumulativeDistribution(x) : normalDistribution.Density(x));
}