% Function computes the conditional expectation of 'a'
function outvalue = quad_afa4(vk,vf,sd)
outvalue=quad(@(a)a.*exp((-1)*a/sd)./(sd*(1+exp((-1)*a/sd)).^2.),vk,vf);
end