§ 8. Вычисление значения числа π методом Монте-Карло

8.5. Создание компьютерной расчетной модели (этап 3в)

В системе PascalABC.NET создадим программу montekarlo. В ней объявим переменные n и k типа integer для хранения числа точек-песчинок на базовом квадрате и на круге соответственно, а также переменные s, x и y типа real для хранения значений площади круга и координат точек-песчинок соответственно.

В основном разделе программы, задаем начальные значения и организуем цикл (пример 8.5). Далее подсчитываем результат, выводим на экран результат и точное значение числа :

s := 4 * k / n;
writeln('Результат pi = ',s);
writeln('Точно pi = ',pi);

Пример 8.5. Задаем начальные значения:

n := 1000;
k := 0;

В цикле for с начальным значением переменной цикла 1 и конечным значением n присваиваем случайные значения координатам очередной точки:

x := 2 * random();
y := 2 * random();

С помощью условного оператора if организуем подсчет числа k точек, которые попали в круг:

if sqr(x-1)+sqr(y-1)<=1 then k:=k+1;

Осталось подсчитать площадь круга по основной формуле метода.