)
h=bar(xr,Px); set(h,<span class="string">'facecolor'</span>,[1 0 0])
axis([min(xr) max(xr) 0 max(Px)*1.05])
title(<span class="string">'expected (marginal) posterior for (elements of) x'</span>,<span class="string">'fontsize'</span>,14)
<span class="comment">% do a histogram closeup</span>
subplot(212)
h=bar(xr,Px); grid <span class="string">on</span>; set(h,<span class="string">'facecolor'</span>,[1 0 0])
axis([min(xr) max(xr) 0 2*max(Px(1:fix(length(xr)*.49)))])

Pw = zeros(1,length(xr));
<span class="k