В Паскале вывести все числа, имеющие один делитель, не считая 1 и

Var a,b,i,j,k:integer;
begin
write('a = ');
readln(a);
write('b = ');
readln(b);
for i:=a to b do
 begin
 k:=0;
 for j:=2 to i div 2 do
  if i mod j=0 then k:=k+1;
 if k=1 then write(i,' ');
 end;
end.

Пример:
a = 1
b = 500
4 9 25 49 121 169 289 361 

PS. Разыскиваемые числа - это квадраты обычных чисел, потому можно и более эффективную программу написать.