demo_break_adaptive_simpsons.m

function y = demo_break_adaptive_simpsons(x)
% we are y = 2*x except for set of measure zero.
% But Adaptive Simpson will think we are y = 10.

  set = [-1 -1/2 0 1/2 1];

  % Lie if we are close to the nodes.
  for j = 1:length(set)
    if abs(x - set(j)) < 1e-15
      y = 10;
      return
    end
  end

  % otherwise, tell the truth :)
  y = 2*x;