¹⁴A faster solution could be obtained by, for each element in $aVec$ , computing $ς v⇁(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.