function x = gaussienne_perturbee(q);

x = randn(1000,1);
bid = rand(1000,1);
x = x + 5*rand(1000,1).*(bid<q);
