function val=g_t1(theta,par)

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

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