function fu=fun(th)
fu=(exp(-(th-pi).^2)-3*exp(-2*(th-3*pi/2).^2));
