procedure sdiff(var x1,x2,x3: scalar);
begin { sdiff }
  underflow:=false;
  r_diff(down,x1.l,x2.u,x3.l);
  r_diff(up,x1.u,x2.l,x3.u);
  if underflow then snormalize(x3);
end { sdiff };

