Определить наибольший угол между прямой и заданными плоскостями при помощи процедуры пользователя

type vector = array ['x'..'z'] of double;
 
procedure angle(a, n: vector; var alp: double);
var aan, aa, an: double;
    i: char;
begin
  aan := 0;
  aa := 0;
  an := 0;
  for i := 'x' to 'z' do
    begin
      aan := aan + a[i] * n[i];
      aa := aa + sqr(a[i]);
      an := an + sqr(a[i]);
    end;
  alp := aan / sqrt(aa) / sqrt(an)
end;
 
procedure read_q(s: string; var q: double);
begin
  write(s, ' = ');
  readln(q)
end;
 
procedure read_v(s, t: string; var v: vector);
var i: char;
begin
  repeat
    for i := 'x' to 'z' do read_q(s[byte(i) - byte('w')], v[i]);
    if (v['x'] = 0) and (v['y'] = 0) and (v['z'] = 0) then writeln('Это не уравнение ', t, ', повторите ввод')
  until (v['x'] <> 0) or (v['y'] <> 0) or (v['z'] <> 0)
end;
 
var i: integer;
    plane: array [1..2] of vector;
    line: vector;
    alpha: array [1..2] of double;
 
begin
  writeln('Прямая задана параметрически: x=x0+kt, y=y0+mt, z=z0+nt');
  read_v('kmn', 'прямой', line);
  writeln('Плоскости заданы в общей форме: Ax+By+Cz+D=0');
  for i := 1 to 2 do
    begin
      writeln('Введите коэффициенты для ', i, ' плоскости:');
      read_v('ABC', 'плоскости', plane[i])
    end;
  for i := 1 to 2 do angle(line, plane[i], alpha[i]);
  if alpha[1] > alpha[2]
    then i := 1
    else if alpha[1] < alpha[2]
      then i := 2
      else i := 3;
  if i = 3
    then write('Углы между прямой и каждой из плоскостей одинаковы и равны ', alpha[1]:0:5)
    else write('Наибольший угол: между прямой и плоскостью ', i, ', равен ', alpha[i]:0:5);
  write(' радиан');
  readln
end.