Фактический параметр должен быть переменной

program lalala;
const
     n=5;
Type
    mas = array[1..n,1..n] of real;
    vek = array[1..n] of real;
const
      a:mas =( (1,2,0,0,0), (3,4,5,0,0), (0,6,7,8,0), (0,0,9,10,11), (0,0,0,12,13) );
      b:vek = (6,13,19,25,31);
Var
   x,k,m,t,p,q:vek;
   tmp,ev:real;
 
 
Procedure InitMas(Var a:mas; Var b:vek);
Var i,j:integer;
begin
    write('b(',i,')=');
    readln(b[i]);
 end;
 
 
Procedure Prog(Var a:mas; Var k,m,t:vek);
Var i,j:integer;
Begin
for i:=1 to n do
    begin
    if i=1 then
       k[i]:=0
    else
        k[i]:=a[i,i-1];
        m[i]:=-a[i,i];
    if i=n then
        t[i]:=0
    else
        t[i]:=a[i,i+1];
    end;
 end;
 
Procedure Vec(Var a:mas; Var q,p,k,m,x:vek);
Var i,j:integer;
Begin
p[1]:=t[1]/m[1];
q[1]:=-b[1]/m[1];
for i:=2 to n do
    begin
    p[i]:=-t[i]/(k[i]*p[i-1]-m[i]);
    q[i]:=(b[i]-k[i]*q[i-1])/(k[i]*p[i-1]-m[i]);
    end;
x[n]:=(b[n]-k[n]*q[n-1])/(k[n]*p[n-1]-m[n]);
for i:=n-1 downto 1 do
   x[i]:=p[i]*x[i+1]+q[i];
end;
 
Procedure PrintB(Var a:mas; Var x:vek);
Var i,j:integer;
begin
for i:=1 to n do
    writeln('x(',i,')=',x[i]);
end;
 
Procedure Evklid(Var a:mas);
Var i,j:integer;
Var ev:real;
begin
ev:=0;
 for i:=0 to n do
        begin
            for j:=0 to n do
             begin
             ev:=ev + sqrt(sqr(a[i,j]));
             end;
        end;
         writeln('ev=',ev);
 end;
 
Procedure Proverka(Var a:mas);
Var i,j:integer;
f1:boolean;
Begin
     for i:=2 to n do
    Begin
        for j:=0 to  i-1 do
        Begin
            if ((a[i,j] <> 0) or (a[j,i] <> 0)) then
            Begin
            f1:= false;
            end;
        end;
    end;
 
if (f1) then
    Begin
        writeln('Matrix \n\n');
        Prog(a,k,m,t);
     end
     else
      Begin
      writeln('Error');
     end;
 end;
 
begin
     InitMas(a,b);
     Prog(a,k,m,t);
     Vec(a,p,q,k,b,x);
     PrintB(a,x);
     Evklid(a);
 readln();
end.

begin
     InitMas(a,b);          //Вот тут, но я подозреваю что будет еще
     Prog(a,k,m,t);        //Тут
     Vec(a,p,q,k,b,x);    //Тут
     PrintB(a,x);           //Тут 
     Evklid(a);             //А вот тут пишет  Ошибка 0: - выход за границы диапазона изменения диапазона 1..5, строка 74
 readln();