function y = BR_MH_fun1(x)
%

y = exp(-x.^2).*(2+sin(x*5) + sin(x*2));
