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;