Дано число N. Найти сумму чисел в интервале от 1 до N, взаимно простых с N — Pascal(Паскаль)

uses crt;
var n,i,j:integer;
    s:longint;
    f:boolean;
begin
clrscr;
write('n=');
readln(n);
writeln('Числа, взаимно простые с ',n);
s:=0;
for i:=1 to n do
 begin
  f:=true;
  j:=2;
  while f and(j<=i) do
   begin
    if(i mod j=0)and(n mod j=0) then f:=false
    else j:=j+1;
   end;
  if f then
   begin
    write(i,' ');
    s:=s+i;
   end;
 end;
writeln;
write('Их сумма=',s);
readln
end.