示例#1
0
    public static void jacobi_poly_values_test()
    //****************************************************************************80
    //
    //  Purpose:
    //
    //    JACOBI_POLY_VALUES_TEST tests JACOBI_POLY_VALUES.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    19 April 2012
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        double a  = 0;
        double b  = 0;
        double fx = 0;
        int    n  = 0;
        double x  = 0;

        Console.WriteLine("");
        Console.WriteLine("JACOBI_POLY_VALUES_TEST:");
        Console.WriteLine("  JACOBI_POLY_VALUES returns values of");
        Console.WriteLine("  the Jacobi polynomial.");
        Console.WriteLine("");
        Console.WriteLine("       N         A         B      X       J(N,A,B)(X)");
        Console.WriteLine("");
        int n_data = 0;

        for (;;)
        {
            Jacobi.jacobi_poly_values(ref n_data, ref n, ref a, ref b, ref x, ref fx);
            if (n_data == 0)
            {
                break;
            }

            Console.WriteLine("  "
                              + n.ToString(CultureInfo.InvariantCulture).PadLeft(6) + "  "
                              + a.ToString(CultureInfo.InvariantCulture).PadLeft(8) + "  "
                              + b.ToString(CultureInfo.InvariantCulture).PadLeft(8) + "  "
                              + x.ToString(CultureInfo.InvariantCulture).PadLeft(10) + "  "
                              + fx.ToString("0.################").PadLeft(24) + "");
        }
    }