f = @(x) exp(-x.^2/2);
S=[-1 0 1]; % I think these are required to be in the support of the density
x=IA2RMS(f,S,1e4,0);