Результатами вычислений по формуле x^2+x+40 при 0<=x<=40 являются простые числа".Все результаты вывести на экран - Pascal(Паскаль)

program p;
uses crt;
var a,i,j:integer;
f:boolean;
begin
clrscr;
 f:=true;
for i:=0 to 40 do
 begin
   a:=(i*i)+i+40;
   for j:=2 to round(sqrt(a)) do{считать надо только до корня}
   if a mod j =0 then f:=false;
   writeln(i,'^2+',i,'+40=',a);
   if i=20 then readln;{задерка, а то все на экран не входит}
 end;
   if not f then writeln('  Утверждение неверно!')
   else writeln('  Утверждение верно!');
readln
end.