function val=h_nigm(par,hpar)

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

m1=hpar(1); v1=hpar(2); a=hpar(3); b=hpar(4);

val=-.5*(mu-m1).^2/v1-a*lsig2-b*exp(-lsig2);

