procedure num_heap_sort(var num:real_array;n:integer);
var i,j,r,l:integer;
    temp:real;
    not_ready:boolean;
begin
  if n>0 then
  begin
    l:=n div 2+1;
    r:=n;
    repeat
      if l>1 then l:=l-1
      else if r>1 then begin
                         temp:=num[l];
                         num[l]:=num[r];
                         num[r]:=temp;
                         r:=r-1;
                       end;
      i:=l;
      j:=2*i;
      temp:=num[i];
      not_ready:=(j<=r);
      while not_ready do
        begin
          if j<r then if num[j]<num[j+1] then j:=j+1;
          not_ready:=false;
          if j<=r then not_ready:=(temp<num[j]);
          if not_ready then begin
                              num[i]:=num[j];
                              i:=j;
                              j:=2*i;
                            end;
        end;
      num[i]:=temp;
    until r=1;
  end
  else writeln('<num_heap_sort>: n <  1 ');
end;


                                    