26Alternatively, the rootfinding operation would be 0 = ˆ𝔳ς(at,ςt)
    t  , where the interpolation function of 𝔳ς(at,ςt)
 t  is used instead. However, the results obtained (especially ˆςt(at)  ) are much less satisfactory.