Определить площадь круга с использованием метода "Монте-Карло" - Pascal ABC

Формулировка задачи:

определить площадь круга с использованием метода "Монте- Карло"

Решение задачи: «Определить площадь круга с использованием метода "Монте-Карло"»

Листинг программы

uses crt;
var r,x,y,m,n,i:integer;
    s:real;
begin
randomize;
write('R=');
readln(r);
m:=0;
n:=10000;
for i:=1 to n do
 begin
  x:=random(2*r+1)-r;
  y:=random(2*r+1)-r;
  if x*x+y*y<r*r then m:=m+1
 end;
s:=4*r*r*m/n;
writeln('Вычисленная площадь=',s:0:2);
write('Аналитически=',pi*r*r:0:2);
end.

Объяснение кода листинга программы

Используется библиотека CRT для работы с консолью.
Задаются переменные: r - радиус круга, x, y - координаты точек, m - количество точек внутри круга, n - общее количество точек, i - номер текущей точки.
Переменная s инициализируется значением площади круга с радиусом r.
Выводится запрос на ввод радиуса круга.
Производится инициализация переменных m и n.
В цикле происходит генерация случайных точек внутри квадрата со стороной 2r+1 и вычисление их количества, удовлетворяющих условию xx+yy<r*r.
Значение переменной m используется для вычисления площади круга по методу Монте-Карло.
Выводится результат вычисления площади круга.
Выводится аналитическое значение площади круга (πrr).
Программа завершается.