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