一种计算π的方法

更新时间：2022-09-15 16:52:57

假设点可以均匀的扔到正方形中。如果一共扔N次，其中M次扔进内切圆内，则可以推导出计算π的公式：

一种计算π的方法

随机计算π

大体思路：获得随机数（x,y），判断x*x + y*y是否在单位圆内，可以得到M的个数。

程序：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 1000
int main()
{
    int i, M=0;
    double x, y, pi;
    srand((int)time(NULL));
    for (i=0; i<N; i++)
    {
        x = rand() / (double)RAND_MAX;
        y = rand() / (double)RAND_MAX;
        if (x*x + y*y < 1)
            M++;
    }
    pi = (4.0 * M) / N;
    printf("pi:%f\n", pi);
    return 0;
}

注意：有关随机数的获取可以参考：http://www.cnblogs.com/kaituorensheng/archive/2013/03/05/2944008.html

一组结果：一种计算π的方法

画出网格计算π

大体思路：

程序：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define n 10000000
int main()
{
    int i;
    double sum = 0;
    for (i=0; i<n; i++)
        sum += sqrt((double)n*n - (double)i*i);
    printf("%f\n", 4.0*sum /n/n);
}

一组结果：

一种计算π的方法

本文转自jihite博客园博客，原文链接：http://www.cnblogs.com/kaituorensheng/p/3145090.html，如需转载请自行联系原作者

上一篇 : ：找出一堆数中个数超过一半的数下一篇 : Apache+Subversion+TortoiseSVN

一种计算π的方法

相关阅读

推荐文章