for i=1:30
  h(i) = 10^(-i);
  y = d_right(1,h(i));
  err(i) = y + sin(1);
end
