public static void test47()
//****************************************************************************80
//
// Purpose:
//
// TEST47 tests TORUS_SQUARE_5C2 and TORUS_SQUARE_14C.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 03 April 2008
//
// Author:
//
// John Burkardt
//
{
int i;
string name = "";
const double r1 = 1.0;
const double r2 = 0.125;
Console.WriteLine("");
Console.WriteLine("TEST47");
Console.WriteLine(" For integrals inside a torus with square cross-section:");
Console.WriteLine(" TORUS_SQUARE_5C2 approximates the integral;");
Console.WriteLine(" TORUS_SQUARE_14C approximates the integral.");
Console.WriteLine("");
Console.WriteLine(" Inner radius = " + r1 + "");
Console.WriteLine(" Outer radius = " + r2 + "");
Console.WriteLine(" Volume = " + Torus.torus_square_volume_3d(r1, r2) + "");
Console.WriteLine("");
Console.WriteLine(" F(X) 5C2 14C");
Console.WriteLine("");
int num = functions.function_3d_num();
for (i = 1; i <= num; i++)
{
int function_3d_index = i;
functions.function_3d_name(function_3d_index, ref name);
double result1 = Torus.torus_square_5c2(function_3d_index, functions.function_3d, r1, r2);
double result2 = Torus.torus_square_14c(function_3d_index, functions.function_3d, r1, r2);
Console.WriteLine(" " + name
+ " " + result1.ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + result2.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
