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