using System;

using Accord.Statistics.Distributions.Univariate;

public class Program

{

    public static void Main()

    {

        // Create a new Poisson distribution with 

        var dist = new PoissonDistribution(lambda: 2.0);

        // Cumulative distribution functions

        double cdf1 = dist.DistributionFunction(k: 2); // 0.676676416183064

        double cdf2 = dist.DistributionFunction(k: 4); // 0.947346982656289

        double cdf3 = dist.DistributionFunction(k: 7); // 0.998903281032141     

        // Quantile function

        int icdf1 = dist.InverseDistributionFunction(p: cdf1); // 2

        int icdf2 = dist.InverseDistributionFunction(p: cdf2); // 4

        int icdf3 = dist.InverseDistributionFunction(p: cdf3); // 7

        Console.WriteLine("2 = " + icdf1);

        Console.WriteLine("4 = " + icdf2);

        Console.WriteLine("7 = " + icdf3);

    }

}