procedure shrink(var s: scalar);
var
  r: hreal;
begin { shrink }
  rsum(s.l,s.u,r);
  irdiv(2,r);
  s.l:=r;
  s.u:=r;
end { shrink };

