楕円の円周

例えば、長半径0.5、短半径0.5、すなわち、直径1の真円として計算してみてください。そうすれば, 3.14158335634999245734 と出るはずです。つまり、円周率です。円周率の真値は,3.1415926535・・・だそうですので、精度としては、（長半径1の時）、 0.0001の位までだと分かります。

なぜ、1000箇所に区切ったのかというと、それ以上は、タイムアウトで計算できなかったなどの問題があったからです。もし、ご自分の計算機で計算されるなら、もっと精度の高い値が得られるかもしれません。

＿＿＿＿＿

以下に、プログラムを示します。精度=0.001の部分を改変し、もっと小さい値にして、かつ、2行目のn<=1000を、1÷精度の値になるように調整してください。

k=0;精度=0.001;縦横比率=短半径/長半径;
for(n=1;n<=1000;n=n+1)
{ j=(sqrt((精度^2)+(縦横比率^2*((sqrt(1-(精度*(n-1))^2)-sqrt(1-(精度*n)^2))^2))));
k=k+j;
答え=k*長半径*4;}
print(答え);

＿＿＿＿＿

あるいは、他の計算機なら、こう記述するかもしれません。

Σ (n=1,1/★) {√((★^2)+(▲^2*((√(1-(★*(n-1))^2)-√(1-(★*n)^2))^2)))}

★の部分が精度（例えば、0.001）で、▲が短半径/長半径の値です。ただし、この計算結果に、×長半径×4としたものが、その楕円の円周です。