procedure ssqrt(var x1,x2: scalar);
begin { ssqrt }
  rsqrt(x1.l,x2.l);
  r_sqrt(up,x1.u,x2.u);
end { ssqrt };

