Проверить гипотезу Эйлера: "Простое число вида 4n+1 является суммой двух квадратов"

function prost(x:integer):boolean;
var  y:integer;
     f:boolean;
begin
x:=abs(x);
if x<2 then f:=false{0,1 не простое}
else if x=2 then f:=true{2 простое}
else if x mod 2=0 then f:=false{четные больше 2 не простые}
else
 begin
  f:=true;
  y:=3;
  while(y*y<=x)and f do
  if x mod y=0 then f:=false
  else inc(y,2);
 end;
prost:=f
end;
procedure sqr2(n:integer;var a,b:integer;var f:boolean);
begin
f:=false;
a:=1;
while(a<=trunc(sqrt(n)))and not f do
 begin
  b:=a;
  while(b<=trunc(sqrt(n)))and not f do
  if (a*a+b*b=n) then f:=true
  else inc(b);
  if not f then inc(a);
 end;
end;
var n,a,b,k:integer;
    f,t:boolean;
begin
n:=5;
t:=true;
k:=0;
while n<=2000 do
 begin
  if prost(n) then
   begin
    sqr2(n,a,b,f);
    if not f then t:=false
    else
     begin
      writeln(n,'=',(n-1)div 4,'*4+1=',a,'^2+',b,'^2');
      inc(k);
      if k mod 22=0 then
       begin
        write('Нажмите Enter');
        readln
       end;
     end;
   end;
  inc(n,4);
 end;
if not t then write('В диапазоне 1..2000 эта гипотеза не верна');
end.