I'm checking two different procedures that will eventually be used with regular floating point arithmetic. I need to estimate numerical error arising from floating point arithmetic.
Is there a way to force computation on floating point numbers? Model problem could be error estimation for:
Erf[x] - Erf[x+h]
for small vaues of $h$ and large values of $x$.