function posterior = prob7(y,xs)
ys = repmat(y,1,length(xs));
posterior = exp(-0.5*(ys-xs).^2 - abs(xs));
