Циклический сдвиг элементов квадратной матрицы вправо (оформить в виде процедуры) - Pascal

const nmax=20;
type matr=array[1..nmax,1..nmax] of integer;
procedure Sdvig(var a:matr;n,pr,k:integer);
var p,q,i,j,j1,s,x,y:integer;
begin
p:=n-2*(pr-1);
s:=4*p-4;
k:=k mod s;
y:=pr-1;
for i:=1 to k do
 begin
  x:=a[pr+1,pr];
  for j:=1 to s-1 do
  if j<p-1 then a[j+y+1,pr]:=a[j+y+2,pr]{первый столбец}
  else if j in [p-1..2*p-3] then a[p+y,j+y-p+2]:=a[p+y,j+y-p+3]{последняя строка}
  else if j in [2*p-2..3*p-4] then
  a[3*p-j+y-2,p+y]:=a[3*p-j+y-3,p+y]{последний столбец}
  else if j in [3*p-3..s-1] then a[pr,s-j+y+1]:=a[pr,s-j+y];{первая строка}
  a[pr,pr]:=x;
 end;
end;
var a:matr;
    n,m,i,j,r:integer;
begin
repeat
write('Размер матрицы от 2 до ',nmax,' n=');
readln(n);
until n in [2..nmax];
writeln('Матрица A:');
for i:=1 to n do
 begin
   for j:=1 to n do
    begin
     a[i,j]:=(i-1)*n+j;
     write(a[i,j]:4);
    end;
   writeln;
 end;
writeln;
write('На сколько позиций сдвинуть элнменты в слоях r=');
readln(r);
i:=1;
repeat
Sdvig(a,n,i,r);
i:=i+1;
until i>n div 2;
writeln('Сдвиг по часовой стрелке на ',r,' позиций:');
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:4);
  writeln;
 end;
end.