Бинарный поиск: элемент, равный среднему геометрическому элементов массива

const
  N = 100;
 
type
  arr = array[1..N] of real;
 
var
  a: arr;
  i: byte;
  p, q: real;
  e: real;
 
procedure fillarr(var b: arr; min, max: integer; var sr_geom: real);
var
  k: byte;
begin
  sr_geom := 1;
  for k := 1 to N do 
  begin
    b[k] := random * (max - min) + min;
    sr_geom := sr_geom * b[k];
  end;
  sr_geom := power(sr_geom, 1 / n);
end;
 
procedure search(var c: arr; elem: real);
var
  m, i, j: integer;
begin
  m := N div 2;
  i := 1;
  j := N;
  while (c[m] <> elem) and (i <= j) do 
  begin
    if elem > c[m] then 
      i := m + 1
    else 
      j := m - 1;
    m := (i + j) div 2;
  end;
  if i > j then 
    writeln('Нет такого числа в массиве')
  else 
    writeln('Есть такое число в массиве');
end;
 
procedure sort(c: arr; elem: real);
var
  j, m, id: byte;
  k: real;
begin
  m := N;
  while m > 1 do 
  begin
    id := 1;
    for j := 2 to m do
      if c[j] > c[id] then 
        id := j;
    k := c[id];
    c[id] := c[m];
    c[m] := k;
    m := m - 1;
  end;
  search(c, elem);
end;
 
begin
  randomize;
  fillarr(a, 1, n, e);
  writeln('Среднее геометрическое: ', e);
  sort(a, e);
  for i := 1 to N do 
    write(a[i]:16:4);
  writeln;
end.