% Function computes the conditional expectation of 'a'

function f = afa(a)
global mu sig
f = a.*lognpdf(a,mu,sig);