円周率が出ない。。。

[山本]

#include <stdio.h>
#include <stdlib.h>

int main()
{
double pi;
printf("pi=%f\n",pi);
}

double calc_pi()
{double x,y,i,pi;
int count;
for(i=0;i<=1000;i++){
x=(double)rand() / (RAND_MAX + 1);
y=(double)rand() / (RAND_MAX + 1);
if((x*x+y*y)<1.0)
{count+=1.0;}
}
pi=count/1000*4.0;
return (int)pi;
}

これで円周率が出てくれないのはなぜ？？

[fatens]

>これで円周率が出てくれないのはなぜ？？

・calc_pi 関数を呼びだしていない
・calc_pi 関数のプロトタイプ宣言がない
・countを0で初期化していない
・pi=count/1000*4.0　count/1000はint型同士の割り算なので結果もint型になってしまう
・return (int)pi int型にキャストしたのでは、小数点以下の値が帰らない

など

[初級者]

プログラムは、思ったとおりには動きません。
書いたとおりに動きます。

[パコネコ]

実際にうごいてるのは
/*
double pi;
printf("pi=%f\n",pi);
*/
これだけ…

＞int count;
int（整数）でいいんですか？
どちらかというと「i」のほうが整数っぽい

[初級者]

円の内側にある点の「個数」を数えるための変数ですから、
「当然」整数です。

[パコネコ]

そうですか…
計算が楽になるかな～と思ったのですが…
それに1.0ずつ足していたので整数型ではないと考えたんですが、
ダメみたいですね。すいませんでした。

[たいちう]

> それに1.0ずつ足していたので整数型ではないと考えたんですが、
> ダメみたいですね。すいませんでした。

違和感を感じる着眼点は正しいですが、結論が逆でしたね。

1.0ずつ足しているのに整数型、ではなく、
整数型なのに1.0ずつ足している、です。

[QERA]

#include <stdio.h>
#include <stdlib.h>

double calc_pi(void);

int main()
{
double pi;
pi = calc_pi();
printf("pi=%lf\n",pi);
}

double calc_pi()
{double x,y,pi,count=0;
int i;
for(i=0;i<=1000;i++){
x=(double)rand() / (RAND_MAX + 1);
y=(double)rand() / (RAND_MAX + 1);
if((x*x+y*y)<1.0)
{count+=1.0;}
}
pi=count/1000*4.0;
return pi;
}

こうですかね？

[初級者]

1001回の試行で得られたcountの値を1000で割っているため、
やや厳密性に欠けるでしょう。
まあ、モンテカルロ法自体がそうなんですけどね。

[フリオ]

　
ごめんなさい。勘違いです。
countはdoubleですね。