Задача Бюффона о бросании иглы

[Young_programmer]

Как решать задачу Бюффона с помощью программирования? Помогите, пожалуйста, с алгоритмом.

[type_Oleg]

Эту что ли? Бю
Что значит решить? Смоделировать бросание иглы?

Если это, то вот, разбирайтесь. На Delphi, без верхушки, там понятно что.

const
 HH=100.0;  // высота поля, куда бросают

implementation
uses Math;
{$R *.dfm}
procedure NeedleThrow(L:Double;var h1,h2:Double);
begin
 h1:=HH*Random;
 h2:=h1+L*Sin(2*PI*Random);
end;
function NumberOfIntersections(h1,h2,r,L:Double):Integer;
var s,nr,i:Integer;
    t:Double;
begin
 s:=0;
 nr:=Floor((HH+L)/r);
 if h2<h1 then
  begin
   t:=h1;
   h1:=h2;
   h2:=t;
  end;
 for i:=0 to nr do
  if (i*r>h1)and(i*r<h2) then Inc(s);
 Result:=s; 
end;
procedure TForm1.Button1Click(Sender: TObject);
var h1,h2,r,L:Double;
   N,k,i:Integer;
begin
 r:=10;          // расстоянии между линиями
 L:=9;           //  длина иглы
 N:=1000000; //  число бросков
 k:=0;           //  число пересечений
 for i:=1 to N do
  begin
   NeedleThrow(L,h1,h2);
   k:=k+NumberOfIntersections(h1,h2,r,L);
  end;
 Edit1.Text:=FloatToStr(k/N);
end;

initialization
 Randomize;
end.

[type_Oleg]

Нашел у себя ошибку. В функции NumberOfIntersections надо так:

for i:=-Round(L/r) to nr do

Это сказывается при длинной игле.

[Young_programmer]

А можно пошагово объяснить каким образом моделировали?

[type_Oleg]

Сообщение от Young_programmer

пошагово

NeedleThrow: Моделиуем ( генерируем ) координаты концов иглы. Горизонтальные координаты не имеют никакого значения, их не рассматриаем.
Random - возвр. случ. велич., равномерно распределенную 0...1.
h1 - координата одного конца иглы.
h2 - координата другого конца иглы. Кстати, может быть h1<h2, а может быть h1>h2.

NumberOfIntersections - подсчет количества пересечений. Их может быть и > 1, если длина иглы L > r (расст. между линиями).