「Buffonの針」記事へのコメント

そりゃ出ませんよ

[nekopon (1483)]

まずは以下のように直してみてください

harisaki=sin((double)rand()/RAND_MAX*M_PI/2)*0.5+tyuuten;

針先の角度(Θ)が一様分布するのであって、針先と中点のある軸方向に対する距離(0.5sin Θ)が一様分布するわけじゃあないわけですんで。

残念ながら昨日提出しました

[yakusouX5 (8222)]

一応、中点を決めてそこで回転した針のX座標をX軸に下ろしたつもりなんですが
更に、ここには書いていませんが課題の条件で角度やπを使うことは禁止となっているので、そこで困ってしまいました

--
うすっぺらいコメントがあらわれた！ ▼

Re:残念ながら昨日提出しました

[nekopon (1483)]

気がつくのが遅くてごめんなさい。m(__)m それとなかなか厳しい条件ですね。(Buffon の針なら、角度が一様分布するのが本質的な条件になるので) 難しそうなので、しばらく考えさせてください。

別の方法で乱数からpiを近似する方法ならご紹介できます:

double monte_carlo(int repeats)
{
  int count, i;
  double x, y;

  for (i = 0, count = 0; i < repeats; i++) {
    x = rand() / RAND_MAX;
    y = rand() / RAND_MAX;
    if (x * x + y * y < 1.0)
       count++;
  }
  return 4.0 * count / repeats;
}

Re:残念ながら昨日提出しました

[yakusouX5 (8222)]

いえいえそんなに悩まないでください。先生も分かっていないでその条件にしている可能性も否定できませんので。ちょっとボケが入っていそうです
夏休み明けを一週間間違えて初回の講義が休講になったり。しかも旅行に行っていたとか

別の方法の乱数からπを求める方法は前回の課題で既に…

--
うすっぺらいコメントがあらわれた！ ▼