function y = h(x)
  n = 99;
  poly_val = 1+x;
  for j = 1:n
      poly_val = 1+x*poly_val;
  end
  y = poly_val;
  