Реализация алгоритма "Решето Сундарама" для поиска простых чисел(Delphi)

[inviter]

Реализовываю код для вычисления простых чисел на заданном отрезке.

Program prj2;

{$APPTYPE CONSOLE}

uses SysUtils;

var
n, i, j, k, q: integer;
A: array of Integer;

Begin
  write('Please enter value of N=');
  readln(N);

  setlength(A,N);   {создание массива}
  for i:=0 to N do a[i]:=1; {и инициализация его всеми единицами}


Begin   {цикл прохода по всему массиву с 1}
 i:=1;
 j:=1;
 While (3*i+1)<n Do
 Begin                      {вычеркивание всех составных чисел}
   k:=i+j+2*i*j;
   inc(i);
   While ((k < n) and (j <= i)) Do
   begin
     inc(j);
     a[k]:=0;
   end;
 end;
end;

Begin    {вывод всех простых чисел}
  q:=0;
  i:=1;
  While i<n Do
  Begin
    inc(i);
    if a[i]=1 Then
      for i:=1 to n do
      begin
        write((2*i+1):4);
        inc(q);
        if  (q mod 5)=0 then
        writeln;
      end;
  end;
End;

A:=NIL;
readln;
End.

Разумеется, ничего не работает. Как мне кажется, проблема кроется в той части, где из массива исключаются составные числа. Подскажите в чём ошибка и как её исправить. Спасибо!