function [ h_list, fd_list ] = fd_table ( f, x, h, alpha, n )

  p = 0 : n - 1;
  h_list = h * alpha .^ p;
  fd_list = fd ( f, x, h_list );

  return
end


