Для 3-их 4-мерных векторов a, b, c найти их попарное скалярное произведение

 uses crt;
Type k=array[1..4] of integer;
 
function f(ar:k;br:k):integer;
         var s,i:integer;
begin
     i:=1;
     while (i>=1) and (i<=4) do
       i:=i+1;
       S:=0;
       S:=S+ar[i]*br[i];
       f:=S;
end;
var i,j,z:integer; sp1,sp2,sp3:real;
    a:k; c:k; b:k;
    begin
    clrscr;
    randomize;
    writeln('vektor a');
    for i:=1 to 4 do begin
      a[i]:=random(10);
      write(a[i]:4);
    end;
    writeln;
      writeln('vector b');
    for i:=1 to 4 do begin
      b[i]:=random(21);
      write(b[i]:4);
    end;
    writeln;
      writeln('vektor c');
    for i:=1 to 4 do begin
      c[i]:=random(12);
      write(c[i]:4);
    end;
    writeln;
    sp1:=f(a,b);
    sp2:=f(b,c);
    sp3:=f(a,c);
writeln(' a*b sp1= ' ,sp1:0:5);
writeln(' b*c sp2= ' ,sp2:0:5);
writeln(' a*c sp3= ' ,sp3:0:5);
readln;

function f(ar, br:k):integer;

while (i>=1) and (i<=4) do
       i:=i+1;

S:=0;

a, b, c:k;

uses
    Crt;
const
    n = 4;
type
    TMass = Array[1..n] Of Integer;
function Scalar(const ar, br : TMass) : LongInt;
var
    i : Byte;
    Result : LongInt;
begin
    Result := 0;
    for i := 1 to n do
        Result := Result + ar[i] * br[i];
    Scalar := Result;
end;
procedure InitMass(var arr : TMass; const c : Char);
var
    i : Byte;
begin
    WriteLn;
    WriteLn('vector ', c, ': ');
    for i := 1 to n do
    begin
        arr[i] := Random(19) - 9;
        Write(arr[i]:3);
    end;
    WriteLn;
end;
 
var
    a, b, c : TMass;
begin
    InitMass(a, 'a'); InitMass(b, 'b'); InitMass(c, 'c');
    WriteLn('a*b: ', Scalar(a, b));
    WriteLn('a*c: ', Scalar(a, c));
    WriteLn('b*c: ', Scalar(b, c));
end.