Из множества чисел [1..n] выделить подмножество простых чисел вида 4k+1 — Pascal(Паскаль)

uses crt;
function Prost(n:longint):boolean;
var i:longint;
begin
for i:=2 to round(sqrt(n)) do
if n mod i=0 then
Prost:=false;
end;
var n,i,k:longint;
begin
clrscr;
write('n=');readln(n);
k:=0;
for i:=1to n do
if Prost(i)and(i mod 4=1)then
 begin
  writeln(i);
  k:=k+1;
  if k mod 25=0 then readln
 end;
readln
end.