function posterior = expectedprob7(y,xs)
% "mean", or expected posterior
ys = repmat(y,1,length(xs));
scalefactor = sum(exp(-0.5*(ys-xs).^2 - abs(xs)));
posterior = sum(xs/scalefactor.*exp(-0.5*(ys-xs).^2 - abs(xs)));
