In the program we set m = 1/2. Thus k = 1/2 and from sinak we find a = /4. // elliptic.cpp #include #include using namespace std; int main(void) { double pi 3.14159265; double m 0.5%; } double a0 = 1.0, b0 = sqrt (1.0-m), c0 = sqrt (m); double eps = 0.0001; double al, b1, c1%; while (fabs (c0) > eps) { a1 = (a+b0)/2.0; b1 = sqrt(a0b0); c1 = (a0b0)/2.0; a0a1; b0b1; c0 = c1; } double Kalpha = pi/(2.0*a1); cout < < "Kalpha return 0; = " < < Kalpha;