function val=h_normal(par,hpar)

mu=par(:,1); lsig2=par(:,2);

m1=hpar(1); v1=hpar(2); m2=hpar(3); v2=hpar(4);

val=-.5*(mu-m1).^2/v1-.5*(lsig2-m2).^2/v2;

