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