n = 1000;
delta = 0.1;
x = [0:n-1]'/(n-1);
F = exp (- 100*(x - 1/5).^2) + exp (-  500*(x - 2/5).^2) + ...
    exp (-2500*(x - 3/5).^2) + exp (-12500*(x - 4/5).^2);
randn ('seed', 0);
f = F + delta * randn (size (x));

