bayes2
[x,y]=bg_gauss(0.1,0.1);
sample=sample_rho(x,y,30);
a=linspace(0.01,0.5,60);
b=linspace(0.01,0.5,60);
lv=log_bayesv(a,b,30,sample);
[c,h]=contourf(a,b,lv,20);
clabel (c,h);
colorbar ();

sample300=sample_rho(x,y,300);
plot(x,sample,"*", x, y*30);
plot(x,sample300,"*", x, y*300);

##note review event definition and generation in sample to avoid request a,b<<1


