h = [.3 .1 .03 .01 10^(-5) 10^(-8) 10^(-10) ];
[one, n] = size(h);
x = 1;
for i=1:n
  hh = h(i);
  fp_hat = ( ( x+hh )*exp( x+hh ) - x*exp(x) ) /hh;
  fp     = exp(x) + x*exp(x);
  e      = fp_hat - fp;
  printLine = sprintf( ' %e %e  %e',hh, fp_hat, e);
  disp( printLine)
  clear printline
end
