Найти произведение натуральных чисел от 1 до n

pascal abc;
var a,n,pr: longint;
begin
pr:=1;
writeln('введите число n ');
readln(n);
for a:=1 to n do begin
pr:=a*pr;
end;
writeln('произведение чисел от 1 до ',n,' = ',pr);
end.