Найти угол между векторами

const m = 100;
 
type vec = array [1..m] of double;
 
procedure read_v(s: string; var x: vec; k: integer);
var i: integer;
begin
  writeln('Введите координаты вектора ', s);
  for i := 1 to k do
    begin
      write(s, '[', i, '] = ');
      readln(x[i])
    end
end;
 
procedure prn_v(s: string; var x: vec; k: integer);
var i: integer;
begin
  writeln('Вектор ', s);
  for i := 1 to k do write(' ', x[i]:19);
  writeln
end;
 
function p(var x, y: vec; k, mode: integer): double;
var i: integer;
    z: double;
begin
  z := 0;
  for i := 1 to k do z := z + x[i] * y[i];
  if mode = 1 then z := sqrt(z);
  p := z
end;
 
function arccos(a: double): double;
begin
  if a = -1
    then arccos := pi
    else if a = 0
      then arccos := pi / 2
      else if a = 1
        then arccos := 0
        else arccos := arctan(sqrt(1 - a * a) / a)
end;
 
var n: integer;
    a, b: vec;
 
begin
  repeat
    write('Введтите размерность пространства ( не менее 1 ): ');
    readln(n)
  until n > 1;
  read_v('A', a, n);
  read_v('B', b, n);
  writeln('Исходные векторы:');
  prn_v('A', a, n);
  prn_v('B', b, n);
  writeln;
  write('Угол между векторами: ',
         arccos(p(a, b, n, 0) / p(a, a, n, 1) / p(b, b, n, 1)):0:11,
         ' радиан.');
  readln
end.

         arccos(p(a, b, n, 0) / p(a, a, n, 1) / p(b, b, n, 1)) / pi * 180:0:11,