25A faster solution could be obtained by, for each element in aVec, computing 
of a grid of
values of 
, and then using an approximating interpolating function (rather than the full expectation) in the
FindRoot command. The associated speed improvement is fairly modest, however, so this route was not
pursued.