예제 #1
0
    public static void elliptic_inc_pik_values_test()
    //****************************************************************************80
    //
    //  Purpose:
    //
    //    ELLIPTIC_INC_PIK_VALUES_TEST tests ELLIPTIC_INC_PIK_VALUES.
    //
    //  Licensing:
    //
    //    This code is distributed under the GNU LGPL license.
    //
    //  Modified:
    //
    //    24 June 2018
    //
    //  Author:
    //
    //    John Burkardt
    //
    {
        double k   = 0;
        double n   = 0;
        double phi = 0;
        double pik = 0;

        Console.WriteLine("");
        Console.WriteLine("ELLIPTIC_INC_PIK_VALUES_TEST:");
        Console.WriteLine("  ELLIPTIC_INC_PIK_VALUES stores values of");
        Console.WriteLine("  the incomplete elliptic integral of the third");
        Console.WriteLine("  kind, with parameters PHI, N and K.");
        Console.WriteLine("");
        Console.WriteLine("    PHI           N             K            Pi(PHI,N,K)");
        Console.WriteLine("");
        int n_data = 0;

        while (true)
        {
            Elliptic.elliptic_inc_pik_values(ref n_data, ref phi, ref n, ref k, ref pik);
            if (n_data == 0)
            {
                break;
            }
            Console.WriteLine("  " + phi.ToString("0.########").PadLeft(12)
                              + "  " + n.ToString("0.########").PadLeft(12)
                              + "  " + k.ToString("0.########").PadLeft(12)
                              + "  " + pik.ToString("0.################").PadLeft(24) + "");
        }
    }