25A faster solution could be obtained by, for each element in aVec, computing 𝔳ςt(mt − ct,ς)  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.