#include <cmath>
#include <iostream>

int main()
{
  const auto pi = std::acos(-1.0);
  
  std::cout << "# k=0,n=0 k=0,n=0.5 k=1,n=0.5\n";
  
  for (int i = 0; i < 100; ++i)
  {
    double phi = 0.01*i*pi/2;
    std::cout << phi 
      << '\t' << std::ellint_3(0.0, 0.0, phi)
      << '\t' << std::ellint_3(0, 0.5, phi)
      << '\t' << std::ellint_3(1, 0.5, phi)
      << '\n';
  }
}