function val=g_normal(theta,par)

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

val=-.5*lsig2-.5*(theta-mu).^2./exp(lsig2);

