f=mexMultiNomBM
f.PXX=[0.8 0.2;
    0.3 0.7];
f.PYX=[0.4 0.6;
    0.2 0.8];
f.p0=[0.5;0.5];
 
f=f.validate
 

f.apost_pdf.mean

for t=1:10
   f=f.bayes(mod(t,2)+1);
   f.apost_pdf.mean
end