procedure lvho(imin: integer; var s,sho: scalar; var v: vector);
var
  i: integer;
  sf,sp,st: scalar;
begin { lvho }
  squot(sone,s,sf);
  sf.l.si:=0;
  sabs0(sho,sp);
  for i:=imin to vdim-1 do begin
     ssymm(sp,v[i]);
     sprod(sf,sp,st);
     sp:=st;
     isdiv(i+1,sp);
  end;
  sabs0(sp,v[vdim]);
end { lvho };