我正在尝试在 C 中编写 pi 的准蒙特卡罗近似。我还不太精通它,并且正在尝试翻译我的基于 python 的技能,所以我可能只是忽略了一些东西。结果我一直得到0,我不知道为什么。我应该如何解决这个问题?另外,我在最后两个调用中收到错误,printf
说它们是 typedouble *
而不是double
。无论如何它都会编译,这有关系吗?
#include <stdio.h>
/*
Tristen Wentling
montepithon.c
October 31, 2013
*/
int main(void)
{
float i,j,x;
float count=0,counter=0;
printf("Please enter the desired grid division size (n=?)");
scanf("%f", &x);
float y=x*x;
if(x==0){
printf("goodbye");
}
else if(x!=0){
for(i=0;i<=x;i++){
for(j=0;j<=x;j++){
float check=((i*i)*(1/y))+((j*j)*(1/y));
/*printf("%f\n", check);*/
if(check<=1){
count+=1;
}
else{
counter+=1;
}
}
}
}
else{
printf("error");
}
float prsum=count/y;
float ptsum=(1-counter)*(1/y);
double pirprox=4*prsum;
double pitprox=4*ptsum;
printf("%f\n", &pirprox);
printf("%f\n", &pitprox);
getchar();
}