function x = exp_arg (j, md, e, K, sj, oj)

  sj = 1./(sqrt (2)*sj);
  ls = log (sj);
  xj = j.*sj + md.*sj + e*K.*sj;
  x = ls - (xj - oj.*sj).^2;

endfunction;

