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) + ""); } }