#include <cmath>
#include <iostream>

int main()
{
  std::cout << "# k n=-3 n=-0.5 n=0 n=0.5\n";
  
  for (int i = 0; i < 100; ++i)
  {
    double k = 0.01*i;
    std::cout << k 
      << '\t' << std::comp_ellint_3(k, -3.0) 
//       << '\t' << std::comp_ellint_3(k, -1.0)  // inf
      << '\t' << std::comp_ellint_3(k, -0.5) 
      << '\t' << std::comp_ellint_3(k, 0.0) 
      << '\t' << std::comp_ellint_3(k, 0.5)
//      << '\t' << std::comp_ellint_3(k, 1.0) // domain error
//      << '\t' << std::comp_ellint_3(k, 3.0) // domain error
      << '\n';
  }
}