Ejemplo n.º 1
0
    public static void triangular_cdf_test()
    //****************************************************************************80
    //
    //  Purpose:
    //
    //    TRIANGULAR_CDF_TEST tests TRIANGULAR_CDF.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    31 March 2016
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        int i;
        int seed = 123456789;

        Console.WriteLine("");
        Console.WriteLine("TRIANGULAR_CDF_TEST");
        Console.WriteLine("  TRIANGULAR_CDF evaluates the Triangular CDF;");
        Console.WriteLine("  TRIANGULAR_CDF_INV inverts the Triangular CDF.");
        Console.WriteLine("  TRIANGULAR_PDF evaluates the Triangular PDF;");

        double a = 1.0;
        double b = 10.0;

        Console.WriteLine("");
        Console.WriteLine("  PDF parameter A =      " + a + "");
        Console.WriteLine("  PDF parameter B =      " + b + "");

        if (!Triangular.triangular_check(a, b))
        {
            Console.WriteLine("");
            Console.WriteLine("TRIANGULAR_CDF_TEST - Fatal error!");
            Console.WriteLine("  The parameters are not legal.");
            return;
        }

        Console.WriteLine("");
        Console.WriteLine("       X            PDF           CDF            CDF_INV");
        Console.WriteLine("");

        for (i = 1; i <= 10; i++)
        {
            double x   = Triangular.triangular_sample(a, b, ref seed);
            double pdf = Triangular.triangular_pdf(x, a, b);
            double cdf = Triangular.triangular_cdf(x, a, b);
            double x2  = Triangular.triangular_cdf_inv(cdf, a, b);

            Console.WriteLine("  "
                              + x.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "  "
                              + pdf.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "  "
                              + cdf.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "  "
                              + x2.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "");
        }
    }