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

    //****************************************************************************80
    //
    //  Purpose:
    //
    //    TEST0105 tests BESSEL_I0.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    26 August 2006
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        double fx = 0;
        double x  = 0;

        Console.WriteLine("");
        Console.WriteLine("BESSEL_I0_TEST:");
        Console.WriteLine("  BESSEL_I0 evaluates the Bessel function of the");
        Console.WriteLine("  first kind and order 0;");
        Console.WriteLine("");
        Console.WriteLine("      X       Exact F       BESSEL_I0(X)");
        Console.WriteLine("");

        int n_data = 0;

        for (;;)
        {
            Bessel.bessel_i0_values(ref n_data, ref x, ref fx);

            if (n_data == 0)
            {
                break;
            }

            double fx2 = Bessel.bessel_i0(x);

            Console.WriteLine("  "
                              + x.ToString(CultureInfo.InvariantCulture).PadLeft(8) + "  "
                              + fx.ToString(CultureInfo.InvariantCulture).PadLeft(16) + "  "
                              + fx2.ToString(CultureInfo.InvariantCulture).PadLeft(16) + "");
        }
    }
Exemplo n.º 2
0
    public static void bessel_i0_values_test()
    //****************************************************************************80
    //
    //  Purpose:
    //
    //    BESSEL_I0_VALUES_TEST tests BESSEL_I0_VALUES.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    07 February 2007
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        double fx = 0;
        double x  = 0;

        Console.WriteLine("");
        Console.WriteLine("BESSEL_I0_VALUES_TEST:");
        Console.WriteLine("  BESSEL_I0_VALUES stores values of ");
        Console.WriteLine("  the Bessel I0 function.");
        Console.WriteLine("");
        Console.WriteLine("      X         I0(X)");
        Console.WriteLine("");
        int n_data = 0;

        for (;;)
        {
            Bessel.bessel_i0_values(ref n_data, ref x, ref fx);
            if (n_data == 0)
            {
                break;
            }

            Console.WriteLine("  "
                              + x.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "  "
                              + fx.ToString(CultureInfo.InvariantCulture).PadLeft(12) + "");
        }
    }