円周率を求めるC言語のプログラム

質問

int a=10000,b,c=8400,d,e,f[8401],g;main(){for(;b-c;)f[b++]=a/5;for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);

というのが本（「π魅惑の数」）に載っていたので
#include <stdio.h>

int main(void){
int a=10000,b,c=8400,d,e,f[8401],g;

for(;b-c;)f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
return 0;
}
としたのですが
円周率を表示してくれません。
（ｂ＝0にすると一応は出てくるのですが微妙に間違ってます（よく分からないですけどそのくらいの誤差のあるプログラムなのでしょうか？)）

回答No.1

liar_adan · ベストアンサー

おそろしく読みにくくてわけのわからないプログラムですね。
(けなしてるのではなくて、書いた人に感心しています)
「こんなんでπが出るのか！？」と思ったけれど、
出るようで驚きました。

私には計算内容は理解できないので、
プログラム上の事だけ話します。
下の方のプログラムで変数を宣言しているところで、
b=0, ... ,e=0, ..
とすると、上のプログラムと同じ数値を出します。

main()などの関数の中で宣言されている変数(自動変数)は、
値が不定です。つまり最初どういう値になっているかわかりません。

対して、関数の外側で宣言されている変数(静的変数)は、
0で初期化されることになっています。
この違いのため、バグが出ることがあります。
ここではeの初期値が0でなかったので動作が違ってしまったのでしょう。

お礼

>関数の外側で宣言されている変数(静的変数)は、
0で初期化されることになっています。

始めて知りました。

意味は全く分らないのですが
これはオイラーが見つけた級数を使っているらしいです。
　　　∞ 　　n
π＝２Σ　―――――
　　ｎ＝0 (2n+1)

回答ありがとうございました。