CassiusT 2019-04-25 10:01
浏览 293

用蒙特卡洛方法算浦丰投针问题,为什么算出来的结果pi是在分子上的?

public class randomPoints {
double x;
double y;
static int intersect = 0;
static int total = 0;

public double getx() {
    return this.x;
}

public double gety() {
    return this.y;
}

randomPoints() {
    x = (Math.random() );
    y = (Math.random() );
    randomPoints points = sndP(new randomPoints(x, y));
    double a = points.getx();
    double b = points.gety();
    if (b>1&&y>0) {
        intersect++;
    }
    total++;
}

randomPoints(double xy, double yx) {
    x = xy;
    y = yx;
}

public static randomPoints sndP(randomPoints p) {
    double a = Math.random() + p.x - 0.5;
    boolean c = Math.random() - 0.5 > 0;
    double b;
    if (c) {
        b = p.y + (double) Math.sqrt(0.25 - (a - p.x) * (a - p.x));
    } else {
        b = p.y - (double) Math.sqrt(0.25 - (a - p.x) * (a - p.x));
    }
    return new randomPoints(a, b);
}

public static void main(String[] args) {
    for (int i = 0; i <= 999999999; i++) {
        randomPoints points = new randomPoints();

// System.out.print(points.x+" "+points.y+"\n" );
}
float c = (float)intersect / (float)total;
System.out.print(4*c);
}
}

  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 如何在scanpy上做差异基因和通路富集?
    • ¥20 关于#硬件工程#的问题,请各位专家解答!
    • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
    • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
    • ¥30 截图中的mathematics程序转换成matlab
    • ¥15 动力学代码报错,维度不匹配
    • ¥15 Power query添加列问题
    • ¥50 Kubernetes&Fission&Eleasticsearch
    • ¥15 報錯:Person is not mapped,如何解決?
    • ¥15 c++头文件不能识别CDialog