Exemplo n.º 1
0
    private static void hermite_polynomial_test04()

    //****************************************************************************80
    //
    //  Purpose:
    //
    //    HERMITE_POLYNOMIAL_TEST04 tests H_POLYNOMIAL_ZEROS.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    19 October 2014
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        int degree;

        double[] hz;
        string   title;

        double[] z;

        Console.WriteLine("");
        Console.WriteLine("HERMITE_POLYNOMIAL_TEST04:");
        Console.WriteLine("  H_POLYNOMIAL_ZEROS computes the zeros of H(n,x)");
        Console.WriteLine("  Check by calling H_POLYNOMIAL there.");

        for (degree = 1; degree <= 5; degree++)
        {
            z     = Hermite.h_polynomial_zeros(degree);
            title = "  Computed zeros for H(" + degree + ",z):";
            typeMethods.r8vec_print(degree, z, title);

            hz    = Hermite.h_polynomial_value(degree, degree, z);
            title = "  Evaluate H(" + degree + ",z):";
            typeMethods.r8vec_print(degree, hz, title, aIndex: +degree * degree);
        }
    }