% Script roundoff3.m

% this version of sqrt(1 + x) - 1  has no subtractions

f2 = @(x) x / (sqrt(1 + x) + 1);

num2hex(f2(0.123456789012345E-1))
num2hex(f2(0.123456789012345E-5))
num2hex(f2(0.123456789012345E-10))
num2hex(f2(0.123456789012345E-12))
num2hex(f2(0.123456789012345E-15))