Сравнить скорость сходимости разных разложений Pi - Pascal ABC

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

Сравнить скорость сходимости (число слагаемых для достижения заданной точности ε) следующих разложений π:

Решение задачи: «Сравнить скорость сходимости разных разложений Pi»

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

var n1,n2,n3,z:integer;
    s,s1,t,e:real;
begin
repeat
write('Введите точность 0.00001<=e<1 e=');
readln(e);
until(e>=0.00001)and(e<1);
n1:=0;
t:=4;
s:=t;
z:=1;
while abs(s-pi)>e do
 begin
  n1:=n1+1;
  z:=-z;
  t:=4/(2*n1+1);
  s:=s+z*t;
 end;
writeln('s1=',s:0:6,' n1=',n1+1);
n2:=1;
t:=4/24;
s:=3+t;
z:=1;
while abs(s-pi)>e do
 begin
  n2:=n2+1;
  z:=-z;
  t:=4/(2*n2)/(2*n2+1)/(2*n2+2);
  s:=s+z*t;
 end;
writeln('s2=',s:0:6,' n2=',n2+1);
n3:=1;
t:=6;
s1:=t;
s:=sqrt(s1);
while abs(s-pi)>e do
 begin
  n3:=n3+1;
  t:=6/n3/n3;
  s1:=s1+t;
  s:=sqrt(s1);
 end;
writeln('s3=',s:0:6,' n3=',n3+1);
end.

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

Переменные: n1, n2, n3, z (целые числа), s, s1, t, e (действительные числа).
Пользователю предлагается ввести точность e (с помощью цикла repeat-until).
Вычисляются начальные значения s, t, z.
Цикл while (абсолютное значение разности s и pi больше e) выполняет итерационный процесс.
В каждой итерации значение z меняется на противоположное, значение t вычисляется как 4/(2n1+1), а значение s обновляется как сумма zt.
Значения n2, t, s вычисляются заново после каждой итерации цикла while.
Значение t вычисляется как 4/(2n2)/(2n2+1)/(2*n2+2).
Значения n3, t, s1 вычисляются заново после каждой итерации цикла while.
Значение t вычисляется как 6/n3/n3.
Значение s1 обновляется как сумма t.
Значение s вычисляется как квадратный корень из s1.
Цикл while завершается, когда абсолютное значение разности s и pi меньше или равно e.
Выводятся значения s, n1, n2, n3.
Отображается результат s3, полученный с использованием метода Монте-Карло для вычисления числа Пи.
Отображается результат s2, полученный с использованием метода Монте-Карло для вычисления числа Пи.
Отображается результат s1, полученный с использованием метода Монте-Карло для вычисления числа Пи.
Отображается результат s, полученный с использованием метода Монте-Карло для вычисления числа Пи.
Конец программы.