Найти натуральное число в интервале от 1 до 1000 с максимальной суммой делителей — Pascal(Паскаль)

uses crt;
var i,j,s,mx,imx:integer;
begin
clrscr;
mx:=1;imx:=1;
for i:=1 to 1000 do
 begin
  s:=0;
  for j:=1 to i do
  if i mod j=0 then s:=s+j;
  if s>mx then
   begin
    mx:=s;
    imx:=i;
   end;
 end;
write('Максимальную сумму делителей=',s,' имеет число ',imx);
readln
end.