Используя оператор цикла с постусловием, вычислить сумму ряда с точностью (эпсилон)=10^(-3) — Pascal(Паскаль)

uses crt;

var
  n: integer;
  e, t, s: real;

begin
  clrscr;
  e := 0.001;
  s := 1 / 3;
  t := 1 / 3;
  n := 1;
  repeat
    n := n + 1;
    t := -t * 2 / (2 * n * (2 * n + 1));
    s := s + t;
  until abs(t) < e;
  write('S=', s:0:4, ' n=', n); { n-для проверки, можно не выводить }
  readln

end.