※上の続き
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
int main()
{
double i,imax,n;
double x,y,pi[10],ave,sum;
int a;
for(a=0; a<10; a++)
{
srand((unsigned)time(NULL)%32768);
n=0.0;
imax=10000.0; //乱数発生回数
for(i=0; i<=imax; i++)
{
x=rand()/(RAND_MAX+1.0); //範囲(0,1)の乱数発生
y=rand()/(RAND_MAX+1.0);
if((x*x+y*y)<=1.0) //0<=x,y<=1の範囲にある中心原点１の４分円内に入ればカウント
{
n+=1.0;
}
}
pi[a]=n/imax*4.0; //半径１の円の面積、すなわち円周率
}
for(a=0, sum=0; a<10; a++) //結果の合計
sum += pi[a];
ave=sum/10.0; //結果の平均
printf("π＝%f\n",ave);
return 0;
}