где тут ошибка? уже голову сломала( код вроде простой

Program Chmet;
{
   Используем 2 метода, погрешность вычисляем
   через разницу между найденным значением и значением, полученным по формуле        Ньютона-Лейбница:
Определенный интеграл = F(b)-F(a), где f(x)=F'}
 
Uses Crt;
 
  function fx(k:real):real;
  begin
  fx:=1/(1+k*k);
  end;
 
const e=0.001; { точность. Константа }
 
label Z1,Z2; { раздел меток }
 
var   a,b:real;  { наши пределы }
      n,i:integer; { начальное кол-во отрезков и счетчик }
      d:real;     { дельта-икс - "длинна отрезка" }
      S:real;    { сумма площадей - т.е. наш определенный интеграл }
      metod:char;{ выбираем метод }
 
begin
clrscr;
writeln ('Введите пределы интегрирования  A,B и начальное количество отрезков N:');
write('a=');
readln(a);
write('b=');
readln(b);
write('n=');
readln(n);
writeln;
writeln('Выберите метод интегрирования:');
writeln('1. Центральных прямоугольников');
writeln('2. Трапеции');
 
repeat
metod:=readkey
until metod in ['1','2'];
writeln;
writeln('Вы выбрали метод ',metod);
goto Z1;
Z2:
n:=round(n*1.1); { округление }
S:=0;
d:=(b-a)/n;
Z1:
if metod='1' then
   begin
   for i:=1 to n-1 do
   S:=S+fx(a+d/2+d*(i-1))*d;
   end
else
   begin
   for i:=1 to n-1 do
   S:=S+(fx(a+d*(i-1))+fx(a+d*i))*d/2;
   end;
 
if abs(arctan(b)-arctan(a)-S)>e then goto Z2;
writeln('По методу ',metod,' S= ',S:1:10);
writeln('По формуле Ньютона-Лейбница S=',arctan(b)-arctan(a):1:10);
readkey
end.