using System;

using MathNet.Numerics.Distributions;

​

public class Program

{

    public static void Main()

    {

        // Parameters for the binomial distribution

        int n = 100;        // Number of trials

        double p = 1.0 / 4; // Probability of success

​

        // Mean and standard deviation for the normal approximation

        double mean = n * p;

        double stddev = Math.Sqrt(n * p * (1 - p));

​

        // Create a Normal distribution object for approximation

        var normalDist = new Normal(mean, stddev);

​

        // Generate a random sample from the normal distribution

        double sample = normalDist.Sample(); // No continuity correction applied, instead rounding to the nearest integer

​

        // Round the result to the nearest integer

        int roundedSample = (int)Math.Round(sample);

​

        // Enforce the range between 0 and n

        roundedSample = Math.Max(0, Math.Min(n, roundedSample));

​

        // Print the result

        Console.WriteLine("Random sample from Normal approximation (rounded and bounded): {0}", roundedSample);

    }

}

​