Вычисление арктангенса числа -Pascal(Паскаль)

program c;

const e=0.001;
var a1,a0:real;
    i:integer;
begin a1:=1;i:=1;
      writeln(a1);
      repeat a0:=a1;
             a1:=arctan(a0)+1;
             writeln(a1);
             i:=i+1;

      until abs(a1-a0)<e;
      writeln(i);
   end.